Image2D.pm (PDL-2.080) | : | Image2D.pm (PDL-2.081) | ||
---|---|---|---|---|

skipping to change at line 17 | skipping to change at line 17 | |||

our %EXPORT_TAGS = (Func=>\@EXPORT_OK); | our %EXPORT_TAGS = (Func=>\@EXPORT_OK); | |||

use PDL::Core; | use PDL::Core; | |||

use PDL::Exporter; | use PDL::Exporter; | |||

use DynaLoader; | use DynaLoader; | |||

our @ISA = ( 'PDL::Exporter','DynaLoader' ); | our @ISA = ( 'PDL::Exporter','DynaLoader' ); | |||

push @PDL::Core::PP, __PACKAGE__; | push @PDL::Core::PP, __PACKAGE__; | |||

bootstrap PDL::Image2D ; | bootstrap PDL::Image2D ; | |||

#line 4 "image2d.pd" | #line 5 "image2d.pd" | |||

use strict; | use strict; | |||

use warnings; | use warnings; | |||

=head1 NAME | =head1 NAME | |||

PDL::Image2D - Miscellaneous 2D image processing functions | PDL::Image2D - Miscellaneous 2D image processing functions | |||

=head1 DESCRIPTION | =head1 DESCRIPTION | |||

skipping to change at line 46 | skipping to change at line 46 | |||

use PDL; # ensure qsort routine available | use PDL; # ensure qsort routine available | |||

use PDL::Math; | use PDL::Math; | |||

use Carp; | use Carp; | |||

#line 48 "Image2D.pm" | #line 48 "Image2D.pm" | |||

=head1 FUNCTIONS | =head1 FUNCTIONS | |||

=cut | =cut | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

#line 63 "Image2D.pm" | #line 63 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

#line 68 "Image2D.pm" | #line 68 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

#line 73 "Image2D.pm" | #line 73 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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); int opt) | Signature: (a(m,n); kern(p,q); [o]b(m,n); 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 96 | |||

=> 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 120 "Image2D.pm" | #line 120 "Image2D.pm" | |||

#line 949 "../../blib/lib/PDL/PP.pm" | #line 950 "/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, | |||

(!(defined $opt && exists $$opt{Boundary}))?0: | (!(defined $opt && exists $$opt{Boundary}))?0: | |||

(($$opt{Boundary} eq "Reflect") + | (($$opt{Boundary} eq "Reflect") + | |||

2*($$opt{Boundary} eq "Truncate") + | 2*($$opt{Boundary} eq "Truncate") + | |||

3*($$opt{Boundary} eq "Replicate"))); | 3*($$opt{Boundary} eq "Replicate"))); | |||

return $c; | return $c; | |||

} | } | |||

#line 140 "Image2D.pm" | #line 140 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*conv2d = \&PDL::conv2d; | *conv2d = \&PDL::conv2d; | |||

#line 147 "Image2D.pm" | #line 147 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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); int opt) | Signature: (a(m,n); kern(p,q); [o]b(m,n); double+ [t]tmp(pq); 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) | |||

Note: only points in the kernel E<gt>0 are included in the median, other | Note: only points in the kernel E<gt>0 are included in the median, other | |||

points are weighted by the kernel value (medianing lots of zeroes | points are weighted by the kernel value (medianing lots of zeroes | |||

is rather pointless) | is rather pointless) | |||

=for usage | =for usage | |||

skipping to change at line 159 | skipping to change at line 159 | |||

=> 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 194 "Image2D.pm" | |||

#line 949 "../../blib/lib/PDL/PP.pm" | #line 950 "/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, | |||

(!(defined $opt && exists $$opt{Boundary}))?0: | (!(defined $opt && exists $$opt{Boundary}))?0: | |||

(($$opt{Boundary} eq "Reflect") + | (($$opt{Boundary} eq "Reflect") + | |||

2*($$opt{Boundary} eq "Truncate") + | 2*($$opt{Boundary} eq "Truncate") + | |||

3*($$opt{Boundary} eq "Replicate"))); | 3*($$opt{Boundary} eq "Replicate"))); | |||

return $c; | return $c; | |||

} | } | |||

#line 215 "Image2D.pm" | #line 215 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*med2d = \&PDL::med2d; | *med2d = \&PDL::med2d; | |||

#line 222 "Image2D.pm" | #line 222 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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); int __p_size; int __q_size; int opt) | Signature: (a(m,n); [o]b(m,n); int __p_size; int __q_size; 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 224 | skipping to change at line 224 | |||

=> 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 271 "Image2D.pm" | #line 271 "Image2D.pm" | |||

#line 949 "../../blib/lib/PDL/PP.pm" | #line 950 "/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, | |||

(!(defined $opt && exists $$opt{Boundary}))?0: | (!(defined $opt && exists $$opt{Boundary}))?0: | |||

