Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  734 1732 (N28E43 L=189)   416     84        0 ( 715/1693- 756/1770) -63.4 E45-E40 N24-N30
      2  424 1646 (N15E62 L=170)   436     87        0 ( 409/1603- 444/1688) -45.0 E64-E58 N12-N18
      3  405 1386 (S01E49 L=184)   342     93        1 ( 366/1347- 435/1415) -14.0 E52-E44 S02-N00
      4  244 1123 (S19E55 L=177)  1164    167        2 ( 182/1076- 319/1160)  26.6 E62-E48 S23-S16
      5 1241  969 (S03W13 L=245)   306    119        1 (1227/ 910-1254/1024) -71.6 W11-W15 S06-N00
      6  614  849 (S27E20 L=212)   409    133        0 ( 591/ 787- 648/ 891)  81.9 E20-E18 S30-S23
      7  550  524 (S49E19 L=213)   254     50        0 ( 532/ 495- 569/ 552) -36.9 E20-E16 S51-S47
      8  700  467 (S48E04 L=228)   756    153        1 ( 622/ 436- 772/ 494)  13.0 E11-W02 S51-S45
      9 1378  426 (S28W44 L=276)   363     68        1 (1342/ 411-1416/ 443)  18.4 W40-W46 S29-S26
     10 1555  344 (S26W65 L=297)   424     26        0 (1535/ 323-1575/ 362) -71.6 W62-W67 S27-S24
     11 1174  217 (S48W50 L=282)   813    124        1 (1130/ 172-1225/ 255)  -6.0 W42-W61 S48-S45
  

    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