From: Alan W. I. <ir...@be...> - 2005-03-25 16:10:21
|
On 2005-03-25 13:21+0100 Arjen Markus wrote: > Curtis Cooper wrote: >> >>> Hi, >>> >>> I'm using PLplot to produce pictures of longitude, latitude temperature >>> maps for my planetary atmospheres research. I first apply an orthographic >>> projection to the lon, lat coordinates and use the grid so computed as the >>> argument into plshades along with pltr2: >>> >>> plshades(temperature, clevels, fill_width, 0, pltr2, X, Y, 2) >>> >>> Temperature here is known at all longitudes and latitudes for which I have >>> computed the X,Y projection matrices, i.e., T(lon, lat). >>> >>> This is just like example 16. Now, the routine plots the map projection >>> beautifully. However, for further computations, I would really like to be >>> able to extract the projected data, i.e., T in the X,Y coordinate system. >>> What would be involved in modifying the PLplot source code to allow me to >>> do that? >>> > > Nothing: > > The function pltr2() is doing the actual coordinate transformation. You > can > call it yourself like this: > > for (i = 0; i < n; i++) { > pltr2 (x[i], y[i], &tx, &ty, pltr_data); > x[i] = tx; > y[i] = ty; > } > > where pltr_data is the grid matrix. It may not be much further help beyond what Arjen has already said, but sparse documentation of pltr2 actually exists at http://plplot.sourceforge.net/docbook-manual/plplot-html-5.3.1/pltr2.html . Alan __________________________ Alan W. Irwin email: ir...@be... phone: 250-727-2902 Astronomical research affiliation with Department of Physics and Astronomy, University of Victoria (astrowww.phys.uvic.ca). Programming affiliations with the FreeEOS equation-of-state implementation for stellar interiors (freeeos.sf.net); PLplot scientific plotting software package (plplot.org); the Yorick front-end to PLplot (yplot.sf.net); the Loads of Linux Links project (loll.sf.net); and the Linux Brochure Project (lbproject.sf.net). __________________________ Linux-powered Science __________________________ |