Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1213 1471 (N24W09 L=218)   307     61        0 (1170/1450-1249/1486) -11.3 W06-W11 N22-N25
      2  297 1343 (N12E53 L=156)   386     54        2 ( 277/1316- 315/1366)  56.3 E55-E51 N10-N13
      3 1657 1147 (N07W40 L=249)   898    222        1 (1582/1066-1737/1241) -31.6 W34-W47 N03-N12
      4 1756 1213 (N13W49 L=257)   362     77        0 (1721/1188-1780/1232) -21.8 W45-W51 N11-N13
      5  850 1079 (S02E11 L=198)   512    120        1 ( 829/1009- 870/1141) -74.1 E12-E09 S06-N01
      6  834  977 (S08E11 L=197)   201     45        0 ( 804/ 959- 864/ 994)  33.7 E12-E09 S09-S07
      7 1396  914 (S09W23 L=232)   194     74        0 (1368/ 880-1429/ 947)  53.1 W22-W25 S10-S06
      8 1867  851 (S08W63 L=272)   756    109        3 (1833/ 807-1895/ 896)  56.3 W59-W67 S10-S04
      9 1284  860 (S13W17 L=225)   243    112        0 (1228/ 839-1328/ 877)  -9.5 W13-W19 S13-S11
     10 1572  728 (S18W38 L=247)  1992    359        1 (1408/ 619-1731/ 790)  26.6 W28-W51 S26-S13
     11  345  654 (S30E48 L=160)   605    157        1 ( 310/ 607- 380/ 712)  40.6 E52-E43 S32-S26
     12 1243  492 (S36W20 L=229)   510    187        0 (1148/ 449-1338/ 539)  30.3 W13-W25 S40-S32
  

    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