[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