- Messages
- 295
- Country

I have FWTools2.4.7 in Windows 10. I can gdalwarp my 10 VT jp2 files with no problem.
C:\FWTools2.4.7\bin\gdalwarp.exe -t_srs "+proj=latlong +datum=WGS84" -of GTiff -co "INTERLEAVE=PIXEL" -r cubic %1.jp2 %1WARP.TIF
C:\FWTools2.4.7\bin\listgeo.exe %1WARP.tif > %1WARP.gtf
I reloaded the Geo data to get a GeoTiff file with this:
C:\FWTools2.4.7\bin\geotifcp.exe -g %1warp.gtf %1warp.tif %1GeoTiff.tif
Output of my resample:
Resample Info File: 'VT.INF'
[Source] Type = 'Multisource'
[Source] NumberOfSources = 10
[Source1] Type = 'GeoTiff'
[Source1] SourceDir = '.'
[Source1] SourceFile = 'VT_44642197_20180422Geotiff.tif'
[Source1] Layer = 'Imagery'
[Source1] NullValue = '0,0,0'
[Source1] Variation = 'Day'
Pixel width and height must be greater than zero.
Unable to open data source 'C:\FWTools2.4.7\bin\VT_44642197_20180422Geotiff.tif'.
Failed to create data source.
The gdalinfo of the JPEG2000 file "looks" OK to me:
C:\FWTools2.4.7\bin>gdalinfo VT_44642197_20180422.jp2
Driver: JP2ECW/ERMapper JPEG2000
Files: VT_44642197_20180422.jp2
Size is 6000, 6000
Coordinate System is:
PROJCS["NAD_1983_2011_StatePlane_Vermont_FIPS_4400",
GEOGCS["GCS_NAD_1983_2011",
DATUM["NAD_1983_2011",
SPHEROID["GRS_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",42.5],
PARAMETER["central_meridian",-72.5],
PARAMETER["scale_factor",0.9999642857142858],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (446000.000000000000000,220200.000000000000000)
Pixel Size = (0.150000000000000,-0.150000000000000)
Corner Coordinates:
Upper Left ( 446000.000, 220200.000) ( 73d10'43.63"W, 44d28'48.04"N)
Lower Left ( 446000.000, 219300.000) ( 73d10'43.29"W, 44d28'18.88"N)
Upper Right ( 446900.000, 220200.000) ( 73d10'2.90"W, 44d28'48.28"N)
Lower Right ( 446900.000, 219300.000) ( 73d10'2.57"W, 44d28'19.12"N)
Center ( 446450.000, 219750.000) ( 73d10'23.10"W, 44d28'33.58"N)
Band 1 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 2 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 3 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 4 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
My GTF Files "look" ok to me.
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
-73.178784856304 44.4800773404146 0
ModelPixelScaleTag (1,3):
1.64394339422765e-0061.64394339422765e-0060
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeGeographic
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GeographicTypeGeoKey (Short,1): GCS_WGS_84
GeogCitationGeoKey (Ascii,7): "WGS 84"
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
GeogSemiMajorAxisGeoKey (Double,1): 6378137
GeogInvFlatteningGeoKey (Double,1): 298.257223563
End_Of_Keys.
End_Of_Geotiff.
The gdalinfo of the GeoTiff file "looks" OK:
Example:
GCS: 4326/WGS 84
Datum: 6326/World Geodetic System 1984
Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Corner Coordinates:
Upper Left ( 73d10'43.63"W, 44d28'48.28"N)
Lower Left ( 73d10'43.63"W, 44d28'18.88"N)
Upper Right ( 73d10' 2.57"W, 44d28'48.28"N)
Lower Right ( 73d10' 2.57"W, 44d28'18.88"N)
Center ( 73d10'23.10"W, 44d28'33.58"N)
GCP Projection =
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"]]
GCP[ 0]: Id=1, Info=
(0,0) -> (-73.178784856304,44.4800773404146,0)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 4967.0)
Upper Right ( 6937.0, 0.0)
Lower Right ( 6937.0, 4967.0)
Center ( 3468.5, 2483.5)
Band 1 Block=6937x1 Type=Byte, ColorInterp=Red
Mask Flags: PER_DATASET ALPHA
Band 2 Block=6937x1 Type=Byte, ColorInterp=Green
Mask Flags: PER_DATASET ALPHA
Band 3 Block=6937x1 Type=Byte, ColorInterp=Blue
Mask Flags: PER_DATASET ALPHA
Band 4 Block=6937x1 Type=Byte, ColorInterp=Alpha
I looked at the posts about this problem but saw nothing that would help me.
I don't have the app to shrink this image in Irfanview. I attached the zipped inf file.
C:\FWTools2.4.7\bin\gdalwarp.exe -t_srs "+proj=latlong +datum=WGS84" -of GTiff -co "INTERLEAVE=PIXEL" -r cubic %1.jp2 %1WARP.TIF
C:\FWTools2.4.7\bin\listgeo.exe %1WARP.tif > %1WARP.gtf
I reloaded the Geo data to get a GeoTiff file with this:
C:\FWTools2.4.7\bin\geotifcp.exe -g %1warp.gtf %1warp.tif %1GeoTiff.tif
Output of my resample:
Resample Info File: 'VT.INF'
[Source] Type = 'Multisource'
[Source] NumberOfSources = 10
[Source1] Type = 'GeoTiff'
[Source1] SourceDir = '.'
[Source1] SourceFile = 'VT_44642197_20180422Geotiff.tif'
[Source1] Layer = 'Imagery'
[Source1] NullValue = '0,0,0'
[Source1] Variation = 'Day'
Pixel width and height must be greater than zero.
Unable to open data source 'C:\FWTools2.4.7\bin\VT_44642197_20180422Geotiff.tif'.
Failed to create data source.
The gdalinfo of the JPEG2000 file "looks" OK to me:
C:\FWTools2.4.7\bin>gdalinfo VT_44642197_20180422.jp2
Driver: JP2ECW/ERMapper JPEG2000
Files: VT_44642197_20180422.jp2
Size is 6000, 6000
Coordinate System is:
PROJCS["NAD_1983_2011_StatePlane_Vermont_FIPS_4400",
GEOGCS["GCS_NAD_1983_2011",
DATUM["NAD_1983_2011",
SPHEROID["GRS_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",42.5],
PARAMETER["central_meridian",-72.5],
PARAMETER["scale_factor",0.9999642857142858],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (446000.000000000000000,220200.000000000000000)
Pixel Size = (0.150000000000000,-0.150000000000000)
Corner Coordinates:
Upper Left ( 446000.000, 220200.000) ( 73d10'43.63"W, 44d28'48.04"N)
Lower Left ( 446000.000, 219300.000) ( 73d10'43.29"W, 44d28'18.88"N)
Upper Right ( 446900.000, 220200.000) ( 73d10'2.90"W, 44d28'48.28"N)
Lower Right ( 446900.000, 219300.000) ( 73d10'2.57"W, 44d28'19.12"N)
Center ( 446450.000, 219750.000) ( 73d10'23.10"W, 44d28'33.58"N)
Band 1 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 2 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 3 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
Band 4 Block=6000x1 Type=Byte, ColorInterp=Undefined
Overviews: arbitrary
My GTF Files "look" ok to me.
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
-73.178784856304 44.4800773404146 0
ModelPixelScaleTag (1,3):
1.64394339422765e-0061.64394339422765e-0060
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeGeographic
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GeographicTypeGeoKey (Short,1): GCS_WGS_84
GeogCitationGeoKey (Ascii,7): "WGS 84"
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
GeogSemiMajorAxisGeoKey (Double,1): 6378137
GeogInvFlatteningGeoKey (Double,1): 298.257223563
End_Of_Keys.
End_Of_Geotiff.
The gdalinfo of the GeoTiff file "looks" OK:
Example:
GCS: 4326/WGS 84
Datum: 6326/World Geodetic System 1984
Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Corner Coordinates:
Upper Left ( 73d10'43.63"W, 44d28'48.28"N)
Lower Left ( 73d10'43.63"W, 44d28'18.88"N)
Upper Right ( 73d10' 2.57"W, 44d28'48.28"N)
Lower Right ( 73d10' 2.57"W, 44d28'18.88"N)
Center ( 73d10'23.10"W, 44d28'33.58"N)
GCP Projection =
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"]]
GCP[ 0]: Id=1, Info=
(0,0) -> (-73.178784856304,44.4800773404146,0)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 4967.0)
Upper Right ( 6937.0, 0.0)
Lower Right ( 6937.0, 4967.0)
Center ( 3468.5, 2483.5)
Band 1 Block=6937x1 Type=Byte, ColorInterp=Red
Mask Flags: PER_DATASET ALPHA
Band 2 Block=6937x1 Type=Byte, ColorInterp=Green
Mask Flags: PER_DATASET ALPHA
Band 3 Block=6937x1 Type=Byte, ColorInterp=Blue
Mask Flags: PER_DATASET ALPHA
Band 4 Block=6937x1 Type=Byte, ColorInterp=Alpha
I looked at the posts about this problem but saw nothing that would help me.
I don't have the app to shrink this image in Irfanview. I attached the zipped inf file.
Attachments
Last edited by a moderator:

