Странное считывание из бинарного файла

#include <fcntl.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <stdio.h>
#include <math.h>
 
#define LCH 245
#define RCH 255
#define sqr(x) pow(x, 2)
 
struct line {
    double k;
    double b;
};
 
struct line background() /* Фон */ 
{
    int i, c, in;
    long int lb[3], rb[3];
    float avlb, avrb, s;
    struct line co;
    
    in = open("EU152_1K.SPK", O_RDONLY);
    
    for (i = -2; i <= 0; i++) { 
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &c, 4);
        lb[2+i] = c;
    }
    
    for (i = 0; i <= 2; i++) {
        lseek(in, 512 + 4*(i+RCH), SEEK_SET);
        read(in, &c, 4);
        rb[i] = c;
    }
    
    avlb = (lb[0]+lb[1]+lb[2])/3;
    avrb = (rb[0]+rb[1]+rb[2])/3;
    
    s = (avrb+avlb)*(RCH-LCH+1)/2;
    co.k = (avrb-avlb)/(RCH-LCH+1);
    co.b = lb[1] - co.k*(LCH-1);
    close(in);
    
    return co;
}
 
void area() /* Вычисляет площадь и ошибку S */
{
    int i, in;
    long int N, c;
    struct line co;
    float s_err = 0, s = 0;
    
    co = background();
    in = open("EU152_1K.SPK", O_RDONLY);
    for (i = 0; i <= (RCH - LCH); i++) {
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &N, 4);
        c = N;
        s = s + c - ((LCH+i)*co.k + co.b);
        s_err = s_err + c + co.k*(LCH+i)+co.b;
    }
    close(in);
    printf("S = %.2f; S_ERR = %.2f\n", s, sqrt(s_err));
}
 
void p() /* Считаем центр тяжести и ошибку p */
{
    int i, in;
    long int c;
    struct line co;
    float p = 0, s = 0, p_err = 0;
    
    co = background();
    in = open("EU152_1K.SPK", O_RDONLY);
    
    for (i = 0; i <= (RCH - LCH); i++) {
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &c, 4);
        s = s + c - ((LCH+i)*co.k + co.b);
        p = p + (c - ((LCH+i)*co.k + co.b))*(LCH + i);
    }
    
    for (i = 0; i <= (RCH - LCH); i++) {
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &c, 4);
        p_err = p_err + (c + co.k*(LCH+i)+co.b)/sqr(c - co.k*(LCH+i)-co.b) + sqr(227.51)/sqr(38533);
    }
    close(in);
    
    printf("P = %.2f; P_ERR = %.4f\n", p/s, sqrt(p_err));
}
 
void delta() /* Считаем delta */
{
    int i, in;
    long int c;
    struct line co;
    float d = 0, s = 0, p = 0;
    
    co = background();
    in = open("EU152_1K.SPK", O_RDONLY);
    
    for (i = 0; i <= (RCH - LCH); i++) {
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &c, 4);
        s = s + c - ((LCH+i)*co.k + co.b);
        p = p + (c - ((LCH+i)*co.k + co.b))*(LCH + i);
    }
    
    for (i = 0; i <= (RCH - LCH); i++) {
        lseek(in, 512 + 4*(i+LCH), SEEK_SET);
        read(in, &c, 4);
        printf("c = %ld\n", c);
        d = d + (c - ((LCH+i)*co.k + co.b))*sqr(LCH+i - p);
    }
    
    printf("d = %.2f; \n", sqrt(d)/s);
}
 
main()
{
    p();
    area();
    delta();
}

lseek(in, 512 + 4*(i+LCH), SEEK_SET);

read(in, &c, 4);

printf("%d\n", c);

long int c;
read(in, &c, 4);

read(in, &c, sizeof c);

long int c;
...
read(in, &c, 4);

printf("%ld", (long) c)