(Perl Data Language) aims to turn perl into an efficient numerical language for scientific computing (similar to IDL and MatLab).
| Image2D.pm (PDL-2.082) | : | Image2D.pm (PDL-2.083) |
|---|
| | |
| skipping to change at line 42 | | skipping to change at line 42 |
|---|
| | |
| use PDL::Image2D; | | use PDL::Image2D; |
| | |
| =cut | | =cut |
| | |
| use PDL; # ensure qsort routine available | | use PDL; # ensure qsort routine available |
| use PDL::Math; | | use PDL::Math; |
| use Carp; | | use Carp; |
| | |
| my %boundary2value = (Reflect=>1, Truncate=>2, Replicate=>3); | | my %boundary2value = (Reflect=>1, Truncate=>2, Replicate=>3); |
|
| #line 50 "Image2D.pm" | | #line 51 "Image2D.pm" |
| | |
| =head1 FUNCTIONS | | =head1 FUNCTIONS |
| | |
| =cut | | =cut |
| | |
|
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| #line 65 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| #line 70 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| #line 75 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | | |
| =head2 conv2d | | =head2 conv2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); kern(p,q); [o]b(m,n); indx [t]mapi(isize); indx [t]mapj(js
ize); int opt) | | Signature: (a(m,n); kern(p,q); [o]b(m,n); indx [t]mapi(isize); indx [t]mapj(js
ize); int opt) |
| | |
| =for ref | | =for ref |
| | |
| 2D convolution of an array with a kernel (smoothing) | | 2D convolution of an array with a kernel (smoothing) |
| | |
| | |
| skipping to change at line 96 | | skipping to change at line 85 |
|---|
| (i.e. wrap around axis) | | (i.e. wrap around axis) |
| => Reflect - reflect at boundary | | => Reflect - reflect at boundary |
| => Truncate - truncate at boundary | | => Truncate - truncate at boundary |
| => Replicate - repeat boundary pixel values | | => Replicate - repeat boundary pixel values |
| | |
| =for bad | | =for bad |
| | |
| Unlike the FFT routines, conv2d is able to process bad values. | | Unlike the FFT routines, conv2d is able to process bad values. |
| | |
| =cut | | =cut |
|
| #line 122 "Image2D.pm" | | |
| | |
| #line 959 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| sub PDL::conv2d { | | sub PDL::conv2d { |
| my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; | | my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; |
| die 'Usage: conv2d( a(m,n), kern(p,q), [o]b(m,n), {Options} )' | | die 'Usage: conv2d( a(m,n), kern(p,q), [o]b(m,n), {Options} )' |
| if $#_<1 || $#_>2; | | if $#_<1 || $#_>2; |
| my($x,$kern) = @_; | | my($x,$kern) = @_; |
| my $c = $#_ == 2 ? $_[2] : $x->nullcreate; | | my $c = $#_ == 2 ? $_[2] : $x->nullcreate; |
| PDL::_conv2d_int($x,$kern,$c, | | PDL::_conv2d_int($x,$kern,$c, |
| (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} | | (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} |
| ); | | ); |
| return $c; | | return $c; |
| } | | } |
|
| #line 140 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *conv2d = \&PDL::conv2d; | | *conv2d = \&PDL::conv2d; |
|
| #line 147 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 med2d | | =head2 med2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); kern(p,q); [o]b(m,n); double+ [t]tmp(pq); indx [t]mapi(isi
ze); indx [t]mapj(jsize); int opt) | | Signature: (a(m,n); kern(p,q); [o]b(m,n); double+ [t]tmp(pq); indx [t]mapi(isi
ze); indx [t]mapj(jsize); int opt) |
| | |
| =for ref | | =for ref |
| | |
| 2D median-convolution of an array with a kernel (smoothing) | | 2D median-convolution of an array with a kernel (smoothing) |
| | |
| skipping to change at line 157 | | skipping to change at line 137 |
|---|
| => Reflect - reflect at boundary | | => Reflect - reflect at boundary |
| => Truncate - truncate at boundary | | => Truncate - truncate at boundary |
| => Replicate - repeat boundary pixel values | | => Replicate - repeat boundary pixel values |
| | |
| =for bad | | =for bad |
| | |
| Bad values are ignored in the calculation. If all elements within the | | Bad values are ignored in the calculation. If all elements within the |
| kernel are bad, the output is set bad. | | kernel are bad, the output is set bad. |
| | |
| =cut | | =cut |
|
| #line 194 "Image2D.pm" | | |
| | |
| #line 959 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| sub PDL::med2d { | | sub PDL::med2d { |
| my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; | | my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; |
| die 'Usage: med2d( a(m,n), kern(p,q), [o]b(m,n), {Options} )' | | die 'Usage: med2d( a(m,n), kern(p,q), [o]b(m,n), {Options} )' |
| if $#_<1 || $#_>2; | | if $#_<1 || $#_>2; |
| my($x,$kern) = @_; | | my($x,$kern) = @_; |
| croak "med2d: kernel must contain some positive elements.\n" | | croak "med2d: kernel must contain some positive elements.\n" |
| if all( $kern <= 0 ); | | if all( $kern <= 0 ); |
| my $c = $#_ == 2 ? $_[2] : $x->nullcreate; | | my $c = $#_ == 2 ? $_[2] : $x->nullcreate; |
| PDL::_med2d_int($x,$kern,$c, | | PDL::_med2d_int($x,$kern,$c, |
| (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} | | (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} |
| ); | | ); |
| return $c; | | return $c; |
| } | | } |
|
| #line 213 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *med2d = \&PDL::med2d; | | *med2d = \&PDL::med2d; |
|
| #line 220 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 med2df | | =head2 med2df |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); [o]b(m,n); indx [t]mapi(isize); indx [t]mapj(jsize); int p
_size=>p; int q_size=>q; int opt) | | Signature: (a(m,n); [o]b(m,n); indx [t]mapi(isize); indx [t]mapj(jsize); int p
_size=>p; int q_size=>q; int opt) |
| | |
| =for ref | | =for ref |
| | |
| 2D median-convolution of an array in a pxq window (smoothing) | | 2D median-convolution of an array in a pxq window (smoothing) |
| | |
| skipping to change at line 220 | | skipping to change at line 191 |
|---|
| => Reflect - reflect at boundary | | => Reflect - reflect at boundary |
| => Truncate - truncate at boundary | | => Truncate - truncate at boundary |
| => Replicate - repeat boundary pixel values | | => Replicate - repeat boundary pixel values |
| | |
| =for bad | | =for bad |
| | |
| med2df does not process bad values. | | med2df does not process bad values. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 269 "Image2D.pm" | | |
| | |
| #line 959 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| sub PDL::med2df { | | sub PDL::med2df { |
| my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; | | my $opt; $opt = pop @_ if ref($_[$#_]) eq 'HASH'; |
| die 'Usage: med2df( a(m,n), [o]b(m,n), p, q, {Options} )' | | die 'Usage: med2df( a(m,n), [o]b(m,n), p, q, {Options} )' |
| if $#_<2 || $#_>3; | | if $#_<2 || $#_>3; |
| my($x,$p,$q) = @_; | | my($x,$p,$q) = @_; |
| croak "med2df: kernel must contain some positive elements.\n" | | croak "med2df: kernel must contain some positive elements.\n" |
| if $p == 0 && $q == 0; | | if $p == 0 && $q == 0; |
| my $c = $#_ == 3 ? $_[3] : $x->nullcreate; | | my $c = $#_ == 3 ? $_[3] : $x->nullcreate; |
| &PDL::_med2df_int($x,$c,$p,$q, | | &PDL::_med2df_int($x,$c,$p,$q, |
| (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} | | (!($opt && exists $$opt{Boundary}))?0:$boundary2value{$$opt{Boundary}} |
| ); | | ); |
| return $c; | | return $c; |
| } | | } |
|
| #line 288 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *med2df = \&PDL::med2df; | | *med2df = \&PDL::med2df; |
|
| #line 295 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 box2d | | =head2 box2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(n,m); [o] b(n,m); int wx; int wy; int edgezero) | | Signature: (a(n,m); [o] b(n,m); int wx; int wy; int edgezero) |
| | |
| =for ref | | =for ref |
| | |
| fast 2D boxcar average | | fast 2D boxcar average |
| | |
| skipping to change at line 277 | | skipping to change at line 239 |
|---|
| better filters are around (e.g., use L</conv2d> with the appropriate | | better filters are around (e.g., use L</conv2d> with the appropriate |
| kernel). On the other hand it is fast and computational cost grows only | | kernel). On the other hand it is fast and computational cost grows only |
| approximately linearly with window size. | | approximately linearly with window size. |
| | |
| =for bad | | =for bad |
| | |
| box2d does not process bad values. | | box2d does not process bad values. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 338 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *box2d = \&PDL::box2d; | | *box2d = \&PDL::box2d; |
|
| #line 345 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 patch2d | | =head2 patch2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); int bad(m,n); [o]b(m,n)) | | Signature: (a(m,n); int bad(m,n); [o]b(m,n)) |
| | |
| =for ref | | =for ref |
| | |
| patch bad pixels out of 2D images using a mask | | patch bad pixels out of 2D images using a mask |
| | |
| skipping to change at line 310 | | skipping to change at line 266 |
|---|
| C<$bad> is a 2D mask array where 1=bad pixel 0=good pixel. | | C<$bad> is a 2D mask array where 1=bad pixel 0=good pixel. |
| Pixels are replaced by the average of their non-bad neighbours; | | Pixels are replaced by the average of their non-bad neighbours; |
| if all neighbours are bad, the original data value is | | if all neighbours are bad, the original data value is |
| copied across. | | copied across. |
| | |
| =for bad | | =for bad |
| | |
| This routine does not handle bad values - use L</patchbad2d> instead | | This routine does not handle bad values - use L</patchbad2d> instead |
| | |
| =cut | | =cut |
|
| #line 380 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *patch2d = \&PDL::patch2d; | | *patch2d = \&PDL::patch2d; |
|
| #line 387 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 patchbad2d | | =head2 patchbad2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); [o]b(m,n)) | | Signature: (a(m,n); [o]b(m,n)) |
| | |
| =for ref | | =for ref |
| | |
| patch bad pixels out of 2D images containing bad values | | patch bad pixels out of 2D images containing bad values |
| | |
| skipping to change at line 344 | | skipping to change at line 294 |
|---|
| if all neighbours are bad, the output is set bad. | | if all neighbours are bad, the output is set bad. |
| If the input ndarray contains I<no> bad values, then a straight copy | | If the input ndarray contains I<no> bad values, then a straight copy |
| is performed (see L</patch2d>). | | is performed (see L</patch2d>). |
| | |
| =for bad | | =for bad |
| | |
| patchbad2d handles bad values. The output ndarray I<may> contain | | patchbad2d handles bad values. The output ndarray I<may> contain |
| bad values, depending on the pattern of bad values in the input ndarray. | | bad values, depending on the pattern of bad values in the input ndarray. |
| | |
| =cut | | =cut |
|
| #line 421 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *patchbad2d = \&PDL::patchbad2d; | | *patchbad2d = \&PDL::patchbad2d; |
|
| #line 428 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 max2d_ind | | =head2 max2d_ind |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); [o]val(); int [o]x(); int[o]y()) | | Signature: (a(m,n); [o]val(); int [o]x(); int[o]y()) |
| | |
| =for ref | | =for ref |
| | |
| Return value/position of maximum value in 2D image | | Return value/position of maximum value in 2D image |
| | |
| Contributed by Tim Jenness | | Contributed by Tim Jenness |
| | |
| =for bad | | =for bad |
| | |
| Bad values are excluded from the search. If all pixels | | Bad values are excluded from the search. If all pixels |
| are bad then the output is set bad. | | are bad then the output is set bad. |
| | |
| =cut | | =cut |
|
| #line 459 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *max2d_ind = \&PDL::max2d_ind; | | *max2d_ind = \&PDL::max2d_ind; |
|
| #line 466 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 centroid2d | | =head2 centroid2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (im(m,n); x(); y(); box(); [o]xcen(); [o]ycen()) | | Signature: (im(m,n); x(); y(); box(); [o]xcen(); [o]ycen()) |
| | |
| =for ref | | =for ref |
| | |
| Refine a list of object positions in 2D image by centroiding in a box | | Refine a list of object positions in 2D image by centroiding in a box |
| | |
| skipping to change at line 400 | | skipping to change at line 338 |
|---|
| C<$box> is the full-width of the box, i.e. the window | | C<$box> is the full-width of the box, i.e. the window |
| is C<+/- $box/2>. | | is C<+/- $box/2>. |
| | |
| =for bad | | =for bad |
| | |
| Bad pixels are excluded from the centroid calculation. If all elements are | | Bad pixels are excluded from the centroid calculation. If all elements are |
| bad (or the pixel sum is 0 - but why would you be centroiding | | bad (or the pixel sum is 0 - but why would you be centroiding |
| something with negatives in...) then the output values are set bad. | | something with negatives in...) then the output values are set bad. |
| | |
| =cut | | =cut |
|
| #line 497 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *centroid2d = \&PDL::centroid2d; | | *centroid2d = \&PDL::centroid2d; |
|
| #line 504 "Image2D.pm" | | |
| | |
| #line 798 "image2d.pd" | | #line 798 "image2d.pd" |
| | |
| =head2 crop | | =head2 crop |
| | |
| =for ref | | =for ref |
| | |
| Return bounding box of given mask in an C<indx> ndarray, so it can broadcast. | | Return bounding box of given mask in an C<indx> ndarray, so it can broadcast. |
| Use other operations (such as L<PDL::Bad/isgood>, or | | Use other operations (such as L<PDL::Bad/isgood>, or |
| L<PDL::Primitive/eqvec> with a colour vector) to create a mask suitable | | L<PDL::Primitive/eqvec> with a colour vector) to create a mask suitable |
| | |
| skipping to change at line 435 | | skipping to change at line 369 |
|---|
| sub PDL::crop { | | sub PDL::crop { |
| my ($mask) = @_; | | my ($mask) = @_; |
| $mask->xchg(0,1)->orover->_which_int(my $out = null, null); | | $mask->xchg(0,1)->orover->_which_int(my $out = null, null); |
| $out->badflag(1); $out->badvalue(-1); | | $out->badflag(1); $out->badvalue(-1); |
| my ($x1, $x2) = $out->minmaximum; | | my ($x1, $x2) = $out->minmaximum; |
| $mask->orover->_which_int($out = null, null); | | $mask->orover->_which_int($out = null, null); |
| $out->badflag(1); $out->badvalue(-1); | | $out->badflag(1); $out->badvalue(-1); |
| my ($y1, $y2) = $out->minmaximum; | | my ($y1, $y2) = $out->minmaximum; |
| $x1->cat($x2, $y1, $y2)->mv(-1,0); | | $x1->cat($x2, $y1, $y2)->mv(-1,0); |
| } | | } |
|
| #line 536 "Image2D.pm" | | |
| | |
| #line 828 "image2d.pd" | | #line 828 "image2d.pd" |
|
| | | |
| =head2 cc8compt | | =head2 cc8compt |
| | |
| =for ref | | =for ref |
| | |
| Connected 8-component labeling of a binary image. | | Connected 8-component labeling of a binary image. |
| | |
| Connected 8-component labeling of 0,1 image - i.e. find separate | | Connected 8-component labeling of 0,1 image - i.e. find separate |
| segmented objects and fill object pixels with object number. | | segmented objects and fill object pixels with object number. |
| 8-component labeling includes all neighboring pixels. | | 8-component labeling includes all neighboring pixels. |
| This is just a front-end to ccNcompt. See also L</cc4compt>. | | This is just a front-end to ccNcompt. See also L</cc4compt>. |
| | |
| skipping to change at line 480 | | skipping to change at line 412 |
|---|
| | |
| sub PDL::cc8compt{ | | sub PDL::cc8compt{ |
| return ccNcompt(shift,8); | | return ccNcompt(shift,8); |
| } | | } |
| *cc8compt = \&PDL::cc8compt; | | *cc8compt = \&PDL::cc8compt; |
| | |
| sub PDL::cc4compt{ | | sub PDL::cc4compt{ |
| return ccNcompt(shift,4); | | return ccNcompt(shift,4); |
| } | | } |
| *cc4compt = \&PDL::cc4compt; | | *cc4compt = \&PDL::cc4compt; |
|
| #line 584 "Image2D.pm" | | #line 508 "Image2D.pm" |
| | | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 ccNcompt | | =head2 ccNcompt |
| | |
| =for sig | | =for sig |
| | |
| Signature: (a(m,n); int+ [o]b(m,n); int con) | | Signature: (a(m,n); int+ [o]b(m,n); int con) |
| | |
| =for ref | | =for ref |
| | |
| Connected component labeling of a binary image. | | Connected component labeling of a binary image. |
| | |
| skipping to change at line 514 | | skipping to change at line 444 |
|---|
| $segmented2 = ccNcompt( $image > $threshold, 8); | | $segmented2 = ccNcompt( $image > $threshold, 8); |
| | |
| where the second parameter specifies the connectivity (4 or 8) of the labeling. | | where the second parameter specifies the connectivity (4 or 8) of the labeling. |
| | |
| =for bad | | =for bad |
| | |
| ccNcompt ignores the bad-value flag of the input ndarrays. | | ccNcompt ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 627 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *ccNcompt = \&PDL::ccNcompt; | | *ccNcompt = \&PDL::ccNcompt; |
|
| #line 634 "Image2D.pm" | | |
| | |
| #line 997 "image2d.pd" | | #line 997 "image2d.pd" |
| | |
| =head2 polyfill | | =head2 polyfill |
| | |
| =for ref | | =for ref |
| | |
| fill the area of the given polygon with the given colour. | | fill the area of the given polygon with the given colour. |
| | |
| This function works inplace, i.e. modifies C<im>. | | This function works inplace, i.e. modifies C<im>. |
| | |
| skipping to change at line 575 | | skipping to change at line 501 |
|---|
| } | | } |
| else | | else |
| { | | { |
| PDL::polyfill_pp($im,$ps,$colour); | | PDL::polyfill_pp($im,$ps,$colour); |
| } | | } |
| return $im; | | return $im; |
| | |
| } | | } |
| | |
| *polyfill = \&PDL::polyfill; | | *polyfill = \&PDL::polyfill; |
|
| #line 692 "Image2D.pm" | | |
| | |
| #line 1054 "image2d.pd" | | #line 1054 "image2d.pd" |
|
| | | |
| =head2 pnpoly | | =head2 pnpoly |
| | |
| =for ref | | =for ref |
| | |
| 'points in a polygon' selection from a 2-D ndarray | | 'points in a polygon' selection from a 2-D ndarray |
| | |
| =for usage | | =for usage |
| | |
| $mask = $img->pnpoly($ps); | | $mask = $img->pnpoly($ps); |
| | |
| | |
| skipping to change at line 656 | | skipping to change at line 580 |
|---|
| my $testy = PDL::Core::topdl($ty)->dummy(0); | | my $testy = PDL::Core::topdl($ty)->dummy(0); |
| my $vertxj = PDL::Core::topdl($vertx)->rotate(1); | | my $vertxj = PDL::Core::topdl($vertx)->rotate(1); |
| my $vertyj = PDL::Core::topdl($verty)->rotate(1); | | my $vertyj = PDL::Core::topdl($verty)->rotate(1); |
| my $mask = ( ($verty>$testy) == ($vertyj>$testy) ); | | my $mask = ( ($verty>$testy) == ($vertyj>$testy) ); |
| my $c = sumover( ! $mask & ( $testx < ($vertxj-$vertx) * ($testy-$verty) | | my $c = sumover( ! $mask & ( $testx < ($vertxj-$vertx) * ($testy-$verty) |
| / ($vertyj-$verty+$mask) + $vertx) ) % 2; | | / ($vertyj-$verty+$mask) + $vertx) ) % 2; |
| return $c; | | return $c; |
| } | | } |
| | |
| *pnpoly = \&PDL::pnpoly; | | *pnpoly = \&PDL::pnpoly; |
|
| #line 776 "Image2D.pm" | | |
| | |
| #line 1137 "image2d.pd" | | #line 1137 "image2d.pd" |
|
| | | |
| =head2 polyfillv | | =head2 polyfillv |
| | |
| =for ref | | =for ref |
| | |
| return the (dataflowed) area of an image described by a polygon | | return the (dataflowed) area of an image described by a polygon |
| | |
| =for usage | | =for usage |
| | |
| polyfillv($im,$ps,[\%options]); | | polyfillv($im,$ps,[\%options]); |
| | |
| | |
| skipping to change at line 710 | | skipping to change at line 632 |
|---|
| my $parsed = PDL::Options->new({'Method' => undef}); | | my $parsed = PDL::Options->new({'Method' => undef}); |
| $parsed->options($opt); | | $parsed->options($opt); |
| return $im->where(PDL::pnpoly_pp($im, $ps)) if $parsed->current->
{'Method'} eq 'pnpoly'; | | return $im->where(PDL::pnpoly_pp($im, $ps)) if $parsed->current->
{'Method'} eq 'pnpoly'; |
| } | | } |
| | |
| my $msk = zeroes(long,$im->dims); | | my $msk = zeroes(long,$im->dims); |
| PDL::polyfill_pp($msk, $ps, 1); | | PDL::polyfill_pp($msk, $ps, 1); |
| return $im->where($msk); | | return $im->where($msk); |
| } | | } |
| *polyfillv = \&PDL::polyfillv; | | *polyfillv = \&PDL::polyfillv; |
|
| #line 833 "Image2D.pm" | | #line 735 "Image2D.pm" |
| | | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 rot2d | | =head2 rot2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (im(m,n); float angle(); bg(); int aa(); [o] om(p,q)) | | Signature: (im(m,n); float angle(); bg(); int aa(); [o] om(p,q)) |
| | |
| =for ref | | =for ref |
| | |
| rotate an image by given C<angle> | | rotate an image by given C<angle> |
| | |
| skipping to change at line 752 | | skipping to change at line 672 |
|---|
| L<PDL::Transform> offers a more general interface to | | L<PDL::Transform> offers a more general interface to |
| distortions, including rotation, with various types of sampling; but | | distortions, including rotation, with various types of sampling; but |
| rot2d is faster. | | rot2d is faster. |
| | |
| =for bad | | =for bad |
| | |
| rot2d ignores the bad-value flag of the input ndarrays. | | rot2d ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 883 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *rot2d = \&PDL::rot2d; | | *rot2d = \&PDL::rot2d; |
|
| #line 890 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 bilin2d | | =head2 bilin2d |
| | |
| =for sig | | =for sig |
| | |
|
| Signature: (Int(n,m); O(q,p)) | | Signature: (Int(n,m); [io] O(q,p)) |
| | |
| =for ref | | =for ref |
| | |
| Bilinearly maps the first ndarray in the second. The | | Bilinearly maps the first ndarray in the second. The |
| interpolated values are actually added to the second | | interpolated values are actually added to the second |
| ndarray which is supposed to be larger than the first one. | | ndarray which is supposed to be larger than the first one. |
| | |
| =for bad | | =for bad |
| | |
| bilin2d ignores the bad-value flag of the input ndarrays. | | bilin2d ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 920 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *bilin2d = \&PDL::bilin2d; | | *bilin2d = \&PDL::bilin2d; |
|
| #line 927 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 rescale2d | | =head2 rescale2d |
| | |
| =for sig | | =for sig |
| | |
|
| Signature: (Int(m,n); O(p,q)) | | Signature: (Int(m,n); [io] O(p,q)) |
| | |
| =for ref | | =for ref |
| | |
| The first ndarray is rescaled to the dimensions of the second | | The first ndarray is rescaled to the dimensions of the second |
| (expanding or meaning values as needed) and then added to it in place. | | (expanding or meaning values as needed) and then added to it in place. |
| Nothing useful is returned. | | Nothing useful is returned. |
| | |
| If you want photometric accuracy or automatic FITS header metadata | | If you want photometric accuracy or automatic FITS header metadata |
| tracking, consider using L<PDL::Transform::map|PDL::Transform/map> | | tracking, consider using L<PDL::Transform::map|PDL::Transform/map> |
| instead: it does these things, at some speed penalty compared to | | instead: it does these things, at some speed penalty compared to |
| rescale2d. | | rescale2d. |
| | |
| =for bad | | =for bad |
| | |
| rescale2d ignores the bad-value flag of the input ndarrays. | | rescale2d ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 962 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *rescale2d = \&PDL::rescale2d; | | *rescale2d = \&PDL::rescale2d; |
|
| #line 969 "Image2D.pm" | | |
| | |
|
| #line 1454 "image2d.pd" | | #line 1442 "image2d.pd" |
| | |
| =head2 fitwarp2d | | =head2 fitwarp2d |
| | |
| =for ref | | =for ref |
| | |
| Find the best-fit 2D polynomial to describe | | Find the best-fit 2D polynomial to describe |
| a coordinate transformation. | | a coordinate transformation. |
| | |
| =for usage | | =for usage |
| | |
| | |
| skipping to change at line 1241 | | skipping to change at line 1145 |
|---|
| # | | # |
| my $mat = _mkbasis( ones(byte,$nf,$nf), $npts, $u, $v ); | | my $mat = _mkbasis( ones(byte,$nf,$nf), $npts, $u, $v ); |
| | |
| my $x = reshape( $mat x $px->clump(-1)->transpose(), $npts ); | | my $x = reshape( $mat x $px->clump(-1)->transpose(), $npts ); |
| my $y = reshape( $mat x $py->clump(-1)->transpose(), $npts ); | | my $y = reshape( $mat x $py->clump(-1)->transpose(), $npts ); |
| return ( $x, $y ); | | return ( $x, $y ); |
| | |
| } # sub: applywarp2d | | } # sub: applywarp2d |
| | |
| *applywarp2d = \&PDL::applywarp2d; | | *applywarp2d = \&PDL::applywarp2d; |
|
| #line 1398 "Image2D.pm" | | #line 1271 "Image2D.pm" |
| | | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 warp2d | | =head2 warp2d |
| | |
| =for sig | | =for sig |
| | |
| Signature: (img(m,n); ldouble px(np,np); ldouble py(np,np); [o] warp(m,n); ldo
uble [t] poly(np); ldouble [t] kernel(ns); char *kernel_type; double noval; int
nsamples => ns) | | Signature: (img(m,n); ldouble px(np,np); ldouble py(np,np); [o] warp(m,n); ldo
uble [t] poly(np); ldouble [t] kernel(ns); char *kernel_type; double noval; int
nsamples => ns) |
| | |
| =for ref | | =for ref |
| | |
| Warp a 2D image given a polynomial describing the I<reverse> mapping. | | Warp a 2D image given a polynomial describing the I<reverse> mapping. |
| | |
| skipping to change at line 1365 | | skipping to change at line 1267 |
|---|
| | |
| C<NOVAL> gives the value used to indicate that a pixel in the | | C<NOVAL> gives the value used to indicate that a pixel in the |
| output image does not map onto one in the input image. | | output image does not map onto one in the input image. |
| | |
| =for bad | | =for bad |
| | |
| warp2d ignores the bad-value flag of the input ndarrays. | | warp2d ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 1528 "Image2D.pm" | | |
| | |
| #line 959 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| # support routine | | # support routine |
| { | | { |
| my %warp2d = map { ($_,1) } qw( tanh sinc sinc2 lanczos hamming hann ); | | my %warp2d = map { ($_,1) } qw( tanh sinc sinc2 lanczos hamming hann ); |
| # note: convert to lower case | | # note: convert to lower case |
| sub _check_kernel ($$) { | | sub _check_kernel ($$) { |
| my $kernel = lc shift; | | my $kernel = lc shift; |
| my $code = shift; | | my $code = shift; |
| barf "Unknown kernel $kernel sent to $code\n" . | | barf "Unknown kernel $kernel sent to $code\n" . |
| "\tmust be one of [" . join(',',keys %warp2d) . "]\n" | | "\tmust be one of [" . join(',',keys %warp2d) . "]\n" |
| | |
| skipping to change at line 1398 | | skipping to change at line 1297 |
|---|
| my $img = shift; | | my $img = shift; |
| my $px = shift; | | my $px = shift; |
| my $py = shift; | | my $py = shift; |
| my $out = $#_ == -1 ? PDL->null() : shift; | | my $out = $#_ == -1 ? PDL->null() : shift; |
| # safety checks | | # safety checks |
| my $copt = $opts->current(); | | my $copt = $opts->current(); |
| my $kernel = _check_kernel( $$copt{KERNEL}, "warp2d" ); | | my $kernel = _check_kernel( $$copt{KERNEL}, "warp2d" ); |
| &PDL::_warp2d_int( $img, $px, $py, $out, $kernel, $$copt{NOVAL}, _get_kernel
_size() ); | | &PDL::_warp2d_int( $img, $px, $py, $out, $kernel, $$copt{NOVAL}, _get_kernel
_size() ); |
| return $out; | | return $out; |
| } | | } |
|
| #line 1563 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *warp2d = \&PDL::warp2d; | | *warp2d = \&PDL::warp2d; |
|
| #line 1570 "Image2D.pm" | | |
| | |
| #line 958 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| =head2 warp2d_kernel | | =head2 warp2d_kernel |
| | |
| =for sig | | =for sig |
| | |
| Signature: ([o] x(n); [o] k(n); ldouble [t] kernel(n); char *name; PDL_Indx ns
ize => n) | | Signature: ([o] x(n); [o] k(n); ldouble [t] kernel(n); char *name; PDL_Indx ns
ize => n) |
| | |
| =for ref | | =for ref |
| | |
| Return the specified kernel, as used by L</warp2d> | | Return the specified kernel, as used by L</warp2d> |
| | |
| skipping to change at line 1434 | | skipping to change at line 1327 |
|---|
| =for example | | =for example |
| | |
| line warp2d_kernel( "hamming" ); | | line warp2d_kernel( "hamming" ); |
| | |
| =for bad | | =for bad |
| | |
| warp2d_kernel ignores the bad-value flag of the input ndarrays. | | warp2d_kernel ignores the bad-value flag of the input ndarrays. |
| It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. | | It will set the bad-value flag of all output ndarrays if the flag is set for any
of the input ndarrays. |
| | |
| =cut | | =cut |
|
| #line 1607 "Image2D.pm" | | |
| | |
| #line 959 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| sub PDL::warp2d_kernel ($) { | | sub PDL::warp2d_kernel ($) { |
| my $kernel = _check_kernel( shift, "warp2d_kernel" ); | | my $kernel = _check_kernel( shift, "warp2d_kernel" ); |
| &PDL::_warp2d_kernel_int( my $x=PDL->null, my $k=PDL->null, $kernel, _get_ke
rnel_size() ); | | &PDL::_warp2d_kernel_int( my $x=PDL->null, my $k=PDL->null, $kernel, _get_ke
rnel_size() ); |
| return ( $x, $k ); | | return ( $x, $k ); |
| } | | } |
|
| #line 1619 "Image2D.pm" | | |
| | |
| #line 960 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | | |
| | |
| *warp2d_kernel = \&PDL::warp2d_kernel; | | *warp2d_kernel = \&PDL::warp2d_kernel; |
|
| #line 1626 "Image2D.pm" | | |
| | |
| #line 31 "image2d.pd" | | #line 31 "image2d.pd" |
| | |
| =head1 AUTHORS | | =head1 AUTHORS |
| | |
| Copyright (C) Karl Glazebrook 1997 with additions by Robin Williams | | Copyright (C) Karl Glazebrook 1997 with additions by Robin Williams |
| (rjrw@ast.leeds.ac.uk), Tim Jenness (timj@jach.hawaii.edu), | | (rjrw@ast.leeds.ac.uk), Tim Jenness (timj@jach.hawaii.edu), |
| and Doug Burke (burke@ifa.hawaii.edu). | | and Doug Burke (burke@ifa.hawaii.edu). |
| | |
| All rights reserved. There is no warranty. You are allowed | | All rights reserved. There is no warranty. You are allowed |
| to redistribute this software / documentation under certain | | to redistribute this software / documentation under certain |
| conditions. For details, see the file COPYING in the PDL | | conditions. For details, see the file COPYING in the PDL |
| distribution. If this file is separated from the PDL distribution, | | distribution. If this file is separated from the PDL distribution, |
| the copyright notice should be included in the file. | | the copyright notice should be included in the file. |
| | |
| =cut | | =cut |
|
| #line 1648 "Image2D.pm" | | #line 1499 "Image2D.pm" |
| | |
| # Exit with OK status | | # Exit with OK status |
| | |
| 1; | | 1; |
| | | |
| End of changes. 48 change blocks. |
|---|
| 122 lines changed or deleted | | 8 lines changed or added |
|---|
"Fossies" - the Fresh Open Source Software Archive 