/*@null@*/ static PyObject*
Wcs_pix2foc(
Wcs* self,
PyObject* args,
PyObject* kwds) {
PyObject* pixcrd_obj = NULL;
int origin = 1;
PyArrayObject* pixcrd = NULL;
PyArrayObject* foccrd = NULL;
int status = -1;
const char* keywords[] = {
"pixcrd", "origin", NULL
};
if (!PyArg_ParseTupleAndKeywords(args, kwds, "Oi:pix2foc", (char **)keywords,
&pixcrd_obj, &origin)) {
return NULL;
}
pixcrd = (PyArrayObject*)PyArray_ContiguousFromAny(pixcrd_obj, PyArray_DOUBLE, 2, 2);
if (pixcrd == NULL) {
return NULL;
}
if (PyArray_DIM(pixcrd, 1) != NAXES) {
PyErr_SetString(PyExc_ValueError, "Pixel array must be an Nx2 array");
goto _exit;
}
foccrd = (PyArrayObject*)PyArray_SimpleNew(2, PyArray_DIMS(pixcrd), PyArray_DOUBLE);
if (foccrd == NULL) {
goto _exit;
}
Py_BEGIN_ALLOW_THREADS
preoffset_array(pixcrd, origin);
status = pipeline_pix2foc(&self->x,
(unsigned int)PyArray_DIM(pixcrd, 0),
(unsigned int)PyArray_DIM(pixcrd, 1),
(double*)PyArray_DATA(pixcrd),
(double*)PyArray_DATA(foccrd));
unoffset_array(pixcrd, origin);
unoffset_array(foccrd, origin);
Py_END_ALLOW_THREADS
_exit:
Py_XDECREF(pixcrd);
if (status == 0) {
return (PyObject*)foccrd;
} else {
Py_XDECREF(foccrd);
if (status == -1) {
/* Exception already set */
return NULL;
} else {
wcserr_to_python_exc(self->x.err);
return NULL;
}
}
}
/*@null@*/ static PyObject*
Wcs_all_pix2world(
Wcs* self,
PyObject* args,
PyObject* kwds) {
int naxis = 2;
PyObject* pixcrd_obj = NULL;
int origin = 1;
PyArrayObject* pixcrd = NULL;
PyArrayObject* world = NULL;
int status = -1;
const char* keywords[] = {
"pixcrd", "origin", NULL
};
if (!PyArg_ParseTupleAndKeywords(
args, kwds, "Oi:all_pix2world", (char **)keywords,
&pixcrd_obj, &origin)) {
return NULL;
}
naxis = self->x.wcs->naxis;
pixcrd = (PyArrayObject*)PyArray_ContiguousFromAny(pixcrd_obj, PyArray_DOUBLE, 2, 2);
if (pixcrd == NULL) {
return NULL;
}
if (PyArray_DIM(pixcrd, 1) < naxis) {
PyErr_Format(
PyExc_RuntimeError,
"Input array must be 2-dimensional, where the second dimension >= %d",
naxis);
goto exit;
}
world = (PyArrayObject*)PyArray_SimpleNew(2, PyArray_DIMS(pixcrd), PyArray_DOUBLE);
if (world == NULL) {
goto exit;
}
/* Make the call */
Py_BEGIN_ALLOW_THREADS
preoffset_array(pixcrd, origin);
wcsprm_python2c(self->x.wcs);
status = pipeline_all_pixel2world(&self->x,
(unsigned int)PyArray_DIM(pixcrd, 0),
(unsigned int)PyArray_DIM(pixcrd, 1),
(double*)PyArray_DATA(pixcrd),
(double*)PyArray_DATA(world));
wcsprm_c2python(self->x.wcs);
unoffset_array(pixcrd, origin);
Py_END_ALLOW_THREADS
/* unoffset_array(world, origin); */
exit:
Py_XDECREF(pixcrd);
if (status == 0 || status == 8) {
return (PyObject*)world;
} else if (status == -1) {
PyErr_SetString(
PyExc_ValueError,
"Wrong number of dimensions in input array. Expected 2.");
return NULL;
} else {
Py_DECREF(world);
if (status == -1) {
/* exception already set */
return NULL;
} else {
wcserr_to_python_exc(self->x.err);
return NULL;
}
}
}
/*@null@*/ static PyObject*
Wcs_det2im(
Wcs* self,
PyObject* args,
PyObject* kwds) {
PyObject* detcrd_obj = NULL;
int origin = 1;
PyArrayObject* detcrd = NULL;
PyArrayObject* imcrd = NULL;
int status = -1;
const char* keywords[] = {
"detcrd", "origin", NULL
};
if (!PyArg_ParseTupleAndKeywords(args, kwds, "Oi:det2im", (char **)keywords,
&detcrd_obj, &origin)) {
return NULL;
}
if (self->x.det2im[0] == NULL && self->x.det2im[1] == NULL) {
Py_INCREF(detcrd_obj);
return detcrd_obj;
}
detcrd = (PyArrayObject*)PyArray_ContiguousFromAny(detcrd_obj, PyArray_DOUBLE, 2, 2);
if (detcrd == NULL) {
return NULL;
}
if (PyArray_DIM(detcrd, 1) != NAXES) {
PyErr_SetString(PyExc_ValueError, "Pixel array must be an Nx2 array");
goto exit;
}
imcrd = (PyArrayObject*)PyArray_SimpleNew(2, PyArray_DIMS(detcrd), PyArray_DOUBLE);
if (imcrd == NULL) {
status = 2;
goto exit;
}
Py_BEGIN_ALLOW_THREADS
preoffset_array(detcrd, origin);
status = p4_pix2foc(2, (void *)self->x.det2im,
(unsigned int)PyArray_DIM(detcrd, 0),
(double*)PyArray_DATA(detcrd),
(double*)PyArray_DATA(imcrd));
unoffset_array(detcrd, origin);
unoffset_array(imcrd, origin);
Py_END_ALLOW_THREADS
exit:
Py_XDECREF(detcrd);
if (status == 0) {
return (PyObject*)imcrd;
} else {
Py_XDECREF(imcrd);
if (status == -1) {
/* Exception already set */
return NULL;
} else {
PyErr_SetString(PyExc_MemoryError, "NULL pointer passed");
return NULL;
}
}
}
/*@null@*/ static PyObject*
PySip_foc2pix(
PySip* self,
PyObject* args,
PyObject* kwds) {
PyObject* foccrd_obj = NULL;
int origin = 1;
PyArrayObject* foccrd = NULL;
PyArrayObject* pixcrd = NULL;
int status = -1;
const char* keywords[] = {
"foccrd", "origin", NULL };
if (!PyArg_ParseTupleAndKeywords(args, kwds, "Oi:foc2pix", (char **)keywords,
&foccrd_obj, &origin)) {
return NULL;
}
if (self->x.ap == NULL || self->x.bp == NULL) {
PyErr_SetString(
PyExc_ValueError,
"SIP object does not have coefficients for foc2pix transformation (AP and BP)");
return NULL;
}
foccrd = (PyArrayObject*)PyArray_ContiguousFromAny(foccrd_obj, PyArray_DOUBLE, 2, 2);
if (foccrd == NULL) {
goto exit;
}
if (PyArray_DIM(foccrd, 1) != 2) {
PyErr_SetString(PyExc_ValueError, "Pixel array must be an Nx2 array");
goto exit;
}
pixcrd = (PyArrayObject*)PyArray_SimpleNew(2, PyArray_DIMS(foccrd),
PyArray_DOUBLE);
if (pixcrd == NULL) {
status = 2;
goto exit;
}
Py_BEGIN_ALLOW_THREADS
preoffset_array(foccrd, origin);
status = sip_foc2pix(&self->x,
(unsigned int)PyArray_DIM(pixcrd, 1),
(unsigned int)PyArray_DIM(pixcrd, 0),
(double*)PyArray_DATA(foccrd),
(double*)PyArray_DATA(pixcrd));
unoffset_array(foccrd, origin);
unoffset_array(pixcrd, origin);
Py_END_ALLOW_THREADS
exit:
Py_XDECREF(foccrd);
if (status == 0) {
return (PyObject*)pixcrd;
} else {
Py_XDECREF(pixcrd);
if (status == -1) {
/* Exception already set */
return NULL;
} else {
wcserr_to_python_exc(self->x.err);
return NULL;
}
}
}
