Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1308 1802 (N53W17 L=262)   665     95        1 (1275/1764-1340/1840)  35.0 W13-W23 N48-N58
      2  587 1679 (N34E41 L=203)   928    127        1 ( 558/1602- 612/1732) -63.4 E45-E37 N28-N38
      3 1524 1480 (N28W31 L=275)   404    112        0 (1473/1450-1567/1497)  20.6 W26-W34 N25-N30
      4  279 1404 (N14E58 L=186)  1079    196        3 ( 243/1315- 322/1483) -36.3 E66-E51 N08-N19
      5 1757 1176 (N12W49 L=293)   500     58        0 (1724/1144-1782/1202) -21.8 W45-W51 N10-N13
      6  136  796 (S23E74 L=171)  1327    119        2 ( 112/ 742- 164/ 860) -69.4 E79-E69 S26-S18
      7 1379  637 (S25W28 L=272)   228     47        1 (1361/ 613-1399/ 659) -45.0 W25-W29 S26-S23
      8 1553  359 (S41W58 L=302)  1064     78        1 (1514/ 317-1578/ 409) -26.6 W52-W65 S43-S36
      9  626  317 (S56E34 L=211)   358     47        1 ( 606/ 293- 641/ 349)  45.0 E37-E30 S58-S53
     10  561  300 (S58E45 L=199)   845     70        1 ( 535/ 266- 589/ 339)  25.0 E55-E38 S61-S54
  

    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