- Timestamp:
- 07/23/09 11:18:49 (15 years ago)
- Location:
- trunk/SRC/Documentation/idldoc_assistant_output/Interpolation
- Files:
-
- 23 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/clickincell.html
r338 r402 107 107 = 'T', 'W', 'U', 'V' or 'F' 108 108 the type of point 109 110 111 109 that is located in the center of the cell which the click is 110 located. default is T type of cell (with corner defined by F 111 points). 112 112 113 113 … … 199 199 200 200 <h3>Examples</h3><pre> 201 IDL> res = clickincell() 202 Click with the left button to select a cell. Clicking one more 203 time in the same cell remove the cell from the selection. 204 Click on the right button to quit. 205 206 IDL> plt, findgen(jpi,jpj),/nodata,map=[90,0,0],/ortho 207 IDL> print, clickincell(/draw,color=150,/xy) 201 IDL> res = clickincell() 202 203 Click with the left button to select a cell. Clicking one more 204 time in the same cell remove the cell from the selection. 205 Click on the right button to quit. 206 207 IDL> plt, findgen(jpi,jpj),/nodata,map=[90,0,0],/ortho 208 IDL> print, clickincell(/draw,color=150,/xy) 208 209 209 210 </pre><h3>Version history</h3> 210 211 211 212 <h4>Version</h4> 212 $Id: clickincell.pro 3 27 2007-12-13 16:22:35Z pinsard $213 $Id: clickincell.pro 372 2008-08-08 12:31:53Z pinsard $ 213 214 214 215 <h4>History</h4> … … 228 229 229 230 <h4>Uses routines</h4> 230 common.pro231 <a href="..//Obsolete/common.html">common</a> 231 232 232 233 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/compute_fromirr_bilinear_weigaddr.html
r338 r402 235 235 236 236 <h4>Version</h4> 237 $Id: compute_fromirr_bilinear_weigaddr.pro 3 25 2007-12-06 10:04:53Z pinsard $237 $Id: compute_fromirr_bilinear_weigaddr.pro 372 2008-08-08 12:31:53Z pinsard $ 238 238 239 239 <h4>History</h4> 240 240 June 2006: Sebastien Masson (smasson@lodyc.jussieu.fr) 241 241 242 242 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/compute_fromreg_bilinear_weigaddr.html
r338 r402 236 236 237 237 <h4>Version</h4> 238 $Id: compute_fromreg_bilinear_weigaddr.pro 3 27 2007-12-13 16:22:35Z pinsard $238 $Id: compute_fromreg_bilinear_weigaddr.pro 372 2008-08-08 12:31:53Z pinsard $ 239 239 240 240 <h4>History</h4> 241 241 November 2005: Sebastien Masson (smasson@lodyc.jussieu.fr) 242 242 243 243 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/compute_fromreg_imoms3_weigaddr.html
r338 r402 236 236 237 237 <h4>Version</h4> 238 $Id: compute_fromreg_imoms3_weigaddr.pro 3 27 2007-12-13 16:22:35Z pinsard $238 $Id: compute_fromreg_imoms3_weigaddr.pro 372 2008-08-08 12:31:53Z pinsard $ 239 239 240 240 <h4>History</h4> 241 242 241 November 2005: Sebastien Masson (smasson@lodyc.jussieu.fr) 242 March 2006: works for rectangular grids 243 243 244 244 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/cutpar.html
r338 r402 275 275 <h3>Examples</h3><pre> 276 276 277 IDL> x0 = [2,6,2]278 IDL> y0 = [0,2,6]279 IDL> x1 = [3,8,4]280 IDL> y1 = [4,4,6]281 IDL> x2 = [1,6,4]282 IDL> y2 = [5,6,8]283 IDL> x3 = [0,4,2]284 IDL> y3 = [1,4,8]285 IDL> n = 4286 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata287 IDL> for i=0,2 do oplot, [x0[i],x1[i],x2[i],x3[i],x0[i]],[y0[i],y1[i],y2[i],y3[i],y0[i]]288 IDL> res=cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n)289 IDL> for i=0,2 do oplot, [res[0,*,i]], [res[1,*,i]], color = 20+10*i, psym = 1, thick = 3277 IDL> x0 = [2,6,2] 278 IDL> y0 = [0,2,6] 279 IDL> x1 = [3,8,4] 280 IDL> y1 = [4,4,6] 281 IDL> x2 = [1,6,4] 282 IDL> y2 = [5,6,8] 283 IDL> x3 = [0,4,2] 284 IDL> y3 = [1,4,8] 285 IDL> n = 4 286 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 287 IDL> for i=0,2 do oplot, [x0[i],x1[i],x2[i],x3[i],x0[i]],[y0[i],y1[i],y2[i],y3[i],y0[i]] 288 IDL> res=cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n) 289 IDL> for i=0,2 do oplot, [res[0,*,i]], [res[1,*,i]], color = 20+10*i, psym = 1, thick = 3 290 290 291 291 </pre><h3>Version history</h3> 292 292 293 293 <h4>Version</h4> 294 $Id: cutpar.pro 3 27 2007-12-13 16:22:35Z pinsard $294 $Id: cutpar.pro 372 2008-08-08 12:31:53Z pinsard $ 295 295 296 296 <h4>History</h4> 297 298 297 S. Masson (smasson@lodyc.jussieu.fr) 298 - July 5th, 2002 299 299 300 300 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/cutsegment.html
r338 r402 210 210 211 211 <h3>Examples</h3><pre> 212 IDL> x0=[2,5] 213 IDL> y0=[5,1] 214 IDL> x1=[9,3] 215 IDL> y1=[1,8] 216 IDL> res=cutsegment(x0, y0, x1, y1, 10) 217 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 218 IDL> oplot, [x0[0], x1[0]], [y0[0], y1[0]] 219 IDL> oplot, [res[0,*,0]], [res[1,*,0]], color = 20, psym = 1, thick = 3 220 IDL> oplot, [x0[1], x1[1]], [y0[1], y1[1]] 221 IDL> oplot, [res[0,*,1]], [res[1,*,1]], color = 40, psym = 1, thick = 3 212 213 IDL> x0=[2,5] 214 IDL> y0=[5,1] 215 IDL> x1=[9,3] 216 IDL> y1=[1,8] 217 IDL> res=cutsegment(x0, y0, x1, y1, 10) 218 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 219 IDL> oplot, [x0[0], x1[0]], [y0[0], y1[0]] 220 IDL> oplot, [res[0,*,0]], [res[1,*,0]], color = 20, psym = 1, thick = 3 221 IDL> oplot, [x0[1], x1[1]], [y0[1], y1[1]] 222 IDL> oplot, [res[0,*,1]], [res[1,*,1]], color = 40, psym = 1, thick = 3 222 223 223 224 </pre><h3>Version history</h3> 224 225 225 226 <h4>Version</h4> 226 $Id: cutsegment.pro 3 27 2007-12-13 16:22:35Z pinsard $227 $Id: cutsegment.pro 371 2008-08-07 09:32:02Z pinsard $ 227 228 228 229 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/extrapolate.html
r338 r402 76 76 </font></h2> 77 77 78 <p><font face="Courier"><i>result = </i>extrapolate(<i><a href="#_extrapolate_keyword_zinput">zinput</a>, <a href="#_extrapolate_keyword_maskinput">maskinput</a>[, <a href="#_extrapolate_keyword_nb_iteration">nb_iteration</a>]</i>, <a href="#_extrapolate_keyword_ X_PERIODIC">X_PERIODIC</a>=<i>scalar, 0 or 1</i>, <a href="#_extrapolate_keyword_MINVAL">MINVAL</a>=<i>scalar</i>, <a href="#_extrapolate_keyword_MAXVAL">MAXVAL</a>=<i>scalar</i>, <a href="#_extrapolate_keyword_GE0">GE0</a>=<i>scalar 0 or 1</i>, <a href="#_extrapolate_keyword__EXTRA">_EXTRA</a>=<i>_EXTRA</i>)</font></p>78 <p><font face="Courier"><i>result = </i>extrapolate(<i><a href="#_extrapolate_keyword_zinput">zinput</a>, <a href="#_extrapolate_keyword_maskinput">maskinput</a>[, <a href="#_extrapolate_keyword_nb_iteration">nb_iteration</a>]</i>, <a href="#_extrapolate_keyword_FILLXDIR">FILLXDIR</a>=<i>scalar, 0 or 1</i>, <a href="#_extrapolate_keyword_FILLYDIR">FILLYDIR</a>=<i>scalar, 0 or 1</i>, <a href="#_extrapolate_keyword_X_PERIODIC">X_PERIODIC</a>=<i>scalar, 0 or 1</i>, <a href="#_extrapolate_keyword_MINVAL">MINVAL</a>=<i>scalar</i>, <a href="#_extrapolate_keyword_MAXVAL">MAXVAL</a>=<i>scalar</i>, <a href="#_extrapolate_keyword_GE0">GE0</a>=<i>scalar 0 or 1</i>, <a href="#_extrapolate_keyword__EXTRA">_EXTRA</a>=<i>_EXTRA</i>)</font></p> 79 79 80 80 … … 147 147 148 148 149 <a name="#_extrapolate_keyword_FILLXDIR"></a> 150 <h4>FILLXDIR 151 152 153 154 155 156 <font size="-1" color="#006633">type:</font> <font size="-1" color="#006633"><i>scalar, 0 or 1</i></font> 157 <font size="-1" color="#006633">default:</font> <font size="-1" color="#006633"><i>0</i></font> 158 159 </h4> 160 161 162 put 1 to specify that filling of the data must be done only along x direction 163 164 165 <a name="#_extrapolate_keyword_FILLYDIR"></a> 166 <h4>FILLYDIR 167 168 169 170 171 172 <font size="-1" color="#006633">type:</font> <font size="-1" color="#006633"><i>scalar, 0 or 1</i></font> 173 <font size="-1" color="#006633">default:</font> <font size="-1" color="#006633"><i>0</i></font> 174 175 </h4> 176 177 178 put 1 to specify that filling of the data must be done only along y direction 179 180 149 181 <a name="#_extrapolate_keyword_X_PERIODIC"></a> 150 182 <h4>X_PERIODIC … … 229 261 230 262 <h3>Examples</h3><pre> 231 IDL> a=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 232 IDL> tvplus, a 233 IDL> tvplus, a*(1-tmask[*,*,0]) 263 264 IDL> a=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 265 IDL> tvplus, a 266 IDL> tvplus, a*(1-tmask[*,*,0]) 267 234 268 get the coastline: 235 IDL> a=extrapolate(tmask[*,*,0],tmask[*,*,0],1,/x_periodic) 236 IDL> tvplus, a-tmask[*,*,0] 269 270 IDL> a=extrapolate(tmask[*,*,0],tmask[*,*,0],1,/x_periodic) 271 IDL> tvplus, a-tmask[*,*,0] 237 272 238 273 </pre><h3>Version history</h3> 239 274 240 275 <h4>Version</h4> 241 $Id: extrapolate.pro 3 27 2007-12-13 16:22:35Z pinsard$276 $Id: extrapolate.pro 384 2008-11-06 02:32:33Z smasson $ 242 277 243 278 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/extrapsmooth.html
r338 r402 211 211 212 212 <h3>Examples</h3><pre> 213 IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 214 IDL> tvplus, a 213 214 IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 215 IDL> tvplus, a 216 215 217 compare to extrapolate result: 216 IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 217 IDL> tvplus, b, window = 1 218 219 IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) 220 IDL> tvplus, b, window = 1 218 221 219 222 </pre><h3>Version history</h3> 220 223 221 224 <h4>Version</h4> 222 $Id: extrapsmooth.pro 3 27 2007-12-13 16:22:35Z pinsard $225 $Id: extrapsmooth.pro 371 2008-08-07 09:32:02Z pinsard $ 223 226 <h4>History</h4> 224 227 January 2007: Sebastien Masson (smasson@lodyc.jussieu.fr) -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/file_interp.html
r363 r402 703 703 <h3>Examples</h3><pre> 704 704 705 IDL> file_interp, filein, fileout, gridout, inxaxisname = 'lo', inyaxisname = 'la', keep = ['lo', 'la', 'cond_sed'] 706 707 IDL> file_interp, in, out, gdout, inuseasmask = 'sst', inmissing_value = -1.00000e+30, missing_value = -1000.00 705 IDL> file_interp, filein, fileout, gridout, inxaxisname = 'lo', inyaxisname = 'la', keep = ['lo', 'la', 'cond_sed'] 706 IDL> file_interp, in, out, gdout, inuseasmask = 'sst', inmissing_value = -1.00000e+30, missing_value = -1000.00 708 707 709 708 </pre><h3>Version history</h3> 710 709 711 710 <h4>Version</h4> 712 $Id: file_interp.pro 3 62 2008-07-29 09:00:03Z smasson $711 $Id: file_interp.pro 383 2008-10-30 13:55:33Z smasson $ 713 712 714 713 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/fromirr.html
r338 r402 63 63 64 64 interpolate data from an irregular 2D grid to any 2D grid. 65 Only 1 method available = bilinear 65 66 Only 1 method available = bilinear 66 67 67 68 … … 268 269 dataout = total(weig*datain[addr], 1) 269 270 dataout = reform(dataout, jpio, jpjo, /over) 271 270 272 In that case, method, lonin, latin, are not used (but are necessary). 271 273 lonout, latout are used only to know the output domain size … … 293 295 To interpolate 1 field: 294 296 295 IDL> tncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout)297 IDL> tncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout) 296 298 297 299 or if you have several fields to interpolate from the same source and target grids … … 300 302 (that must be undefined or equal to 0 before calling fromirr) 301 303 302 IDL> t1ncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout $304 IDL> t1ncep = fromirr('bilinear', topa, glamt, gphit, tmask[*,*,0], lonout, latout, mskout $ 303 305 , WEIG = a, ADDR = b) 304 IDL> help, a, b306 IDL> help, a, b 305 307 306 308 2) use a and b that are now defined to bypass the computation of the weights 307 309 and addresses and speed-up the computation! 308 310 309 IDL> t2ncep = fromirr('bilinear', topa, WEIG = a, ADDR = b)311 IDL> t2ncep = fromirr('bilinear', topa, WEIG = a, ADDR = b) 310 312 311 313 </pre><h3>Version history</h3> 312 314 313 315 <h4>Version</h4> 314 $Id: fromirr.pro 3 27 2007-12-13 16:22:35Z pinsard $316 $Id: fromirr.pro 372 2008-08-08 12:31:53Z pinsard $ 315 317 316 318 <h4>History</h4> 317 319 June 2006: Sebastien Masson (smasson@lodyc.jussieu.fr) 318 320 319 321 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/fromreg.html
r338 r402 63 63 64 64 interpolate data from a "regular/rectangular grid" to any grid. 65 2 methods available: bilinear and imoms3 66 A "regular/rectangular grid" is defined as a grid for which 67 each longitude lines have the same latitude and each latitude columns 65 66 2 methods available: bilinear and imoms3 67 68 A "regular/rectangular grid" is defined as a grid for which 69 70 Each longitude lines have the same latitude and each latitude columns 68 71 have the same longitude. 69 72 … … 236 239 dataout = total(weig*datain[addr], 1) 237 240 dataout = reform(dataout, jpio, jpjo, /over) 241 238 242 In that case, method, lonin, latin, are not used (but are necessary). 239 243 lonout, latout are used only to know the output domain size … … 295 299 To interpolate 1 field: 296 300 297 IDL> topa = fromreg('bilinear', tncep, xncep, yncep, glamt, gphit)301 IDL> topa = fromreg('bilinear', tncep, xncep, yncep, glamt, gphit) 298 302 299 303 or if you have several fields to interpolate from the same source and target grids … … 302 306 (that must be undefined or equal to 0 before calling fromreg 303 307 304 IDL> t1opa = fromreg('bilinear', t1ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b)305 IDL> help, a, b308 IDL> t1opa = fromreg('bilinear', t1ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b) 309 IDL> help, a, b 306 310 307 311 2) use a and b that are now defined to bypass the computation of the weights and addresses 308 312 and speed-up the computation! 309 313 310 IDL> t2opa = fromreg('bilinear', t2ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b)314 IDL> t2opa = fromreg('bilinear', t2ncep, xncep, yncep, glamt, gphit, WEIG = a, ADDR = b) 311 315 312 316 </pre><h3>Version history</h3> 313 317 314 318 <h4>Version</h4> 315 $Id: fromreg.pro 3 27 2007-12-13 16:22:35Z pinsard $319 $Id: fromreg.pro 372 2008-08-08 12:31:53Z pinsard $ 316 320 317 321 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/get_gridparams.html
r338 r402 262 262 263 263 Case 1: 264 IDL> get_gridparams, file name/id, lonname, latname, lon, lat, jpi, jpj, n_dimensions264 IDL> get_gridparams, file name/id, lonname, latname, lon, lat, jpi, jpj, n_dimensions 265 265 266 266 Case 2: 267 IDL> get_gridparams, lon, lat, jpi, jpj, n_dimensions 268 269 270 271 1) IDL> ncdf_get_gridparams, 'coordinates_ORCA_R05.nc', 'glamt', 'gphit' $ 267 IDL> get_gridparams, lon, lat, jpi, jpj, n_dimensions 268 269 270 271 1) 272 273 IDL> ncdf_get_gridparams, 'coordinates_ORCA_R05.nc', 'glamt', 'gphit' $ 272 274 , olon, olat, jpio, jpjo, 2 273 275 274 2) IDL> ncdf_get_gridparams, olon, olat, jpio, jpjo, 2 276 2) 277 278 IDL> ncdf_get_gridparams, olon, olat, jpio, jpjo, 2 275 279 276 280 </pre><h3>Version history</h3> 277 281 278 282 <h4>Version</h4> 279 $Id: get_gridparams.pro 3 27 2007-12-13 16:22:35Z pinsard$283 $Id: get_gridparams.pro 379 2008-09-03 15:32:14Z smasson $ 280 284 281 285 <h4>History</h4> 282 286 November 2005: Sebastien Masson (smasson@lodyc.jussieu.fr) 283 287 284 288 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/inquad.html
r338 r402 339 339 <h3>Examples</h3><pre> 340 340 341 IDL> x = 1.*[1, 2, 6, 7, 3]342 IDL> y = 1.*[1, 3, 3, 4, 7]343 IDL> x1 = 1.*[0,4,2]344 IDL> y1 = 1.*[1,4,8]345 IDL> x2 = 1.*[1,6,4]346 IDL> y2 = 1.*[5,6,8]347 IDL> x3 = 1.*[3,8,4]348 IDL> y3 = 1.*[4,4,6]349 IDL> x4 = 1.*[2,6,2]350 IDL> y4 = 1.*[0,2,6]351 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata352 IDL> for i=0,2 do oplot, [x4[i],x1[i],x2[i],x3[i],x4[i]],[y4[i],y1[i],y2[i],y3[i],y4[i]]353 IDL> oplot, x, y, color = 20, psym = 1, thick = 2354 IDL> print, inquad(x, y, x1, y1, x2, y2, x3, y3, x4, y4)341 IDL> x = 1.*[1, 2, 6, 7, 3] 342 IDL> y = 1.*[1, 3, 3, 4, 7] 343 IDL> x1 = 1.*[0,4,2] 344 IDL> y1 = 1.*[1,4,8] 345 IDL> x2 = 1.*[1,6,4] 346 IDL> y2 = 1.*[5,6,8] 347 IDL> x3 = 1.*[3,8,4] 348 IDL> y3 = 1.*[4,4,6] 349 IDL> x4 = 1.*[2,6,2] 350 IDL> y4 = 1.*[0,2,6] 351 IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata 352 IDL> for i=0,2 do oplot, [x4[i],x1[i],x2[i],x3[i],x4[i]],[y4[i],y1[i],y2[i],y3[i],y4[i]] 353 IDL> oplot, x, y, color = 20, psym = 1, thick = 2 354 IDL> print, inquad(x, y, x1, y1, x2, y2, x3, y3, x4, y4) 355 355 356 356 On a sphere see … … 360 360 361 361 <h4>Version</h4> 362 $Id: inquad.pro 3 27 2007-12-13 16:22:35Z pinsard $362 $Id: inquad.pro 371 2008-08-07 09:32:02Z pinsard $ 363 363 364 364 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/inrecgrid.html
r338 r402 62 62 63 63 64 given - a list of points, (x,y) position 65 - the x and y limits of a rectangular grid 64 given 65 - a list of points, (x,y) position 66 - the x and y limits of a rectangular grid 67 66 68 find in which cell is located each given point. 67 69 … … 201 203 <h3>Examples</h3><pre> 202 204 203 IDL> a=indgen(5)204 IDL> b=indgen(7)205 IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,b)206 IDL> print, r205 IDL> a=indgen(5) 206 IDL> b=indgen(7) 207 IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,b) 208 IDL> print, r 207 209 20 13 7 208 IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,a+1,b,b+1,/output2d)209 IDL> print, r210 IDL> r=inrecgrid([0.25,3.25,2],[4.25,2.8,1.4],a,a+1,b,b+1,/output2d) 211 IDL> print, r 210 212 0.00000 4.00000 211 213 3.00000 2.00000 … … 215 217 216 218 <h4>Version</h4> 217 $Id: inrecgrid.pro 3 27 2007-12-13 16:22:35Z pinsard $219 $Id: inrecgrid.pro 372 2008-08-08 12:31:53Z pinsard $ 218 220 219 221 <h4>History</h4> 220 222 S. Masson (smasson@lodyc.jussieu.fr) 221 July 3rd, 2002222 October 3rd, 2003: use value_locate223 - July 3rd, 2002 224 - October 3rd, 2003: use value_locate 223 225 224 226 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/lbcorca.html
r338 r402 178 178 179 179 <h3>Examples</h3><pre> 180 IDL> @tst_initorca2 181 IDL> a = dist(jpiglo, jpjglo) 182 IDL> print, lbcorca(a, 'T') 183 IDL> dummy = lbcorca(a, 'T', /correction) 184 IDL> print, lbcorca(a, 'T') 180 181 IDL> @tst_initorca2 182 IDL> a = dist(jpiglo, jpjglo) 183 IDL> print, lbcorca(a, 'T') 184 IDL> dummy = lbcorca(a, 'T', /correction) 185 IDL> print, lbcorca(a, 'T') 185 186 186 187 </pre><h3>Version history</h3> 187 188 188 189 <h4>Version</h4> 189 $Id: lbcorca.pro 3 31 2008-02-11 16:46:48Z smasson$190 $Id: lbcorca.pro 371 2008-08-07 09:32:02Z pinsard $ 190 191 191 192 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/ll_narcs_distances.html
r338 r402 193 193 194 194 <h3>Examples</h3><pre> 195 IDL> Lon_lat0 = [1.0, 2.0]; Initial point specified in radians 196 IDL> Arc_Dist = 2.0; Arc distance in radians 197 IDL> Az = 1.0; Azimuth in radians 198 IDL> Result = LL_ARC_DISTANCE(Lon_lat0, Arc_Dist, Az) 199 IDL> PRINT, Result 195 196 IDL> Lon_lat0 = [1.0, 2.0]; Initial point specified in radians 197 IDL> Arc_Dist = 2.0; Arc distance in radians 198 IDL> Az = 1.0; Azimuth in radians 199 IDL> Result = LL_ARC_DISTANCE(Lon_lat0, Arc_Dist, Az) 200 IDL> PRINT, Result 200 201 2.91415 -0.622234 201 202 202 IDL> lon0 = [-10, 20, 100]203 IDL> lat0 = [0, -10, 45]204 IDL> lon1 = [10, 60, 280]205 IDL> lat1 = [0, 10, 45]206 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, azimuth = azi, /two_by_two)207 IDL> earthradius = 6378206.4d0208 IDL> res = ll_narcs_distances(lon0, lat0, dist/earthradius, azi, /degrees)209 IDL> print, reform(res[0, *])203 IDL> lon0 = [-10, 20, 100] 204 IDL> lat0 = [0, -10, 45] 205 IDL> lon1 = [10, 60, 280] 206 IDL> lat1 = [0, 10, 45] 207 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, azimuth = azi, /two_by_two) 208 IDL> earthradius = 6378206.4d0 209 IDL> res = ll_narcs_distances(lon0, lat0, dist/earthradius, azi, /degrees) 210 IDL> print, reform(res[0, *]) 210 211 10.000000 60.000000 280.00000 211 IDL> print, reform(res[1, *])212 IDL> print, reform(res[1, *]) 212 213 1.1999280e-15 10.000000 45.000000 213 214 … … 215 216 216 217 <h4>Version</h4> 217 $Id: ll_narcs_distances.pro 3 27 2007-12-13 16:22:35Z pinsard $218 $Id: ll_narcs_distances.pro 371 2008-08-07 09:32:02Z pinsard $ 218 219 219 220 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/map_npoints.html
r338 r402 253 253 254 254 <h3>Examples</h3><pre> 255 IDL> print, $256 IDL> map_npoints([-105.15,1],[40.02,1],[-0.07,100,50],[51.30,20,0])255 IDL> print, $ 256 IDL> map_npoints([-105.15,1],[40.02,1],[-0.07,100,50],[51.30,20,0]) 257 257 7551369.3 5600334.8 258 258 12864354. 10921254. 259 259 14919237. 5455558.8 260 260 261 IDL> lon0 = [-10, 20, 100]262 IDL> lat0 = [0, -10, 45]263 IDL> lon1 = [10, 60, 280]264 IDL> lat1 = [0, 10, 45]265 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi)266 IDL> help, dist, azi261 IDL> lon0 = [-10, 20, 100] 262 IDL> lat0 = [0, -10, 45] 263 IDL> lon1 = [10, 60, 280] 264 IDL> lat1 = [0, 10, 45] 265 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi) 266 IDL> help, dist, azi 267 267 DIST DOUBLE = Array[3, 3] 268 268 AZI DOUBLE = Array[3, 3] 269 IDL> print, dist[4*lindgen(3)], azi[4*lindgen(3)]269 IDL> print, dist[4*lindgen(3)], azi[4*lindgen(3)] 270 270 2226414.0 4957944.5 10018863. 271 271 90.000000 64.494450 4.9615627e-15 272 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi, /TWO_BY_TWO)273 IDL> help, dist, azi272 IDL> dist = map_npoints(lon0, lat0, lon1, lat1, AZIMUTH = azi, /TWO_BY_TWO) 273 IDL> help, dist, azi 274 274 DIST DOUBLE = Array[3] 275 275 AZI DOUBLE = Array[3] 276 IDL> print, dist, azi276 IDL> print, dist, azi 277 277 2226414.0 4957944.5 10018863. 278 278 90.000000 64.494450 4.9615627e-15 279 IDL> print, map_2points(lon0[0], lat0[0], lon1[0], lat1[0])279 IDL> print, map_2points(lon0[0], lat0[0], lon1[0], lat1[0]) 280 280 20.000000 90.000000 281 IDL> print, map_npoints(lon0[0], lat0[0], lon1[0], lat1[0], AZIMUTH=azi)/6378206.4d0 / !dtor, azi281 IDL> print, map_npoints(lon0[0], lat0[0], lon1[0], lat1[0], AZIMUTH=azi)/6378206.4d0 / !dtor, azi 282 282 20.000000 283 283 90.000000 284 284 285 IDL> lon0 = [-10, 20, 100]286 IDL> lat0 = [0, -10, 45]287 IDL> lon1 = [10, 60, 280]288 IDL> lat1 = [0, 10, 45]289 IDL> mid = map_npoints(lon0, lat0, lon1, lat1, /MIDDLE, /TWO_BY_TWO)290 IDL> print, reform(mid[0,*]), reform(mid[1,*])285 IDL> lon0 = [-10, 20, 100] 286 IDL> lat0 = [0, -10, 45] 287 IDL> lon1 = [10, 60, 280] 288 IDL> lat1 = [0, 10, 45] 289 IDL> mid = map_npoints(lon0, lat0, lon1, lat1, /MIDDLE, /TWO_BY_TWO) 290 IDL> print, reform(mid[0,*]), reform(mid[1,*]) 291 291 0.0000000 40.000000 190.00000 292 292 0.0000000 -1.5902773e-15 90.000000 293 IDL> print, (map_2points(lon0[0], lat0[0], lon1[0], lat1[0], npath = 3))[*, 1]293 IDL> print, (map_2points(lon0[0], lat0[0], lon1[0], lat1[0], npath = 3))[*, 1] 294 294 0.0000000 0.0000000 295 IDL> print, (map_2points(lon0[1], lat0[1], lon1[1], lat1[1], npath = 3))[*, 1]295 IDL> print, (map_2points(lon0[1], lat0[1], lon1[1], lat1[1], npath = 3))[*, 1] 296 296 40.000000 -1.5902773e-15 297 IDL> print, (map_2points(lon0[2], lat0[2], lon1[2], lat1[2], npath = 3))[*, 1]297 IDL> print, (map_2points(lon0[2], lat0[2], lon1[2], lat1[2], npath = 3))[*, 1] 298 298 190.00000 90.000000 299 299 … … 301 301 302 302 <h4>Version</h4> 303 $Id: map_npoints.pro 3 27 2007-12-13 16:22:35Z pinsard $303 $Id: map_npoints.pro 371 2008-08-07 09:32:02Z pinsard $ 304 304 305 305 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/neighbor.html
r338 r402 211 211 212 212 <h3>Examples</h3><pre> 213 IDL> print, neighbor(-105.15,40.02,[-0.07,100,50],[51.30,20,0], $ 214 IDL> distance=dis) 213 214 IDL> print, neighbor(-105.15,40.02,[-0.07,100,50],[51.30,20,0], $ 215 distance=dis) 215 216 0 216 IDL> print, dis217 IDL> print, dis 217 218 105.684 206.125 160.228 218 219 … … 220 221 221 222 <h4>Version</h4> 222 $Id: neighbor.pro 3 27 2007-12-13 16:22:35Z pinsard $223 $Id: neighbor.pro 371 2008-08-07 09:32:02Z pinsard $ 223 224 224 225 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/quadrilateral2square.html
r363 r402 64 64 warm (or map) an arbitrary quadrilateral onto a unit square 65 65 according to the 4-point correspondences: 66 (x0,y0) -> (0,0) 67 (x1,y1) -> (1,0) 68 (x2,y2) -> (1,1) 69 (x3,y3) -> (0,1) 66 - (x0,y0) -> (0,0) 67 - (x1,y1) -> (1,0) 68 - (x2,y2) -> (1,1) 69 - (x3,y3) -> (0,1) 70 70 71 This is the inverse function of <a href="..//Interpolation/square2quadrilateral.html">square2quadrilateral</a>. 71 72 … … 302 303 <h3>Examples</h3><pre> 303 304 304 IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1305 IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20306 IDL> xin = (findgen(11)*.1)#replicate(1, 11)307 IDL> yin = replicate(1, 11)#(findgen(11)*.1)308 IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin)309 IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20310 311 IDL> inorg=quadrilateral2square(2,1,3,0,5,1,2,3,out[0,*],out[1,*])312 IDL> tracegrille, reform(inorg[0,*],11,11), reform(inorg[1,*],11,11),color=indgen(12)*20305 IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1 306 IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20 307 IDL> xin = (findgen(11)*.1)#replicate(1, 11) 308 IDL> yin = replicate(1, 11)#(findgen(11)*.1) 309 IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin) 310 IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20 311 312 IDL> inorg=quadrilateral2square(2,1,3,0,5,1,2,3,out[0,*],out[1,*]) 313 IDL> tracegrille, reform(inorg[0,*],11,11), reform(inorg[1,*],11,11),color=indgen(12)*20 313 314 314 315 </pre><h3>Version history</h3> 315 316 316 317 <h4>Version</h4> 317 $Id: quadrilateral2square.pro 3 58 2008-04-28 09:41:14Z pinsard $318 $Id: quadrilateral2square.pro 372 2008-08-08 12:31:53Z pinsard $ 318 319 319 320 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/spl_fstdrv.html
r338 r402 86 86 87 87 <h3>Return value</h3> 88 88 y2: f'(x2) = y2. 89 89 90 90 … … 167 167 168 168 <h3>Examples</h3><pre> 169 IDL> y2 = spl_fstdrv(x, y, yscd, x2) 169 170 IDL> y2 = spl_fstdrv(x, y, yscd, x2) 170 171 171 172 </pre><h3>Version history</h3> 172 173 173 174 <h4>Version</h4> 174 $Id: spl_fstdrv.pro 3 25 2007-12-06 10:04:53Z pinsard $175 $Id: spl_fstdrv.pro 372 2008-08-08 12:31:53Z pinsard $ 175 176 176 177 <h4>History</h4> 177 178 Sebastien Masson (smasson@lodyc.jussieu.fr): May 2005 178 179 179 180 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/spl_incr.html
r338 r402 213 213 214 214 <h3>Examples</h3><pre> 215 IDL> n = 100L 216 IDL> x = (dindgen(n))^2 217 IDL> y = abs(randomn(0, n)) 218 IDL> y[n/2:n/2+1] = 0. 219 IDL> y[n-n/3] = 0. 220 IDL> y[n-n/6:n-n/6+5] = 0. 221 IDL> y = total(y, /cumulative, /double) 222 IDL> x2 = dindgen((n-1)^2) 223 IDL> n2 = n_elements(x2) 224 IDL> print, min(y[1:n-1]-y[0:n-2]) LT 0 225 IDL> y2 = spl_incr( x, y, x2) 226 IDL> splot, x, y, xstyle = 1, ystyle = 1, ysurx=.25, petit = [1, 2, 1], /land 227 IDL> oplot, x2, y2, color = 100 228 IDL> c = y2[1:n2-1] - y2[0:n2-2] 229 IDL> print, min(c) LT 0 230 IDL> print, min(c, max = ma), ma 231 IDL> splot,c,xstyle=1,ystyle=1, yrange=[-.01,.05], ysurx=.25, petit = [1, 2, 2], /noerase 232 IDL> oplot,[0, n_elements(c)], [0, 0], linestyle = 1 215 216 IDL> n = 100L 217 IDL> x = (dindgen(n))^2 218 IDL> y = abs(randomn(0, n)) 219 IDL> y[n/2:n/2+1] = 0. 220 IDL> y[n-n/3] = 0. 221 IDL> y[n-n/6:n-n/6+5] = 0. 222 IDL> y = total(y, /cumulative, /double) 223 IDL> x2 = dindgen((n-1)^2) 224 IDL> n2 = n_elements(x2) 225 IDL> print, min(y[1:n-1]-y[0:n-2]) LT 0 226 IDL> y2 = spl_incr( x, y, x2) 227 IDL> splot, x, y, xstyle = 1, ystyle = 1, ysurx=.25, petit = [1, 2, 1], /land 228 IDL> oplot, x2, y2, color = 100 229 IDL> c = y2[1:n2-1] - y2[0:n2-2] 230 IDL> print, min(c) LT 0 231 IDL> print, min(c, max = ma), ma 232 IDL> splot,c,xstyle=1,ystyle=1, yrange=[-.01,.05], ysurx=.25, petit = [1, 2, 2], /noerase 233 IDL> oplot,[0, n_elements(c)], [0, 0], linestyle = 1 233 234 234 235 </pre><h3>Version history</h3> 235 236 236 237 <h4>Version</h4> 237 $Id: spl_incr.pro 3 27 2007-12-13 16:22:35Z pinsard $238 $Id: spl_incr.pro 371 2008-08-07 09:32:02Z pinsard $ 238 239 239 240 <h4>History</h4> -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/spl_keep_mean.html
r338 r402 209 209 <h3>Examples</h3><pre> 210 210 211 212 213 IDL> yr1 = 1990214 IDL> yr2 = 1992215 IDL> nyr = yr2-yr1+1216 IDL> n1 = 12*nyr+1217 IDL> x = julday(1+findgen(n1), replicate(1, n1) $218 IDL> , replicate(yr1, n1), fltarr(n1))219 IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 1220 IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $221 IDL> , replicate(yr1, n2), fltarr(n2))222 IDL> y = abs(randomn(0, n1-1))223 IDL> y2 = spl_keep_mean(x, y, x2, /ge0)224 225 IDL> print, min(x, max = ma), ma226 IDL> print, min(x2, max = ma), ma227 IDL> print, vairdate([min(x, max = ma), ma])228 IDL> print, total(y*(x[1:n1-1]-x[0:n1-2]))229 IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2]))211 12 monthly values of precipitations into daily values: 212 213 IDL> yr1 = 1990 214 IDL> yr2 = 1992 215 IDL> nyr = yr2-yr1+1 216 IDL> n1 = 12*nyr+1 217 IDL> x = julday(1+findgen(n1), replicate(1, n1) $ 218 IDL> , replicate(yr1, n1), fltarr(n1)) 219 IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 1 220 IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $ 221 IDL> , replicate(yr1, n2), fltarr(n2)) 222 IDL> y = abs(randomn(0, n1-1)) 223 IDL> y2 = spl_keep_mean(x, y, x2, /ge0) 224 225 IDL> print, min(x, max = ma), ma 226 IDL> print, min(x2, max = ma), ma 227 IDL> print, vairdate([min(x, max = ma), ma]) 228 IDL> print, total(y*(x[1:n1-1]-x[0:n1-2])) 229 IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2])) 230 230 231 231 </pre><h3>Version history</h3> 232 232 233 233 <h4>Version</h4> 234 $Id: spl_keep_mean.pro 3 27 2007-12-13 16:22:35Z pinsard $234 $Id: spl_keep_mean.pro 372 2008-08-08 12:31:53Z pinsard $ 235 235 236 236 <h4>History</h4> 237 237 Sebastien Masson (smasson@lodyc.jussieu.fr): May 2005 238 238 239 239 -
trunk/SRC/Documentation/idldoc_assistant_output/Interpolation/square2quadrilateral.html
r363 r402 64 64 warm (or map) a unit square onto an arbitrary quadrilateral 65 65 according to the 4-point correspondences: 66 (0,0) -> (x0,y0) 67 (1,0) -> (x1,y1) 68 (1,1) -> (x2,y2) 69 (0,1) -> (x3,y3) 66 - (0,0) -> (x0,y0) 67 - (1,0) -> (x1,y1) 68 - (1,1) -> (x2,y2) 69 - (0,1) -> (x3,y3) 70 70 71 The mapping is done using perspective transformation which preserve 71 72 lines in all orientations and permit quadrilateral to quadrilateral … … 284 285 <h3>Examples</h3><pre> 285 286 286 IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1287 IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20288 IDL> xin = (findgen(11)*.1)#replicate(1, 11)289 IDL> yin = replicate(1, 11)#(findgen(11)*.1)290 IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin)291 IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20287 IDL> splot,[0,5],[0,3],/nodata,xstyle=1,ystyle=1 288 IDL> tracegrille, findgen(11)*.1, findgen(11)*.1,color=indgen(12)*20 289 IDL> xin = (findgen(11)*.1)#replicate(1, 11) 290 IDL> yin = replicate(1, 11)#(findgen(11)*.1) 291 IDL> out = square2quadrilateral(2,1,3,0,5,1,2,3, xin, yin) 292 IDL> tracegrille, reform(out[0,*],11,11), reform(out[1,*],11,11),color=indgen(12)*20 292 293 293 294 </pre><h3>Version history</h3> 294 295 295 296 <h4>Version</h4> 296 $Id: square2quadrilateral.pro 3 58 2008-04-28 09:41:14Z pinsard $297 $Id: square2quadrilateral.pro 372 2008-08-08 12:31:53Z pinsard $ 297 298 298 299 <h4>History</h4> 299 Sebastien Masson (smasson@lodyc.jussieu.fr) 300 August 2003 301 Based on "Digital Image Warping" by G. Wolberg 302 IEEE Computer Society Press, Los Alamitos, California 303 Chapter 3, see p 52-56 304 300 Sebastien Masson (smasson@lodyc.jussieu.fr) 301 - August 2003 302 Based on "Digital Image Warping" by G. Wolberg 303 IEEE Computer Society Press, Los Alamitos, California 304 Chapter 3, see p 52-56 305 305 306 306
Note: See TracChangeset
for help on using the changeset viewer.