(($$opt{Boundary} eq "Reflect") + | (($$opt{Boundary} eq "Reflect") + | |||

2*($$opt{Boundary} eq "Truncate") + | 2*($$opt{Boundary} eq "Truncate") + | |||

3*($$opt{Boundary} eq "Replicate"))); | 3*($$opt{Boundary} eq "Replicate"))); | |||

return $c; | return $c; | |||

} | } | |||

#line 292 "Image2D.pm" | #line 292 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*med2df = \&PDL::med2df; | *med2df = \&PDL::med2df; | |||

#line 299 "Image2D.pm" | #line 299 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 283 | skipping to change at line 283 | |||

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. | ||||

=cut | =cut | |||

#line 342 "Image2D.pm" | #line 342 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*box2d = \&PDL::box2d; | *box2d = \&PDL::box2d; | |||

#line 349 "Image2D.pm" | #line 349 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 316 | skipping to change at line 316 | |||

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 384 "Image2D.pm" | #line 384 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*patch2d = \&PDL::patch2d; | *patch2d = \&PDL::patch2d; | |||

#line 391 "Image2D.pm" | #line 391 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 350 | skipping to change at line 350 | |||

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 425 "Image2D.pm" | #line 425 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*patchbad2d = \&PDL::patchbad2d; | *patchbad2d = \&PDL::patchbad2d; | |||

#line 432 "Image2D.pm" | #line 432 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 | |||

skipping to change at line 377 | skipping to change at line 377 | |||

Contributed by Tim Jeness | Contributed by Tim Jeness | |||

=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 463 "Image2D.pm" | #line 463 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*max2d_ind = \&PDL::max2d_ind; | *max2d_ind = \&PDL::max2d_ind; | |||

#line 470 "Image2D.pm" | #line 470 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 406 | skipping to change at line 406 | |||

=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 501 "Image2D.pm" | #line 501 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*centroid2d = \&PDL::centroid2d; | *centroid2d = \&PDL::centroid2d; | |||

#line 508 "Image2D.pm" | #line 508 "Image2D.pm" | |||

#line 797 "image2d.pd" | #line 803 "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 | |||

for your application. | for your application. | |||

skipping to change at line 441 | skipping to change at line 441 | |||

$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 540 "Image2D.pm" | #line 540 "Image2D.pm" | |||

#line 827 "image2d.pd" | #line 833 "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. | |||

skipping to change at line 486 | skipping to change at line 486 | |||

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 588 "Image2D.pm" | #line 588 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 520 | skipping to change at line 520 | |||

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. | |||

=cut | =cut | |||

#line 631 "Image2D.pm" | #line 631 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*ccNcompt = \&PDL::ccNcompt; | *ccNcompt = \&PDL::ccNcompt; | |||

#line 638 "Image2D.pm" | #line 638 "Image2D.pm" | |||

#line 996 "image2d.pd" | #line 1002 "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>. | |||

=for usage | =for usage | |||

skipping to change at line 581 | skipping to change at line 581 | |||

{ | { | |||

PDL::polyfill_pp($im,$ps,$colour); | PDL::polyfill_pp($im,$ps,$colour); | |||

} | } | |||

return $im; | return $im; | |||

} | } | |||

*polyfill = \&PDL::polyfill; | *polyfill = \&PDL::polyfill; | |||

#line 696 "Image2D.pm" | #line 696 "Image2D.pm" | |||

#line 1053 "image2d.pd" | #line 1059 "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 662 | skipping to change at line 662 | |||

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 780 "Image2D.pm" | #line 780 "Image2D.pm" | |||

#line 1136 "image2d.pd" | #line 1142 "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 716 | skipping to change at line 716 | |||

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 837 "Image2D.pm" | #line 837 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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 758 | skipping to change at line 758 | |||

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. | |||

=cut | =cut | |||

#line 887 "Image2D.pm" | #line 887 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*rot2d = \&PDL::rot2d; | *rot2d = \&PDL::rot2d; | |||

#line 894 "Image2D.pm" | #line 894 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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); 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 | |||

skipping to change at line 785 | skipping to change at line 785 | |||

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. | |||

=cut | =cut | |||

#line 924 "Image2D.pm" | #line 924 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*bilin2d = \&PDL::bilin2d; | *bilin2d = \&PDL::bilin2d; | |||

#line 931 "Image2D.pm" | #line 931 "Image2D.pm" | |||

#line 948 "../../blib/lib/PDL/PP.pm" | #line 949 "/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); 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 | |||

skipping to change at line 817 | skipping to change at line 817 | |||

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. | |||

=cut | =cut | |||

#line 966 "Image2D.pm" | #line 966 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*rescale2d = \&PDL::rescale2d; | *rescale2d = \&PDL::rescale2d; | |||

