core
/
eccodes
mirror of https://github.com/ecmwf/eccodes.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

GRIB-438: Error encoding values with mars formula (JPEG data)

This commit is contained in:
Shahram Najm 2013-09-30 14:08:41 +01:00
parent 4baa71d7eb
commit 4e97d036ee
2 changed files with 436 additions and 432 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

582
src/grib_accessor_class_data_jpeg2000_packing.c
View File

@ -166,52 +166,51 @@ static int first = 1;
static void init(grib_accessor* a,const long v, grib_arguments* args)
{
const char * user_lib=NULL;
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
const char * user_lib=NULL;
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
self->type_of_compression_used = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->target_compression_ratio = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->ni = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->nj = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->list_defining_points = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->number_of_data_points = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->scanning_mode = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->edition=2;
a->flags |= GRIB_ACCESSOR_FLAG_DATA;
self->type_of_compression_used = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->target_compression_ratio = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->ni = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->nj = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->list_defining_points = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->number_of_data_points = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->scanning_mode = grib_arguments_get_name(a->parent->h,args,self->carg++);
self->edition=2;
a->flags |= GRIB_ACCESSOR_FLAG_DATA;
#if HAVE_LIBJASPER
self->jpeg_lib=JASPER_LIB;
self->jpeg_lib=JASPER_LIB;
#elif HAVE_LIBOPENJPEG
self->jpeg_lib=OPENJPEG_LIB;
self->jpeg_lib=OPENJPEG_LIB;
#endif
if ((user_lib=getenv("GRIB_JPEG"))!=NULL ) {
if (!strcmp(user_lib,"jasper")) {
self->jpeg_lib=JASPER_LIB;
} else if (!strcmp(user_lib,"openjpeg")) {
self->jpeg_lib=OPENJPEG_LIB;
}
}
if (!strcmp(user_lib,"jasper")) {
self->jpeg_lib=JASPER_LIB;
} else if (!strcmp(user_lib,"openjpeg")) {
self->jpeg_lib=OPENJPEG_LIB;
}
}
self->dump_jpg = getenv("GRIB_DUMP_JPG_FILE");
if(first) {
if(self->dump_jpg)
printf("GRIB JPEG dumping to %s\n",self->dump_jpg);
first = 0;
}
self->dump_jpg = getenv("GRIB_DUMP_JPG_FILE");
if(first) {
if(self->dump_jpg)
printf("GRIB JPEG dumping to %s\n",self->dump_jpg);
first = 0;
}
}
static long value_count(grib_accessor* a)
{
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
long n_vals= 0;
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
long n_vals= 0;
if(grib_get_long_internal(a->parent->h,self->number_of_values,&n_vals) != GRIB_SUCCESS)
return 0;
if(grib_get_long_internal(a->parent->h,self->number_of_values,&n_vals) != GRIB_SUCCESS)
return 0;
return n_vals;
return n_vals;
}
#define EXTRA_BUFFER_SIZE 10240
@ -219,332 +218,341 @@ static long value_count(grib_accessor* a)
#ifdef HAVE_JPEG
static int unpack_double(grib_accessor* a, double* val, size_t *len)
{
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
int err = GRIB_SUCCESS;
int i;
size_t buflen = grib_byte_count(a);
int err = GRIB_SUCCESS;
int i;
size_t buflen = grib_byte_count(a);
double bscale = 0;
double dscale = 0;
unsigned char* buf = NULL;
size_t n_vals = 0;
double bscale = 0;
double dscale = 0;
unsigned char* buf = NULL;
size_t n_vals = 0;
long binary_scale_factor = 0;
long decimal_scale_factor = 0;
double reference_value = 0;
long bits_per_value =0;
double units_factor=1.0;
double units_bias=0.0;
long binary_scale_factor = 0;
long decimal_scale_factor = 0;
double reference_value = 0;
long bits_per_value =0;
double units_factor=1.0;
double units_bias=0.0;
n_vals = grib_value_count(a);
n_vals = grib_value_count(a);
if(self->units_factor)
grib_get_double_internal(a->parent->h,self->units_factor,&units_factor);
if(self->units_factor)
grib_get_double_internal(a->parent->h,self->units_factor,&units_factor);
if(self->units_bias)
grib_get_double_internal(a->parent->h,self->units_bias,&units_bias);
if(self->units_bias)
grib_get_double_internal(a->parent->h,self->units_bias,&units_bias);
if((err = grib_get_long_internal(a->parent->h,self->bits_per_value,&bits_per_value)) != GRIB_SUCCESS)
return err;
if((err = grib_get_double_internal(a->parent->h,self->reference_value, &reference_value)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->binary_scale_factor, &binary_scale_factor)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->decimal_scale_factor, &decimal_scale_factor)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->bits_per_value,&bits_per_value)) != GRIB_SUCCESS)
return err;
if((err = grib_get_double_internal(a->parent->h,self->reference_value, &reference_value)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->binary_scale_factor, &binary_scale_factor)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->decimal_scale_factor, &decimal_scale_factor)) != GRIB_SUCCESS)
return err;
self->dirty=0;
self->dirty=0;
bscale = grib_power(binary_scale_factor,2);
dscale = grib_power(-decimal_scale_factor,10);
bscale = grib_power(binary_scale_factor,2);
dscale = grib_power(-decimal_scale_factor,10);
/* TODO: This should be called upstream */
if(*len < n_vals)
return GRIB_ARRAY_TOO_SMALL;
/* TODO: This should be called upstream */
if(*len < n_vals)
return GRIB_ARRAY_TOO_SMALL;
/* Special case */
/* Special case */
if(bits_per_value == 0)
{
for(i = 0; i < n_vals; i++)
val[i] = reference_value;
*len = n_vals;
return GRIB_SUCCESS;
}
buf = (unsigned char*)a->parent->h->buffer->data;
buf += grib_byte_offset(a);
switch (self->jpeg_lib) {
case OPENJPEG_LIB:
if ((err = grib_openjpeg_decode(a->parent->h->context,buf,&buflen,val,&n_vals)) != GRIB_SUCCESS)
return err;
break;
case JASPER_LIB:
if ((err = grib_jasper_decode(a->parent->h->context,buf,&buflen,val,&n_vals)) != GRIB_SUCCESS)
return err;
break;
}
*len = n_vals;
for (i = 0; i < n_vals; i++) {
val[i] = (val[i] * bscale + reference_value) * dscale;
}
if (units_factor != 1.0) {
if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]=val[i]*units_factor+units_bias;
else
for (i=0;i<n_vals;i++) val[i]*=units_factor;
} else if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]+=units_bias;
return err;
}
static int pack_double(grib_accessor* a, const double* cval, size_t *len)
{
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
grib_accessor_class* super = *(a->cclass->super);
size_t n_vals = *len;
int err = 0;
int i;
double reference_value = 0;
long binary_scale_factor = 0;
long bits_per_value = 0;
long decimal_scale_factor = 0;
double decimal = 1;
size_t simple_packing_size = 0;
unsigned char* buf = NULL;
double divisor = 1;
long width;
long height;
long ni;
long nj;
long target_compression_ratio;
long type_of_compression_used;
long scanning_mode;
long list_defining_points;
long number_of_data_points;
int ret =0;
j2k_encode_helper helper;
double units_factor=1.0;
double units_bias=0.0;
double* val=(double*)cval;
self->dirty=1;
if(*len == 0){
grib_buffer_replace(a, NULL, 0, 1, 1);
return GRIB_SUCCESS;
if(bits_per_value == 0)
{
for(i = 0; i < n_vals; i++)
val[i] = reference_value;
*len = n_vals;
return GRIB_SUCCESS;
}
if(self->units_factor &&
(grib_get_double_internal(a->parent->h,self->units_factor,&units_factor)== GRIB_SUCCESS)) {
grib_set_double_internal(a->parent->h,self->units_factor,1.0);
}
buf = (unsigned char*)a->parent->h->buffer->data;
buf += grib_byte_offset(a);
if(self->units_bias &&
(grib_get_double_internal(a->parent->h,self->units_bias,&units_bias)== GRIB_SUCCESS)) {
grib_set_double_internal(a->parent->h,self->units_bias,0.0);
}
switch (self->jpeg_lib) {
case OPENJPEG_LIB:
if ((err = grib_openjpeg_decode(a->parent->h->context,buf,&buflen,val,&n_vals)) != GRIB_SUCCESS)
return err;
break;
case JASPER_LIB:
if ((err = grib_jasper_decode(a->parent->h->context,buf,&buflen,val,&n_vals)) != GRIB_SUCCESS)
return err;
break;
}
if (units_factor != 1.0) {
*len = n_vals;
for (i = 0; i < n_vals; i++) {
val[i] = (val[i] * bscale + reference_value) * dscale;
}
if (units_factor != 1.0) {
if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]=val[i]*units_factor+units_bias;
else
for (i=0;i<n_vals;i++) val[i]*=units_factor;
} else if (units_bias != 0.0)
} else if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]+=units_bias;
ret = super->pack_double(a,val,len);
switch (ret) {
case GRIB_CONSTANT_FIELD:
return err;
}
static int pack_double(grib_accessor* a, const double* cval, size_t *len)
{
grib_accessor_data_jpeg2000_packing *self =(grib_accessor_data_jpeg2000_packing*)a;
grib_accessor_class* super = *(a->cclass->super);
size_t n_vals = *len;
int err = 0;
int i;
double reference_value = 0;
long binary_scale_factor = 0;
long bits_per_value = 0;
long decimal_scale_factor = 0;
double decimal = 1;
size_t simple_packing_size = 0;
unsigned char* buf = NULL;
double divisor = 1;
long width;
long height;
long ni;
long nj;
long target_compression_ratio;
long type_of_compression_used;
long scanning_mode;
long list_defining_points;
long number_of_data_points;
int ret =0;
j2k_encode_helper helper;
double units_factor=1.0;
double units_bias=0.0;
double* val=(double*)cval;
self->dirty=1;
if(*len == 0){
grib_buffer_replace(a, NULL, 0, 1, 1);
return GRIB_SUCCESS;
}
if(self->units_factor &&
(grib_get_double_internal(a->parent->h,self->units_factor,&units_factor)== GRIB_SUCCESS)) {
grib_set_double_internal(a->parent->h,self->units_factor,1.0);
}
if(self->units_bias &&
(grib_get_double_internal(a->parent->h,self->units_bias,&units_bias)== GRIB_SUCCESS)) {
grib_set_double_internal(a->parent->h,self->units_bias,0.0);
}
if (units_factor != 1.0) {
if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]=val[i]*units_factor+units_bias;
else
for (i=0;i<n_vals;i++) val[i]*=units_factor;
} else if (units_bias != 0.0)
for (i=0;i<n_vals;i++) val[i]+=units_bias;
ret = super->pack_double(a,val,len);
switch (ret) {
case GRIB_CONSTANT_FIELD:
grib_buffer_replace(a, NULL, 0,1,1);
return GRIB_SUCCESS;
break;
case GRIB_SUCCESS:
case GRIB_SUCCESS:
break;
default:
default:
grib_context_log(a->parent->h->context,GRIB_LOG_ERROR,"unable to compute packing parameters\n");
return ret;
}
}
if((ret = grib_get_double_internal(a->parent->h,self->reference_value, &reference_value))
!= GRIB_SUCCESS)
return ret;
if((ret = grib_get_double_internal(a->parent->h,self->reference_value, &reference_value))
!= GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->binary_scale_factor, &binary_scale_factor))
!= GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->binary_scale_factor, &binary_scale_factor))
!= GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->bits_per_value,&bits_per_value)) !=
GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->bits_per_value,&bits_per_value)) !=
GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->decimal_scale_factor, &decimal_scale_factor))
!= GRIB_SUCCESS)
return ret;
if((ret = grib_get_long_internal(a->parent->h,self->decimal_scale_factor, &decimal_scale_factor))
!= GRIB_SUCCESS)
return ret;
decimal = grib_power(decimal_scale_factor,10) ;
divisor = grib_power(-binary_scale_factor,2);
decimal = grib_power(decimal_scale_factor,10) ;
divisor = grib_power(-binary_scale_factor,2);
simple_packing_size = (((bits_per_value*n_vals)+7)/8)*sizeof(unsigned char);
buf = grib_context_malloc_clear(a->parent->h->context,simple_packing_size+EXTRA_BUFFER_SIZE);
if(!buf) {
err = GRIB_OUT_OF_MEMORY;
goto cleanup;
}
simple_packing_size = (((bits_per_value*n_vals)+7)/8)*sizeof(unsigned char);
buf = grib_context_malloc_clear(a->parent->h->context,simple_packing_size+EXTRA_BUFFER_SIZE);
if(!buf) {
err = GRIB_OUT_OF_MEMORY;
goto cleanup;
}
if((err = grib_get_long_internal(a->parent->h,self->ni,&ni)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->nj,&nj)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->type_of_compression_used,&type_of_compression_used)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->target_compression_ratio,&target_compression_ratio)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->scanning_mode,&scanning_mode)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->list_defining_points,&list_defining_points)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->number_of_data_points,&number_of_data_points)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->ni,&ni)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->nj,&nj)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->type_of_compression_used,&type_of_compression_used)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->target_compression_ratio,&target_compression_ratio)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->scanning_mode,&scanning_mode)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->list_defining_points,&list_defining_points)) != GRIB_SUCCESS)
return err;
if((err = grib_get_long_internal(a->parent->h,self->number_of_data_points,&number_of_data_points)) != GRIB_SUCCESS)
return err;
width = ni;
height = nj;
width = ni;
height = nj;
if((scanning_mode & (1<<5)) != 0)
{
long tmp = width;
width = height;
height = tmp;
}
if((scanning_mode & (1<<5)) != 0)
{
long tmp = width;
width = height;
height = tmp;
}
/* The grid is not regular */
if(list_defining_points != 0)
{
width = *len;
height = 1;
}
/* The grid is not regular */
if(list_defining_points != 0)
{
width = *len;
height = 1;
}
/* There is a bitmap */
if(*len != number_of_data_points)
{
width = *len;
height = 1;
}
/* There is a bitmap */
if(*len != number_of_data_points)
{
width = *len;
height = 1;
}
if(width*height != *len)
{
/* fprintf(stderr,"width=%ld height=%ld len=%d\n",(long)width,(long)height,(long)*len); */
Assert(width*height == *len);
}
if(width*height != *len)
{
/* fprintf(stderr,"width=%ld height=%ld len=%d\n",(long)width,(long)height,(long)*len); */
Assert(width*height == *len);
}
switch( type_of_compression_used)
{
case 0:
Assert(target_compression_ratio == 255);
helper.compression = 0;
break;
switch( type_of_compression_used)
{
case 0:
Assert(target_compression_ratio == 255);
helper.compression = 0;
break;
case 1:
Assert(target_compression_ratio != 255);
Assert(target_compression_ratio != 0);
helper.compression = target_compression_ratio;
break;
case 1:
Assert(target_compression_ratio != 255);
Assert(target_compression_ratio != 0);
helper.compression = target_compression_ratio;
break;
default:
err = GRIB_NOT_IMPLEMENTED;
goto cleanup;
}
default:
err = GRIB_NOT_IMPLEMENTED;
goto cleanup;
}
helper.jpeg_buffer = buf;
helper.width = width;
helper.height = height;
helper.bits_per_value = bits_per_value;
helper.jpeg_buffer = buf;
helper.width = width;
helper.height = height;
helper.buffer_size = simple_packing_size + EXTRA_BUFFER_SIZE;
helper.values = val;
helper.no_values = n_vals;
helper.reference_value = reference_value;
helper.divisor = divisor;
helper.decimal = decimal;
helper.jpeg_length = 0;
/* See GRIB-438 */
if (bits_per_value == 0) {
const long bits_per_value_adjusted = 1;
grib_context_log(a->parent->h->context, GRIB_LOG_DEBUG,
"grib_accessor_class_data_jpeg2000_packing(%s) : bits per value was zero, changed to %d",
self->jpeg_lib==OPENJPEG_LIB ? "openjpeg" : "jasper", bits_per_value_adjusted);
bits_per_value = bits_per_value_adjusted;
}
helper.bits_per_value = bits_per_value;
helper.buffer_size = simple_packing_size + EXTRA_BUFFER_SIZE;
helper.values = val;
helper.no_values = n_vals;
helper.reference_value = reference_value;
helper.divisor = divisor;
helper.decimal = decimal;
helper.jpeg_length = 0;
switch (self->jpeg_lib) {
case OPENJPEG_LIB:
if ( (err = grib_openjpeg_encode(a->parent->h->context,&helper)) != GRIB_SUCCESS ) goto cleanup;
break;
case JASPER_LIB:
if ( (err = grib_jasper_encode(a->parent->h->context,&helper)) != GRIB_SUCCESS ) goto cleanup;
break;
}
case OPENJPEG_LIB:
if ( (err = grib_openjpeg_encode(a->parent->h->context,&helper)) != GRIB_SUCCESS ) goto cleanup;
break;
case JASPER_LIB:
if ( (err = grib_jasper_encode(a->parent->h->context,&helper)) != GRIB_SUCCESS ) goto cleanup;
break;
}
if(helper.jpeg_length > simple_packing_size)
grib_context_log(a->parent->h->context, GRIB_LOG_WARNING,
"grib_accessor_data_jpeg2000_packing(%s) : jpeg data (%ld) larger than input data (%ld)",
self->jpeg_lib==OPENJPEG_LIB ? "openjpeg" : "jasper",
helper.jpeg_length, simple_packing_size);
if(helper.jpeg_length > simple_packing_size)
grib_context_log(a->parent->h->context, GRIB_LOG_WARNING,
"grib_accessor_data_jpeg2000_packing(%s) : jpeg data (%ld) larger than input data (%ld)",
self->jpeg_lib==OPENJPEG_LIB ? "openjpeg" : "jasper",
helper.jpeg_length, simple_packing_size);
Assert( helper.jpeg_length <= helper.buffer_size);
Assert( helper.jpeg_length <= helper.buffer_size);
if(self->dump_jpg) {
FILE *f = fopen(self->dump_jpg,"w");
if(f) {
if(fwrite(helper.jpeg_buffer,helper.jpeg_length,1,f) != 1)
perror(self->dump_jpg);
if (fclose(f) != 0)
perror(self->dump_jpg);
}
else perror(self->dump_jpg);
}
if(self->dump_jpg) {
FILE *f = fopen(self->dump_jpg,"w");
if(f) {
if(fwrite(helper.jpeg_buffer,helper.jpeg_length,1,f) != 1)
perror(self->dump_jpg);
if (fclose(f) != 0)
perror(self->dump_jpg);
}
else perror(self->dump_jpg);
}
grib_buffer_replace(a, helper.jpeg_buffer, helper.jpeg_length, 1, 1);
grib_buffer_replace(a, helper.jpeg_buffer, helper.jpeg_length, 1, 1);
cleanup:
cleanup:
grib_context_free(a->parent->h->context,buf);
grib_context_free(a->parent->h->context,buf);
if(err == GRIB_SUCCESS)
err = grib_set_long_internal(a->parent->h,self->number_of_values, *len);
return err;
if(err == GRIB_SUCCESS)
err = grib_set_long_internal(a->parent->h,self->number_of_values, *len);
return err;
}
#else
static int unpack_double(grib_accessor* a, double* val, size_t *len)
{
grib_context_log(a->parent->h->context, GRIB_LOG_ERROR,
"jpeg support not enabled. Please rerun configure with the --with-jasper or --with-openjpeg option.");
return GRIB_NOT_IMPLEMENTED;
grib_context_log(a->parent->h->context, GRIB_LOG_ERROR,
"jpeg support not enabled. Please rerun configure with the --with-jasper or --with-openjpeg option.");
return GRIB_NOT_IMPLEMENTED;
}
static int pack_double(grib_accessor* a, const double* val, size_t *len)
{
grib_context_log(a->parent->h->context, GRIB_LOG_ERROR,
"jpeg support not enabled. Please rerun configure with the --with-jasper or --with-openjpeg option.");
return GRIB_NOT_IMPLEMENTED;
grib_context_log(a->parent->h->context, GRIB_LOG_ERROR,
"jpeg support not enabled. Please rerun configure with the --with-jasper or --with-openjpeg option.");
return GRIB_NOT_IMPLEMENTED;
}
#endif
static int unpack_double_element(grib_accessor* a, size_t idx, double* val)
{
size_t size;
double* values;
int err=0;
size_t size;
double* values;
int err=0;
err=grib_get_size(a->parent->h,"values",&size);
if (err) return err;
if (idx > size) return GRIB_INVALID_NEAREST;
values=grib_context_malloc_clear(a->parent->h->context,size*sizeof(double));
err=grib_get_double_array(a->parent->h,"values",values,&size);
if (err) return err;
*val=values[idx];
grib_context_free(a->parent->h->context,values);
return err;
err=grib_get_size(a->parent->h,"values",&size);
if (err) return err;
if (idx > size) return GRIB_INVALID_NEAREST;
values=grib_context_malloc_clear(a->parent->h->context,size*sizeof(double));
err=grib_get_double_array(a->parent->h,"values",values,&size);
if (err) return err;
*val=values[idx];
grib_context_free(a->parent->h->context,values);
return err;
}

