Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  905 1782 (N60W12 L=287)   396     35        1 ( 883/1756- 925/1806)  33.7 W09-W16 N57-N62
      2 1454 1610 (N33W47 L=321)  1588    216        3 (1367/1570-1555/1665) -45.0 W40-W53 N28-N38
      3  592 1463 (N41E24 L=250)   988    212        1 ( 487/1438- 692/1485) -11.9 E34-E15 N38-N44
      4  249 1457 (N45E69 L=206)  1069    121        2 ( 206/1423- 314/1476)   0.0 E83-E57 N42-N46
      5 1755 1334 (N09W58 L=332)   349     71        0 (1738/1297-1769/1368)  26.6 W54-W60 N07-N10
      6 1745  668 (S31W47 L=321)  2799    364        4 (1628/ 570-1857/ 834)  23.5 W34-W59 S36-S22
      7  236  618 (S08E69 L=206)  2140    133        6 ( 198/ 568- 281/ 670) -40.6 E76-E63 S11-S04
      8 1478  644 (S25W22 L=297)   394     75        0 (1435/ 628-1515/ 661) -18.4 W19-W25 S26-S23
      9  633  461 (S22E38 L=237)   272     31        0 ( 614/ 439- 652/ 483) -90.0 E38-E36 S23-S20
     10 1047  270 (S44E16 L=258)  3523    499        5 ( 853/ 181-1272/ 384)  10.5 E39-W04 S48-S38
  

    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