Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  675 1826 (N40E58 L= 60)  1919    125        2 ( 639/1770- 725/1884) -14.5 E78-E47 N36-N44
      2  905 1657 (N33E25 L= 92)   459     97        0 ( 882/1615- 939/1701) -81.9 E26-E23 N29-N36
      3  340 1537 (N13E63 L= 54)  2686    249        6 ( 256/1461- 402/1643)  81.9 E69-E54 N06-N20
      4  441 1599 (N19E57 L= 60)   837     90        1 ( 406/1567- 475/1635)   0.0 E61-E52 N16-N21
      5  256 1335 (S01E61 L= 56)   531    135        0 ( 225/1271- 272/1381)  90.0 E64-E59 S04-N02
      6 1024 1063 (S01E01 L=116)  1384    350        1 ( 976/ 881-1091/1206) -51.3 E06-W06 S09-N07
      7 1809 1076 (N19W53 L=170)  1821    276        0 (1726/ 964-1873/1182) -25.3 W43-W62 N13-N23
      8  507 1106 (S09E33 L= 84)   478    169        0 ( 433/1086- 574/1155)  33.7 E38-E29 S11-S04
      9 1144  705 (S18W15 L=132)   228     51        0 (1127/ 675-1160/ 734) -56.3 W13-W16 S19-S16
     10 1588  596 (S14W48 L=165)  1706    281        2 (1470/ 531-1714/ 648)  40.9 W41-W56 S20-S07
     11 1593  401 (S25W62 L=179)  2512    183        4 (1504/ 370-1684/ 433)  10.8 W51-W72 S27-S20
     12  944  370 (S45W14 L=131)   650     85        1 ( 896/ 342- 982/ 411)   7.1 W10-W18 S47-S42
  

    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