286
src/grib_jasper_encoding.c
View File

@ -17,202 +17,198 @@
int grib_jasper_decode(grib_context *c,unsigned char *buf, size_t *buflen, double *values, size_t *no_values) {
/*jas_setdbglevel(99999);*/
/*jas_setdbglevel(99999);*/
jas_image_t *image = NULL;
jas_stream_t *jpeg = NULL;
int code = GRIB_SUCCESS;
jas_matrix_t *matrix = NULL;
jas_image_cmpt_t *p;
int i,j,k;
jas_image_t *image = NULL;
jas_stream_t *jpeg = NULL;
int code = GRIB_SUCCESS;
jas_matrix_t *matrix = NULL;
jas_image_cmpt_t *p;
int i,j,k;
jpeg = jas_stream_memopen((void*)buf,*buflen);
if(!jpeg) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
jpeg = jas_stream_memopen((void*)buf,*buflen);
if(!jpeg) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
image = jpc_decode(jpeg,NULL);
if(!image) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
image = jpc_decode(jpeg,NULL);
if(!image) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
p = image->cmpts_[0];
p = image->cmpts_[0];
if (image->numcmpts_ != 1 ) {
/* Image not gray scale */
code = GRIB_DECODING_ERROR;
goto cleanup;
}
if (image->numcmpts_ != 1 ) {
/* Image not gray scale */
code = GRIB_DECODING_ERROR;
goto cleanup;
}
matrix = jas_matrix_create(jas_image_height(image), jas_image_width(image));
matrix = jas_matrix_create(jas_image_height(image), jas_image_width(image));
if(!matrix) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
if(!matrix) {
code = GRIB_DECODING_ERROR;
goto cleanup;
}
jas_image_readcmpt(image,0,0,0,jas_image_width(image), jas_image_height(image),matrix);
jas_image_readcmpt(image,0,0,0,jas_image_width(image), jas_image_height(image),matrix);
Assert(p->height_ * p->width_ == *no_values);
Assert(p->height_ * p->width_ == *no_values);
k=0;
for (i=0;i<p->height_;i++)
for (j=0;j<p->width_;j++) {
values[k++] = matrix->rows_[i][j];
}
k=0;
for (i=0;i<p->height_;i++)
for (j=0;j<p->width_;j++) {
values[k++] = matrix->rows_[i][j];
}
cleanup:
if(matrix) jas_matrix_destroy(matrix);
if(image) jas_image_destroy(image);
if(jpeg) jas_stream_close(jpeg);
if(matrix) jas_matrix_destroy(matrix);
if(image) jas_image_destroy(image);
if(jpeg) jas_stream_close(jpeg);
return code;
return code;
}
int grib_jasper_encode(grib_context *c, j2k_encode_helper *helper) {
int code = GRIB_SUCCESS;
int jaserr;
int code = GRIB_SUCCESS;
int jaserr;
char opts[MAXOPTSSIZE];
double reference_value = helper->reference_value;
double decimal = helper->decimal;
double divisor = helper->divisor;
const double* values = helper->values;
long no_values = helper->no_values;
long bits8;
int i;
char opts[MAXOPTSSIZE];
double reference_value = helper->reference_value;
double decimal = helper->decimal;
double divisor = helper->divisor;
const double* values = helper->values;
long no_values = helper->no_values;
long bits8;
int i;
size_t buflen = 0;
unsigned char *encoded = NULL;
unsigned char *p = NULL;
size_t buflen = 0;
unsigned char *encoded = NULL;
unsigned char *p = NULL;
jas_image_t image = {0,};
jas_stream_t *jpcstream = 0;
jas_stream_t *istream = 0;
jas_image_cmpt_t cmpt = {0,};
jas_image_cmpt_t *pcmpt = 0;
jas_image_t image = {0,};
jas_stream_t *jpcstream = 0;
jas_stream_t *istream = 0;
jas_image_cmpt_t cmpt = {0,};
jas_image_cmpt_t *pcmpt = 0;
image.tlx_ = 0;
image.tly_ = 0;
image.brx_ = helper->width;
image.bry_ = helper->height;
image.numcmpts_ = 1;
image.maxcmpts_ = 1;
image.clrspc_ = JAS_CLRSPC_SGRAY;
image.cmprof_ = 0;
image.inmem_ = 1;
image.tlx_ = 0;
image.tly_ = 0;
image.brx_ = helper->width;
image.bry_ = helper->height;
image.numcmpts_ = 1;
image.maxcmpts_ = 1;
image.clrspc_ = JAS_CLRSPC_SGRAY;
image.cmprof_ = 0;
image.inmem_ = 1;
cmpt.tlx_ = 0;
cmpt.tly_ = 0;
cmpt.hstep_ = 1;
cmpt.vstep_ = 1;
cmpt.width_ = helper->width;
cmpt.height_ = helper->height;
cmpt.type_ = JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_GRAY_Y);
cmpt.prec_ = helper->bits_per_value;
cmpt.sgnd_ = 0;
cmpt.cps_ = (helper->bits_per_value+7)/8;
cmpt.tlx_ = 0;
cmpt.tly_ = 0;
cmpt.hstep_ = 1;
cmpt.vstep_ = 1;
cmpt.width_ = helper->width;
cmpt.height_ = helper->height;
cmpt.type_ = JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_GRAY_Y);
cmpt.prec_ = helper->bits_per_value;
cmpt.sgnd_ = 0;
cmpt.cps_ = (helper->bits_per_value+7)/8;
/* Simple packing encoding */
bits8 = (helper->bits_per_value+7)/8*8;
Assert(bits8 > 0);
encoded = grib_context_malloc_clear(c,bits8/8*no_values);
/* Simple packing encoding */
if(!encoded) {
code = GRIB_OUT_OF_MEMORY;
goto cleanup;
}
buflen = 0;
p = encoded;
bits8 = (helper->bits_per_value+7)/8*8;
encoded = grib_context_malloc_clear(c,bits8/8*no_values);
for(i=0;i< no_values;i++){
long blen = bits8;
unsigned long unsigned_val = (unsigned long)((((values[i]*decimal)-(reference_value))*divisor)+0.5);
while(blen >= 8)
{
blen -= 8;
*p = (unsigned_val >> blen);
p++;
buflen++;
}
}
if(!encoded) {
code = GRIB_OUT_OF_MEMORY;
goto cleanup;
}
/*jas_init();*/
buflen = 0;
p = encoded;
opts[0] = 0;
for(i=0;i< no_values;i++){
long blen = bits8;
unsigned long unsigned_val = (unsigned long)((((values[i]*decimal)-(reference_value))*divisor)+0.5);
while(blen >= 8)
{
blen -= 8;
*p = (unsigned_val >> blen);
p++;
buflen++;
}
}
/*jas_init();*/
opts[0] = 0;
if( helper->compression != 0)
{
/* Lossy */
if( helper->compression != 0)
{
/* Lossy */
#ifndef GRIB_ON_WINDOWS
snprintf (opts, MAXOPTSSIZE, "mode=real\nrate=%f", 1.0/helper->compression);
snprintf (opts, MAXOPTSSIZE, "mode=real\nrate=%f", 1.0/helper->compression);
#else
/* Microsoft Windows Visual Studio support */
_snprintf(opts, MAXOPTSSIZE, "mode=real\nrate=%f", 1.0/helper->compression);
/* Microsoft Windows Visual Studio support */
_snprintf(opts, MAXOPTSSIZE, "mode=real\nrate=%f", 1.0/helper->compression);
#endif
}
}
Assert(cmpt.width_ * cmpt.height_ * cmpt.cps_ == buflen);
Assert(cmpt.width_ * cmpt.height_ * cmpt.cps_ == buflen);
pcmpt = &cmpt;
image.cmpts_ = &pcmpt;
pcmpt = &cmpt;
image.cmpts_ = &pcmpt;
istream = jas_stream_memopen((char *)encoded,buflen);
cmpt.stream_ = istream;
istream = jas_stream_memopen((char *)encoded,buflen);
cmpt.stream_ = istream;
jpcstream = jas_stream_memopen((char*)helper->jpeg_buffer,helper->buffer_size);
jaserr = jpc_encode(&image,jpcstream,opts);
jpcstream = jas_stream_memopen((char*)helper->jpeg_buffer,helper->buffer_size);
jaserr = jpc_encode(&image,jpcstream,opts);
if(jaserr != 0) {
/* increase the number of guard bits */
strcat(opts,"\nnumgbits=4");
grib_context_log(c, GRIB_LOG_ERROR, "JASPER: error %d, increasing the number of guard bits", jaserr);
jaserr=jas_stream_close(istream);
istream = jas_stream_memopen((char *)encoded,buflen);
jaserr=jas_stream_close(jpcstream);
jpcstream = jas_stream_memopen((char*)helper->jpeg_buffer,helper->buffer_size);
jaserr = jpc_encode(&image,jpcstream,opts);
}
if(jaserr != 0) {
/* increase the number of guard bits */
strcat(opts,"\nnumgbits=4");
grib_context_log(c, GRIB_LOG_ERROR, "JASPER: error %d, increasing the number of guard bits", jaserr);
jaserr=jas_stream_close(istream);
istream = jas_stream_memopen((char *)encoded,buflen);
jaserr=jas_stream_close(jpcstream);
jpcstream = jas_stream_memopen((char*)helper->jpeg_buffer,helper->buffer_size);
jaserr = jpc_encode(&image,jpcstream,opts);
}
if ( jaserr != 0 ) {
grib_context_log(c, GRIB_LOG_ERROR, "JASPER: error %d", jaserr);
code = GRIB_ENCODING_ERROR;
goto cleanup;
}
if ( jaserr != 0 ) {
grib_context_log(c, GRIB_LOG_ERROR, "JASPER: error %d", jaserr);
code = GRIB_ENCODING_ERROR;
goto cleanup;
}
helper->jpeg_length = jpcstream->rwcnt_;
jaserr=jas_stream_close(istream); istream = 0;
jaserr=jas_stream_close(jpcstream);jpcstream = 0;
helper->jpeg_length = jpcstream->rwcnt_;
jaserr=jas_stream_close(istream); istream = 0;
jaserr=jas_stream_close(jpcstream);jpcstream = 0;
cleanup:
grib_context_free(c,encoded);
if(istream) jas_stream_close(istream);
if(jpcstream) jas_stream_close(jpcstream);
grib_context_free(c,encoded);
if(istream) jas_stream_close(istream);
if(jpcstream) jas_stream_close(jpcstream);
return code;
return code;
}
#else
int grib_jasper_decode(grib_context *c, unsigned char *buf, size_t *buflen, double *val, size_t *n_vals) {
grib_context_log(c, GRIB_LOG_ERROR,
"grib_accessor_data_jpeg2000_packing : jasper jpeg support not enabled. Please rerun configure with --with-jasper-support");
return GRIB_NOT_IMPLEMENTED;
grib_context_log(c, GRIB_LOG_ERROR,
"grib_accessor_data_jpeg2000_packing : jasper jpeg support not enabled. Please rerun configure with --with-jasper-support");
return GRIB_NOT_IMPLEMENTED;
}
int grib_jasper_encode(grib_context *c, j2k_encode_helper *helper) {
grib_context_log(c, GRIB_LOG_ERROR,
"grib_accessor_data_jpeg2000_packing : jasper jpeg support not enabled. Please rerun configure with --with-jasper-support");
return GRIB_NOT_IMPLEMENTED;
grib_context_log(c, GRIB_LOG_ERROR,
"grib_accessor_data_jpeg2000_packing : jasper jpeg support not enabled. Please rerun configure with --with-jasper-support");
return GRIB_NOT_IMPLEMENTED;
}
#endif