/usr/share/gnudatalanguage/lib/congrid.pro is in libgnudatalanguage0 0.9.7-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 | ; part of GNU Data Language - GDL
;-------------------------------------------------------------
;+
; NAME:
; CONGRID
; PURPOSE:
; Shrink or expands an array (any dimensions) to any size,
; however only the first 3 dims (max) will be resized.
; CATEGORY:
; Array Manipulation.
; CALLING SEQUENCE:
; t2 = CONGRID(t,mx[,my[,mZ]])
; INPUTS:
; t = input array. in
; mx, my, mz = new (output) array size. in
; t2 = resulting array. out
; KEYWORD PARAMETERS:
; /INTERPOLATE means do bilinear interpolation. Kept for
; compatibility but unused: interpolation IS linear.
; CUBIC=cubic: use a cubic interpolation. Currently limited by
; the availability of cubic interpolation in GDL's INTERPOLATE.
; /CENTERED means assume pixels centered. This means
; the pixel at (0,[0,[0]]) is clipped to 1/4 size.
; Default is that pixel start (not center) is at index.'
; MISSING=missing gives a value for undefined pixels.
; /MINUS_ONE: option will be ignored. MISSING can be used instead.
; /HELP gives this help.'
;
; NOTE:
; CONGRID performs a resampling. Does not conserve Fluxes.
; OUTPUTS:
; t2 = resulting array. out
;
; MODIFICATION HISTORY:
; Written by: 2011-12-17 G. Duvert (Initial import)
;
;-
; LICENCE:
; Copyright (C) 2011
; This program is free software; you can redistribute it and/or modify
; it under the terms of the GNU General Public License as published by
; the Free Software Foundation; either version 2 of the License, or
; (at your option) any later version.
;
;-
function CONGRID, t, mx, my, mz, $
CENTER=center, $
INTERP=interp, $
CUBIC = cubic, $
MINUS_ONE=minus_one, $
MISSING=missing, $
HELP=help, test=test
npar=N_PARAMS()
if (npar lt 2) or (npar gt 4) or KEYWORD_SET(help) then begin
print,' Shrink or expands an array (any dimensions) to any size,'
print,' however only the first 3 dims (max) will be resized.'
print,' CALLING SEQUENCE:'
print,' t2 = congrid2(t,mx[,my][,mz])'
print,' t = input array. in'
print,' mx, my, mz = new (output) array size. in'
print,' t2 = resulting array. out'
print,' Keywords:'
print,' /INTERPOLATE means do bilinear interpolation. Kept for'
print,' compatibility but unused: interpolation IS linear.'
print,' /CENTERED means assume pixels centered. This means'
print,' the pixel at (0,[0,[0]]) is clipped to 1/4 size.'
print,' Default is that pixel start (not center) is at index.'
print,' CUBIC=cubic: use a cubic interpolation. Currently limited by'
print,' the availability of cubic interpolation in INTERPOLATE.'
print,' /CENTERED means assume pixels centered. This means'
print,' the pixel at (0,[0,[0]]) is clipped to 1/4 size.'
print,' Default is that pixel start (not center) is at index.'
print,' MISSING=missing gives a value for undefined pixels.'
print,' /MINUS_ONE: option will be ignored. MISSING can be used instead.'
print, ' /HELP gives this help.'
print,' NOTE: CONGRID performs a resampling. Does not conserve Fluxes.'
return,''
endif
;
ON_ERROR, 2 ;Return to caller if error
;
ndim = SIZE(t, /N_DIMENSIONS)
dims = SIZE(t, /DIMENSIONS)
;
if(ndim gt 3) then $
MESSAGE, 'Warning, array dimensions >3, using non-standard CONGRID extension.', /informational
;
cflag = KEYWORD_SET(centered) ; 0=start, 1=centered.
minone = KEYWORD_SET(minus_one)
nnbor = ~KEYWORD_SET(interp) ; nearest neighbor is default for 1 and 2 dim
;
; if (npar gt ndim+1) then message, 'Too many output dimensions.'
;
nopt=0
IF (N_ELEMENTS(mx) NE 0) THEN BEGIN
fx= (minone gt 0)? 1.0D*(dims[0]-1)/(mx-1) : 1.0D*dims[0]/mx
x = (cflag gt 0) ?fx*FINDGEN(mx)+0.5 : fx*FINDGEN(mx)
nopt=1
ENDIF
;
; if 2 dim my MUST be present!
IF (N_ELEMENTS(my) EQ 0 && ndim ge 2) THEN MESSAGE, '2nd dimension not specified.'
IF (N_ELEMENTS(my) NE 0 && ndim ge 2) THEN BEGIN
fy= (minone gt 0)? 1.0D*(dims[1]-1)/(my-1) : 1.0D*dims[1]/my
y = (cflag gt 0) ?fy*FINDGEN(my)+0.5 : fy*FINDGEN(my)
nopt=2
ENDIF
;
; if 3 dim mz MUST be present!
IF (N_ELEMENTS(mz) EQ 0 && ndim ge 3) THEN MESSAGE, '3rd dimension not specified.'
IF (N_ELEMENTS(mz) NE 0 && ndim ge 3) THEN BEGIN
fz= (minone gt 0) ? 1.0D*(dims[2]-1)/(mz-1) : 1.0D*dims[2]/mz
z = (cflag gt 0) ?fz*FINDGEN(mz)+0.5 : fz*FINDGEN(mz)
nopt=3
ENDIF
;
IF (nopt lt ndim) THEN BEGIN
;; create a transposition vector
p=SHIFT(INDGEN(ndim),-nopt)
temp=TRANSPOSE(t,p)
CASE nopt OF
3: temp2 = INTERPOLATE(temp,x,y,z,/grid,cubic=cubic,missing=missing)
2: temp2 = INTERPOLATE(temp,x,y,/grid,nearest_neighbour=nnbor,cubic=cubic,missing=missing)
1: temp2 = INTERPOLATE(temp,x,/grid,nearest_neighbour=nnbor,cubic=cubic,missing=missing)
ENDCASE
p=SHIFT(INDGEN(ndim),nopt)
t2=TRANSPOSE(temp2,p)
ENDIF ELSE BEGIN
CASE nopt OF
3: t2 = INTERPOLATE(t,x,y,z,/grid,cubic=cubic,missing=missing)
2: t2 = INTERPOLATE(t,x,y,/grid,nearest_neighbour=nnbor,cubic=cubic,missing=missing)
1: t2 = INTERPOLATE(t,x,/grid,nearest_neighbour=nnbor,cubic=cubic,missing=missing)
ENDCASE
ENDELSE
;
if KEYWORD_SET(test) then STOP
;
return,t2
;
end
|