[FWTools] gdalwarp newbie
Dave Zuin
d.zuin at wni.com
Wed Apr 23 16:55:52 EDT 2008
Hi All,
I am trying to use gdalwarp to project a polar stereographic image to a
mercator prjection (for google Maps).
i know the source proj.4 string and the destination proj.4 string so i
tried:
gdalwarp -s_srs "+proj=stere +lat_0=90 +lon_0=0 +lat_ts=60 +a=6378.14
+b=6356.75 +x_0=0 y_0=0" -t_srs "+proj=merc +a=6378137 +b=6378137
+lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null
+no_def" ..\borders.png output.tif
which results in this error message:
ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates for ..\borders.png.
There is no affine transformation and no GCPs.
What am i missing?
Kind Regards Dave
Anyone interested in the png, you can find it here:
http://www.xs4all.nl/~zuin/borders.png
More geo info about the file (from gdalinfo from the source hdf5)
Driver: HDF5/Hierarchical Data Format Release 5
Files: ..\RAD_NL25_PCP_NA_200801070600.h5
Size is 512, 512
Coordinate System is `'
Metadata:
map_projection:projection_indication=Y
map_projection:projection_name=STEREOGRAPHIC
map_projection:projection_proj4_params=+proj=stere +lat_0=90 +lon_0=0
+lat_ts=60 +a=6378.14 +b=6356.75 +x_0=0 y_0=0
geographic:geo_par_pixel=X,Y
geographic:geo_dim_pixel=KM,KM
geographic:geo_pixel_def=LU
geographic:geo_number_columns=700
geographic:geo_number_rows=765
geographic:geo_pixel_size_x=1.000001
geographic:geo_pixel_size_y=-1.000005
geographic:geo_column_offset=0.000000
geographic:geo_row_offset=3649.980225
geographic:geo_product_corners=0.000000 49.362061 0.000000 55.973602
10.856429 55.388973 9.009280 48.895302
calibration:calibration_flag=Y
calibration:calibration_formulas=GEO = 0.500000 * PV + -31.500000
calibration:calibration_missing_data=0
calibration:calibration_out_of_image=255
statistics:stat_min_value=-31.000000
statistics:stat_max_value=35.000000
image1:image_product_name=RAD_NL25_PCP_H1.5_NA
image1:image_geo_parameter=REFLECTIVITY_[DBZ]
image1:image_size=535500
image1:image_bytes_per_pixel=1
overview:product_group_name=RAD_NL25_PCP_NA
overview:products_missing=NA
overview:hdftag_version_number=3.5
overview:number_image_groups=1
overview:number_satellite_groups=0
overview:number_radar_groups=2
overview:number_station_groups=0
overview:product_datetime_start=07-JAN-2008;06:00:00.000
overview:product_datetime_end=07-JAN-2008;06:00:00.000
radar1:radar_name=De_Bilt
radar1:radar_location=5.178340 52.101681
radar2:radar_name=Den_Helder
radar2:radar_location=4.799970 52.953339
Subdatasets:
SUBDATASET_0_NAME=HDF5:"..\RAD_NL25_PCP_NA_200801070600.h5"://image1/image_d
ata
SUBDATASET_0_DESC=[765x700] //image1/image_data (8-bit unsigned character)
SUBDATASET_1_NAME=HDF5:"..\RAD_NL25_PCP_NA_200801070600.h5"://visualisation1
/color_palette
SUBDATASET_1_DESC=[256x3] //visualisation1/color_palette (8-bit unsigned
character)
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 512.0)
Upper Right ( 512.0, 0.0)
Lower Right ( 512.0, 512.0)
Center ( 256.0, 256.0)
More information about the FWTools
mailing list