2009-08-13 83 views
0

我使用的是unidata的NetCDF 4 java庫來讀取GRIB2文件,現在我試圖找到一種方法來獲取給定緯度/經度。我有一個北半球的.flt文件,其中包含所有這些信息,我知道如何閱讀它,但我不知道如何正確地進行地圖投影。我能夠使用 LambertConformal(originLat,originLon,stdParallel1,stdParallel2)來進行投影,但我不知道如何在網格中使用投影的x,y值。Lambert共形圓錐映射投影得到x,y值中的一個網格

網格值爲Nx = 1073和Ny = 689,我得到的預測值爲lat = 33.92,lon = -84.33爲x = 995 y = 1035。我知道預測值是正確的,但我不知道如何與他們一起閱讀網格。

也許有人可以幫助我理解下面的值。

:grid_mapping_name = "lambert_conformal_conic"; 
:standard_parallel = 25.0; // double 
:longitude_of_central_meridian = 265.0; // double 
:latitude_of_projection_origin = 25.0; // double 
:earth_shape = "Earth spherical with radius specified by producer"; 
:spherical_earth_radius_meters = 6371200.0; // double 
:GRIB_param_Dx = 5079.4062; // double 
:GRIB_param_Dy = 5079.4062; // double 
:GRIB_param_GDSkey = 1158554943; // int 
:GRIB_param_La1 = 20.192; // double 
:GRIB_param_LaD = 25.0; // double 
:GRIB_param_Latin1 = 25.0; // double 
:GRIB_param_Latin2 = 25.0; // double 
:GRIB_param_Lo1 = 238.446; // double 
:GRIB_param_LoV = 265.0; // double 
:GRIB_param_NpProj = "true"; 
:GRIB_param_Nx = 1073; // int 
:GRIB_param_Ny = 689; // int 
:GRIB_param_ProjFlag = 0; // int 
:GRIB_param_Quasi = "false"; 
:GRIB_param_ResCompFlag = 0; // int 
:GRIB_param_SpLat = -90.0; // double 
:GRIB_param_SpLon = 0.0; // double 
:GRIB_param_Winds = "True"; 
:GRIB_param_grid_name = "Lambert Conformal"; 
:GRIB_param_grid_radius_spherical_earth = 6371200.0; // double 
:GRIB_param_grid_shape = "Earth spherical with radius specified by producer"; 
:GRIB_param_grid_shape_code = 1; // int 
:GRIB_param_grid_type = 30; // int 
:GRIB_param_grid_units = "m"; 

回答

1

您是否嘗試過張貼Unidata的網站上?他們可能會更好地回答這類問題。過去我得到了他們的很好迴應。 -Fred

+0

這似乎是一個合理的建議!我得看看他們的郵件列表,謝謝。 – ravun 2009-08-22 02:01:55