#line 973 "Image2D.pm" | #line 973 "Image2D.pm" | |||

#line 1450 "image2d.pd" | #line 1459 "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 1247 | skipping to change at line 1247 | |||

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 1402 "Image2D.pm" | #line 1402 "Image2D.pm" | |||

#line 1885 "image2d.pd" | #line 1894 "image2d.pd" | |||

=head2 warp2d | =head2 warp2d | |||

=for sig | =for sig | |||

Signature: (img(m,n); double px(np,np); double py(np,np); [o] warp(m,n); { opt ions }) | Signature: (img(m,n); double px(np,np); double py(np,np); [o] warp(m,n); { opt ions }) | |||

=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 1381 | skipping to change at line 1381 | |||

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" | |||

unless exists $warp2d{$kernel}; | unless exists $warp2d{$kernel}; | |||

return $kernel; | return $kernel; | |||

} | } | |||

} | } | |||

#line 1539 "Image2D.pm" | #line 1539 "Image2D.pm" | |||

#line 949 "../../blib/lib/PDL/PP.pm" | #line 950 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

sub PDL::warp2d { | sub PDL::warp2d { | |||

my $opts = PDL::Options->new( { KERNEL => "tanh", NOVAL => 0 } ); | my $opts = PDL::Options->new( { KERNEL => "tanh", NOVAL => 0 } ); | |||

$opts->options( pop(@_) ) if ref($_[$#_]) eq "HASH"; | $opts->options( pop(@_) ) if ref($_[$#_]) eq "HASH"; | |||

die "Usage: warp2d( in(m,n), px(np,np); py(np,np); [o] out(m,n), {Options} ) " | die "Usage: warp2d( in(m,n), px(np,np); py(np,np); [o] out(m,n), {Options} ) " | |||

if $#_<2 || $#_>3; | if $#_<2 || $#_>3; | |||

my $img = shift; | my $img = shift; | |||

my $px = shift; | my $px = shift; | |||

my $py = shift; | my $py = shift; | |||

skipping to change at line 1403 | skipping to change at line 1403 | |||

# 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} ); | &PDL::_warp2d_int( $img, $px, $py, $out, $kernel, $$copt{NOVAL} ); | |||

return $out; | return $out; | |||

} | } | |||

#line 1565 "Image2D.pm" | #line 1565 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*warp2d = \&PDL::warp2d; | *warp2d = \&PDL::warp2d; | |||

#line 1572 "Image2D.pm" | #line 1572 "Image2D.pm" | |||

#line 2199 "image2d.pd" | #line 2199 "image2d.pd" | |||

=head2 warp2d_kernel | =head2 warp2d_kernel | |||

=for ref | =for ref | |||

skipping to change at line 1430 | skipping to change at line 1430 | |||

The valid values for C<$name> are the same as the C<KERNEL> option | The valid values for C<$name> are the same as the C<KERNEL> option | |||

of L<warp2d()|/warp2d>. | of L<warp2d()|/warp2d>. | |||

=for example | =for example | |||

line warp2d_kernel( "hamming" ); | line warp2d_kernel( "hamming" ); | |||

=cut | =cut | |||

#line 1598 "Image2D.pm" | #line 1598 "Image2D.pm" | |||

#line 949 "../../blib/lib/PDL/PP.pm" | #line 950 "/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" ); | |||

my $nelem = _get_kernel_size(); | my $nelem = _get_kernel_size(); | |||

my $x = zeroes( $nelem ); | my $x = zeroes( $nelem ); | |||

my $k = zeroes( $nelem ); | my $k = zeroes( $nelem ); | |||

&PDL::_warp2d_kernel_int( $x, $k, $kernel ); | &PDL::_warp2d_kernel_int( $x, $k, $kernel ); | |||

return ( $x, $k ); | return ( $x, $k ); | |||

# return _get_kernel( $kernel ); | # return _get_kernel( $kernel ); | |||

} | } | |||

*warp2d_kernel = \&PDL::warp2d_kernel; | *warp2d_kernel = \&PDL::warp2d_kernel; | |||

#line 1619 "Image2D.pm" | #line 1619 "Image2D.pm" | |||

#line 950 "../../blib/lib/PDL/PP.pm" | #line 951 "/home/osboxes/pdl-code/blib/lib/PDL/PP.pm" | |||

*warp2d_kernel = \&PDL::warp2d_kernel; | *warp2d_kernel = \&PDL::warp2d_kernel; | |||

#line 1626 "Image2D.pm" | #line 1626 "Image2D.pm" | |||

#line 29 "image2d.pd" | #line 30 "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 Jeness (timj@jach.hawaii.edu), | (rjrw@ast.leeds.ac.uk), Tim Jeness (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 | |||

End of changes. 44 change blocks. | ||||

44 lines changed or deleted | | 44 lines changed or added |