Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1875 (N72W06 L=152)   695     65        2 ( 913/1846- 958/1900) -60.3 W02-W11 N67-N76
      2 1594 1565 (N33W55 L=201)   512     95        0 (1551/1546-1634/1582)  -6.3 W49-W58 N31-N33
      3 1135 1532 (N37W14 L=160)   366     99        0 (1087/1504-1181/1555)  14.0 W09-W17 N35-N38
      4 1215 1521 (N35W20 L=166)   280     77        0 (1197/1484-1238/1551)  53.1 W18-W22 N32-N36
      5 1560 1496 (N29W47 L=193)   624    113        1 (1530/1436-1588/1543)  81.9 W44-W48 N24-N31
      6 1718 1474 (N24W62 L=208)   899    121        1 (1672/1426-1766/1524) -90.0 W55-W67 N20-N27
      7  820 1347 (N27E11 L=135)   529    133        1 ( 786/1291- 865/1406) -58.0 E13-E08 N23-N31
      8  405 1285 (N26E43 L=103)  2612    464        5 ( 297/1174- 577/1347)  11.8 E54-E28 N18-N29
      9 1280 1268 (N18W19 L=165)   384    121        2 (1236/1230-1326/1308) -38.7 W16-W21 N15-N20
     10  414 1075 (N12E41 L=105)   750    153        3 ( 357/1027- 464/1123) -36.9 E46-E36 N09-N15
     11 1751  859 (S13W50 L=196)   583      0       -1 (1723/ 830-1775/ 883)  33.7 W47-W52 S15-S11
     12  693  797 (S06E23 L=123)   270     93        0 ( 647/ 770- 735/ 813) -21.8 E25-E20 S08-S04
     13  257  745 (S08E60 L= 86)   657     59        1 ( 232/ 716- 288/ 771)  45.0 E63-E56 S09-S06
     14 1067  478 (S31E03 L=144)   255    103        1 (1018/ 467-1116/ 491) -18.4 E05-W01 S31-S29
     15  451  403 (S32E56 L= 90)   645    114        0 ( 428/ 355- 483/ 461)  66.0 E58-E53 S36-S27
     16  875  404 (S34E18 L=128)   444    125        1 ( 813/ 389- 943/ 418)  -5.7 E22-E12 S35-S32
     17 1256  250 (S54W13 L=159)   413     37        0 (1233/ 235-1280/ 271) -31.0 W10-W15 S56-S51
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y