Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1478 1768 (N56W22 L=193)   821    142        1 (1432/1744-1529/1793)  24.8 W16-W29 N51-N60
      2  714 1483 (N15E28 L=143)   576    203        2 ( 692/1388- 738/1570) -67.4 E31-E26 N08-N20
      3 1404 1302 (N17W18 L=188)   522    186        0 (1382/1208-1447/1380)  80.5 W16-W19 N10-N22
      4  960 1291 (N07E09 L=162)   386    120        0 ( 921/1235- 998/1355)  82.9 E09-E07 N03-N11
      5  341 1328 (N00E51 L=120)   640    111        1 ( 302/1301- 383/1353)   0.0 E54-E46 S01-N02
      6  221 1198 (S08E60 L=111)   720    151        1 ( 174/1176- 262/1221)   6.3 E64-E54 S09-S06
      7 1344 1197 (N09W15 L=186)   285    101        1 (1311/1166-1394/1215)  11.3 W13-W18 N07-N10
      8  876  797 (S23E05 L=166)   557    171        1 ( 780/ 764- 951/ 843)  35.0 E10-W00 S25-S18
      9 1466  808 (S10W31 L=202)   267     64        1 (1442/ 774-1486/ 841)  90.0 W29-W32 S11-S07
     10 1764  737 (S07W56 L=227)   384    121        0 (1723/ 706-1808/ 755)  35.5 W53-W60 S09-S04
     11 1302  674 (S21W24 L=195)   230     54        1 (1273/ 650-1335/ 695)  53.1 W22-W26 S23-S19
     12 1489  402 (S32W49 L=220)  3975    887        6 (1179/ 351-1812/ 522)  31.9 W28-W74 S43-S15
  

    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