[Thread Prev][Thread Next][Index]
Re: [ferret_users] regriding data from (lon lat) axis to (x_page y_page) axis (stereographic projection)
- To: Cécile Agosta <cecile.agosta@xxxxxxxxx>
- Subject: Re: [ferret_users] regriding data from (lon lat) axis to (x_page y_page) axis (stereographic projection)
- From: Peter Szabo <szabpet83@xxxxxxxxx>
- Date: Wed, 14 May 2008 14:51:37 +0200
- Cc: ferret_users@xxxxxxxx
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:cc:in-reply-to:mime-version:content-type:references; bh=Ao3D2ZuZvrKK95OfM18Iit7xBOaQt+pOsjrHCmCbSbc=; b=keKS1Idw71espQgz4W+Up+A5DpvcsRYzy6/Zld0k3usFrtrv78Tmk//cRI+msp8QC4JXVfKEvchpkDiKk1MU6YGL2nvlbwfsfeAZmx6oLp6M1ZTplbb390ssFj5YAM9olIlOUpBB6ztKTS/rtTW5femcO8M+wGEm8aDJQx4OJwY=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:cc:in-reply-to:mime-version:content-type:references; b=aE6fy1RXhqn2XARX7A167coOLZfDeWQLRTc58TRB+wMXmWimeICp+foip9AmtfnDwiPaDz3BuXhhr0VLqJR+EPwlivpKPctmSebfIAKZtimqVMNQgn14DKIar5JtMDcxLTJX9lSfzKRgbgakCR1lKHIb5rhVKZx9/9Izjtap1wc=
- In-reply-to: <c85e9c350805130923o4c42f0aem66f5d84bbe0d19db@xxxxxxxxxxxxxx>
- References: <c85e9c350805130923o4c42f0aem66f5d84bbe0d19db@xxxxxxxxxxxxxx>
- Sender: owner-ferret_users@xxxxxxxx
Dear Cécile,
As i understand you question you need to have your (rectilinear) data set in a curvilinear grid.
You can have it easily by the RECT_TO_CURV function. You only need two curvilinear grid variables (2D fields of lons and lats describing the curvilinear grid) from an other, curvilinear file or from yours:
go mp_stereographic_south
set grid your_old
!now you have variables like `mp_standard_parallel` and `mp_central_meridian` and x_page and y_page
!the `mp_standard_parallel and `mp_central_meridian` are constants (the middle of your map), the x_page and y_page are variables (they are just differences from the standard_parallel and central_meridian). so somehow you need to add them or sub them to get the correct curvilinear coordinates. i did not have time for further analyse.
let x_st=x_page+`mp_standard_parallel`
let y_st=y_page+`mp_central_meridian`
after this you can use:
let your_new=rect_to_curv(your_old,x_st,y_st,1)
save/file="..." your_new
hope this helps a bit,
Peter
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement