惊了，填色图向经度0纬度0拉扯

两米阳光

※剑桥¤ 发表于 2014-7-22 22:19
你好，非常感谢你的回复。我从文件里面也读取经纬度数据，经纬度数据中没有纬度是0和经度是0的点。那我的 ...

哦 那我不太清楚了，我当时画H2A数据的时候也出现了你这个问题，就是补上了缺省，就正常出图了，我这里没有你这么新的数据，最近的是4月的，我刚看了看还是有（0,0）点的，不知道为啥7月份的没了

longlivehj

两米阳光 发表于 2014-7-23 08:25
哦 那我不太清楚了，我当时画H2A数据的时候也出现了你这个问题，就是补上了缺省，就正常出图了，我这里没 ...

嗨，请教下，H2A是什么数据啊？没有用过，想了解一下。

两米阳光

longlivehj 发表于 2014-7-23 08:33
嗨，请教下，H2A是什么数据啊？没有用过，想了解一下。

海洋2号卫星数据 我国第一颗海洋动力环境卫星~

※剑桥¤

longlivehj 发表于 2014-7-22 23:34
你把ncl_filedump或者ncdump的结果贴上来，问题一定可以解决！

我帖一下ncdump -c 内容如下。
netcdf H2A_SM2B20140714_14020 {
dimensions:
        lat = 1624 ;
        lon = 76 ;
        dim2 = 4 ;
variables:
        short max_likelihood_est(lat, lon, dim2) ;
                max_likelihood_est:unit = "n/a                     " ;
                max_likelihood_est:scale_factor = 0.001f ;
                max_likelihood_est:scale_factor_err = 0.f ;
                max_likelihood_est:add_offset = 0.f ;
                max_likelihood_est:add_offset_err = 0.f ;
                max_likelihood_est:valid_range = -3000.f ;
        short model_dir(lat, lon) ;
                model_dir:unit = "deg                     " ;
                model_dir:scale_factor = 0.1f ;
                model_dir:scale_factor_err = 0.f ;
                model_dir:add_offset = 0.f ;
                model_dir:add_offset_err = 0.f ;
                model_dir:valid_range = 0.f ;
        short model_speed(lat, lon) ;
                model_speed:unit = "m/s                     " ;
                model_speed:scale_factor = 0.01f ;
                model_speed:scale_factor_err = 0.f ;
                model_speed:add_offset = 0.f ;
                model_speed:add_offset_err = 0.f ;
                model_speed:valid_range = 0.f ;
        short num_ambigs(lat, lon) ;
                num_ambigs:unit = "count                   " ;
                num_ambigs:scale_factor = 1.f ;
                num_ambigs:scale_factor_err = 0.f ;
                num_ambigs:add_offset = 0.f ;
                num_ambigs:add_offset_err = 0.f ;
                num_ambigs:valid_range = 0.f ;
        float time_day(lat) ;
        float time_sec(lat) ;
        short wind_dir(lat, lon, dim2) ;
                wind_dir:unit = "deg                     " ;
                wind_dir:scale_factor = 0.1f ;
                wind_dir:scale_factor_err = 0.f ;
                wind_dir:add_offset = 0.f ;
                wind_dir:add_offset_err = 0.f ;
                wind_dir:valid_range = 0.f ;
        short wind_dir_err(lat, lon, dim2) ;
                wind_dir_err:unit = "deg                     " ;
                wind_dir_err:scale_factor = 0.1f ;
                wind_dir_err:scale_factor_err = 0.f ;
                wind_dir_err:add_offset = 0.f ;
                wind_dir_err:add_offset_err = 0.f ;
                wind_dir_err:valid_range = 0.f ;
        short wind_dir_selection(lat, lon) ;
                wind_dir_selection:unit = "deg                     " ;
                wind_dir_selection:scale_factor = 0.1f ;
                wind_dir_selection:scale_factor_err = 0.f ;
                wind_dir_selection:add_offset = 0.f ;
                wind_dir_selection:add_offset_err = 0.f ;
                wind_dir_selection:valid_range = 0.f ;
        short wind_speed(lat, lon, dim2) ;
                wind_speed:unit = "m/s                     " ;
                wind_speed:scale_factor = 0.01f ;
                wind_speed:scale_factor_err = 0.f ;
                wind_speed:add_offset = 0.f ;
                wind_speed:add_offset_err = 0.f ;
                wind_speed:valid_range = 0.f ;
        short wind_speed_selection(lat, lon) ;
                wind_speed_selection:unit = "m/s                     " ;
                wind_speed_selection:scale_factor = 0.01f ;
                wind_speed_selection:scale_factor_err = 0.f ;
                wind_speed_selection:add_offset = 0.f ;
                wind_speed_selection:add_offset_err = 0.f ;
                wind_speed_selection:valid_range = 0.f ;
        float wvc_lat(lat, lon) ;
                wvc_lat:unit = "deg                     " ;
                wvc_lat:scale_factor = 1.f ;
                wvc_lat:scale_factor_err = 0.f ;
                wvc_lat:add_offset = 0.f ;
                wvc_lat:add_offset_err = 0.f ;
                wvc_lat:valid_range = -90.f ;
        float wvc_lon(lat, lon) ;
                wvc_lon:unit = "deg                     " ;
                wvc_lon:scale_factor = 1.f ;
                wvc_lon:scale_factor_err = 0.f ;
                wvc_lon:add_offset = 0.f ;
                wvc_lon:add_offset_err = 0.f ;
                wvc_lon:valid_range = 0.f ;
        short wvc_quality_flag(lat, lon) ;
                wvc_quality_flag:unit = "na                      " ;
                wvc_quality_flag:scale_factor = 1.f ;
                wvc_quality_flag:scale_factor_err = 0.f ;
                wvc_quality_flag:add_offset = 0.f ;
                wvc_quality_flag:add_offset_err = 0.f ;
                wvc_quality_flag:valid_range = 0.f ;
        short wvc_selection(lat, lon) ;
                wvc_selection:unit = "count                   " ;
                wvc_selection:scale_factor = 1.f ;
                wvc_selection:scale_factor_err = 0.f ;
                wvc_selection:add_offset = 0.f ;
                wvc_selection:add_offset_err = 0.f ;
                wvc_selection:valid_range = 0.f ;

// global attributes:
                :longName = "HY-2A Level 2B Ocean Wind Vectors in 25 km Swath Grid                                           " ;
                :ShortName = "HY-2A SCAT-L2B          " ;
                :producer_agency = "CNSA                    " ;
                :producer_institution = "NSOAS                   " ;
                :InstrumentShorName = "HY-2A SCAT              " ;
                :platformLongName = "Haiyang 2A Ocean Observing Satellite                                                            " ;
                :platformShortName = "HY-2A                   " ;
                :platformType = "spacecraft                      " ;
                :project_id = "HY-2A SCAT                      " ;
                :data_format_type = "NCSA HDF                        " ;
                :granulePointer = "H2A_SM2B20140714_14020.h5                                                                                                                                                                                                                                                       " ;
                :inputPointer = "H2A_SM2A20140714_14020.h5                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       " ;
                :ancillary_data_descriptors = "land_mask \n",
                        " ice_mask \n",
                        " NWP wind                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 " ;
                :sis_id = "sis1.0                                  " ;
                :build_id = "V1.0                                    " ;
                :HDF_version_id = "HDF5-1.6.5/2011-8-15                    " ;
                :parameterName = "sea surface wind                        " ;
                :operationMode = "wind observation                        " ;
                :rangeBeginningDate = "20140714                        " ;
                :rangeEndingDate = "20140714                        " ;
                :rangeBeginningTime = "20140714T00:10:20.729           " ;
                :rangeEndingTime = "20140714T01:52:00.978           " ;
                :productionDateTime = "Mon Jul 14 10:11:38 2014                " ;
                :ephemeris_type = "GPS                             " ;
                :l2b_algorithm_descrptor = "Use look up table for ku band GMF/n  Applies median filter technique for ambiguity removal. Ambiguity median filter is based on wind vectors over a 7 by 7 wind vector window                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   " ;
                :median_filter_method = "circular median filter                  " ;
                :l2b_actual_wvc_rows = "1624                    " ;
                :l2b_expected_wvc_rows = "1624                    " ;
                :sigma0_granularity = "whole egg                       " ;
data:
}

※剑桥¤

最终解决情况：
因为经纬度数据的缺省值是0，所以需要将缺省值替换成其他值，最后通过语句解决：
lat2d = where(lat2d.eq.0,-9999,lat2d)
谢谢@两米阳光。

※剑桥¤

也谢谢@longlivehj

风荷1

第一次见这样的问题

mimimomo2010

滿有趣的問題...有時0跟無義值是不同意思的....
尤其是雨量...
where的用法也可以學一下...
把矩陣中的一些值挑掉再做運算....