0

(Matlab で生成された) GeoTIFF ファイルを WorldWind にインポートしようとしていますが、うまくいきません。有用なヒントは大歓迎です。GeoTIFF ファイルは ArcGIS で正常に表示されます (エクスポート時に .tfw ファイルを作成できます) が、WorldWind では次のメッセージが表示されます。

SEVERE: Cannot read raster: C:\Users\Matthias\Desktop\geotiff\fldextent_02- 
Jan-1977(1)renderedno0.tif : gov.nasa.worldwind.formats.tiff.GeotiffImageReader.read(): unable
to decipher image organization
Jul 09, 2013 6:54:33 PM gov.nasa.worldwind.data.CachedDataRaster drawOnTo
SEVERE: C:\Users\Matthias\Desktop\geotiff\fldextent_02-Jan-1977(1)renderedno0.tif : Cannot read
raster: C:\Users\Matthias\Desktop\geotiff\fldextent_02-Jan-1977(1)renderedno0.tif : 
gov.nasa.worldwind.formats.tiff.GeotiffImageReader.read(): unable to decipher image organization
gov.nasa.worldwind.exception.WWRuntimeException: Cannot read raster: C:\Users\Matthias\Desktop
\geotiff\fldextent_02-Jan-1977(1)renderedno0.tif :
gov.nasa.worldwind.formats.tiff.GeotiffImageReader.read(): unable to decipher image organization
   at gov.nasa.worldwind.data.CachedDataRaster.getDataRasters(CachedDataRaster.java:255)
   at gov.nasa.worldwind.data.CachedDataRaster.drawOnTo(CachedDataRaster.java:290)
   at gov.nasa.worldwind.data.TiledRasterProducer.drawDataSources(TiledRasterProducer.java:576) 
[...]

また、FWTools で GeoTIFF ファイルの属性を確認したところ、次の結果が得られました。

C:\Users\Matthias\Desktop\geotiff>gdalinfo fldextent_02-Jan-1977(1)renderedno0.tif
Driver: GTiff/GeoTIFF
Files: fldextent_02-Jan-1977(1)renderedno0.tif
   fldextent_02-Jan-1977(1)renderedno0.tfw
Size is 7200, 7200
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (99.000000000000000,7.000000000000000)
Pixel Size = (0.000833333333333,-0.000833333333333)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  99.0000000,   7.0000000) ( 99d 0'0.00"E,  7d 0'0.00"N)
Lower Left  (  99.0000000,   1.0000000) ( 99d 0'0.00"E,  1d 0'0.00"N)
Upper Right ( 105.0000000,   7.0000000) (105d 0'0.00"E,  7d 0'0.00"N)
Lower Right ( 105.0000000,   1.0000000) (105d 0'0.00"E,  1d 0'0.00"N)
Center      ( 102.0000000,   4.0000000) (102d 0'0.00"E,  4d 0'0.00"N)
Band 1 Block=128x128 Type=Byte, ColorInterp=Gray
    NoData Value=0

.tfw ファイルの内容は次のとおりです。

0.0008333333
0.0000000000
0.0000000000
-0.0008333333
99.0004166667
6.9995833333
4

1 に答える 1

0

私は最終的に問題を発見しました：

重要なことは、Matlab で CLEAN GeoTIFF ファイルを作成することです (透明度のための RGB およびアルファ レイヤー)。ここで、Matlab のガイダンスをいくつか示します。結果の GeoTIFF は、WorldWind に直接インポートできます。

%%% read intensity values Z (2D matrix) - with values of 0 and above 
%%% (we want 0 to be completely transparent in the final geotiff) - 
%%% together with spatialref.GeoRasterReference ss
[Z, ss] = geotiffread('./flddph_1976-01-01.tif');
info_3 = geotiffinfo('./flddph_1976-01-01.tif');

%%% generate indexed image with 0 to 255 (255 equals max. intensity)
indexedimage = gray2ind(Z);
indexedimage = double(indexedimage);

%%% normalize so that everything between 0 and 1
normalizedimg = (indexedimage) / 255;

%%% scaling data and applying colormap
imgscaled = uint8(256*normalizedimg); % scale data

cmp = makeColorMap([1 1 0],[1 0.75 0],[1 0 0],256); 
% 256 element colormap yellow - orange - red 
% (download appropriate function MAKECOLORMAP)
imgrgb = ind2rgb(imgscaled,cmp);

%%% check plot
% subplot(2,1,1)
% imagesc(indexedimage)
% title('indexed image')
% subplot(2,1,2)
% image(img)
% title('rgb image')

%%% generating alpha layer for transparency 
%%% (255 for non-transparent, 0 for transparent)

alpha3 = Z;
alpha3(alpha3>0)=255;
alpha3 = uint8(alpha3);

out2 = cat(3,imgrgb,alpha3);

geotiffwrite('test_rgbhope_flddph.tif',out2,ss);
于 2013-07-11T22:34:03.007 に答える