Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  978 1605 (N32E18 L=295)   351    125        1 ( 928/1572-1024/1645)  69.4 E20-E16 N28-N36
      2 1440 1364 (N27W19 L=333)   727    141        4 (1370/1325-1494/1400)   0.0 W13-W23 N25-N29
      3  899 1289 (N10E13 L=300)   600    254        1 ( 787/1215- 977/1357)  -8.1 E21-E07 N07-N13
      4 1682  864 (N03W46 L=359)   593    130        3 (1640/ 832-1726/ 894)  39.8 W43-W49 N00-N05
      5  669  814 (S22E17 L=296)   340    103        0 ( 614/ 798- 712/ 832)  23.2 E21-E14 S23-S20
      6  187  754 (S37E67 L=247)   660     72        0 ( 163/ 734- 217/ 785) -45.0 E71-E61 S37-S34
      7 1603  463 (S22W58 L= 12)   854    218        2 (1510/ 413-1691/ 509)  39.3 W52-W64 S26-S17
      8 1353  452 (S29W38 L=352)   393    154        0 (1299/ 405-1409/ 500)  69.4 W36-W39 S32-S24
  

    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