Re: hdf plotting error

Jinqing,

In short, you can not just define a region and plot the satellite image
there.

The file is like a pixel image, with only information about its center, and
the height
of the satellite, so there is no simple way of plot a region as you defined=
.


If you really want to plot that region, you have to find a way to
project/transfer the
satellite data to lat,lon grid, and then you can view the region specified.

Regards,

Wei


On Tue, Jun 10, 2014 at 4:54 PM, 刘金卿 <liujq_0912@12=
6.com> wrote:

> Thanks for your patience. I got the right image after adding data(::-1,=
:
> ).
> Can I define a specific domain ? The following sentence do not work for
> me.
>
> res@mpMinLonF = 120
> res@mpMaxLonF = 145
> res@mpMinLatF = 30
> res@mpMaxLatF = 60
>
>
> On Jun 11, 2014, at 1:07, Wei Huang <huangwei@ucar.edu> wrote:
>
> Jinqing,
>
> This is a satellite image, we have to tune to make things look right.
> Attached are a script, and a plot generated.
> You probably need to consult with a satellite image expert to see
> what those parameters mean.
>
> Regards,
>
> Wei
>
>
>
> On Tue, Jun 10, 2014 at 8:50 AM, 刘金卿 <liujq_0912@=
126.com> wrote:
>
>> Hi, Wei,
>> Thanks for your reply, but the picture still seems wrong compared with
>> the image from website quick view (attachment).
>> And I add
>> res@mpDataSetName = "Earth..4"
>> res@mpDataBaseVersion = "MediumRes"
>> res@mpOutlineSpecifiers = (/"China:states","India","Taiwan�=
�/)
>> in the script, but it dose not work. Is there anyway to turn on the
>> boundaries for this data?
>>
>>
>>
>>
>> On Jun 10, 2014, at 0:41, Wei Huang <huangwei@ucar.edu> wrote:
>>
>> Jinqing,
>>
>> My may try this script:
>>
>> ;************************************************
>> fn = "FY2C_TBB_IR1_NOM_20070816_1200.hdf.h5"
>> f = addfile(fn, "r")
>> vNam = getfilevarnames(f) ; all variables names on
>> file
>> print(vNam)
>> ;************************************************
>> ;pltType = "x11"
>> pltType = "png"
>> pltName = "T"
>> ;************************************************
>> ; plotting parameters
>> ;************************************************
>> wks = gsn_open_wks (pltType,pltName) ; open workstation
>> ;gsn_define_colormap(wks,"amwg") ; choose colormap
>> ;gsn_define_colormap(wks,"gsdtol")
>> ;****************************************************
>> res = True ; plot mods desired
>> res@gsnMaximize = True ; make ps/eps/pdf large
>> res@gsnSpreadColors = True
>>
>> res@cnFillOn = True
>> ;res@cnLinesOn = True
>> ;res@cnLineLabelsOn = True
>> res@cnFillMode = "RasterFill" ; faster
>> ;res@cnRasterSmoothingOn = True
>> ;****************************************************
>> x=f->FY2C_TBB_Hourly_Product(:,:)
>>
>> x = where(x.lt.150.0, -999.0, x)
>> x@_FillValue = -999.0
>>
>> printVarSummary(x)
>> plot = gsn_csm_contour(wks, x, res)
>>
>>
>> Wei
>>
>>
>>
>> On Mon, Jun 9, 2014 at 4:56 AM, 刘金卿 <liujq_0912@=
126.com> wrote:
>>
>>> Hi, I try to plot a hdf file, but got the following warnings:
>>> warning:_NhlCreateSplineCoordApprox: Attempt to create spline
>>> approximation for X axis failed: consider adjusting trXTensionF value
>>> warning:IrTransInitialize: error creating spline approximation for
>>> trXCoordPoints; defaulting to linear
>>>
>>> The picture seems wrong too. The data has been transferred on the ftp,
>>> and the script is in attachment. Any suggestions will be appreciate.
>>> ftp> put FY2C_TBB_IR1_NOM_20070816_1200.hdf
>>>
>>> -Liu
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>> [image: 提示图标] 邮件�=
�有附件预览链接，=
若您转发或回复此�=
��件时不希望对方�=
�览附件，建议您手=
动删除链接。
>> 共有 *1* 个附件
>> T.png(458K)极速下载
>> <http://preview.mail.126.com/xdownload?filename=T.png&mid=1tbi1xqlc0=
gY5CQ6YwAAsS&part=3&sign=6631ea7969a28a8329adddcea854b6c2&time=140233=
6372&uid=liujq_0912%40126.com>
>> 在线预览
>> <http://preview.mail.126.com/preview?mid=1tbi1xqlc0gY5CQ6YwAAsS&part=
=3&sign=6631ea7969a28a8329adddcea854b6c2&time=1402336372&uid=liujq_=
0912%40126.com>
>> <T.png>
>>
>>
>>
>>
> [image: 提示图标] 邮件带=
有附件预览链接，�=
��您转发或回复此�=
�件时不希望对方预=
览附件，建议您手�=
��删除链接。
> 共有 *2* 个附件
> T.png(450K) 极速下载
> <http://preview.mail.126.com/xdownload?filename=T.png&mid=1tbiIQSmc04=
gDUzmmQAAsg&part=3&sign=a9406ec56b310ad66bce10f7ae166d88&time=1402437=
915&uid=liujq_0912%40126.com>
> 在线预览
> <http://preview.mail.126.com/preview?mid=1tbiIQSmc04gDUzmmQAAsg&part==
3&sign=a9406ec56b310ad66bce10f7ae166d88&time=1402437915&uid=liujq_091=
2%40126.com>
> fy2c.ncl(4K)极速下载
> <http://preview.mail.126.com/xdownload?filename=fy2c.ncl&mid=1tbiIQSm=
c04gDUzmmQAAsg&part=4&sign=a9406ec56b310ad66bce10f7ae166d88&time=1402=
437915&uid=liujq_0912%40126.com>
> <T.png><fy2c.ncl>
>
>
>

Received on Wed Jun 11 02:37:16 2014