Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  325 1357 (N28E54 L=146)   410     53        0 ( 294/1343- 357/1372)   0.0 E57-E50 N27-N29
      2 1744 1307 (N16W56 L=255)   460    128        0 (1704/1247-1776/1345) -60.3 W52-W57 N11-N18
      3  227 1089 (N11E59 L=140)   861    164        1 ( 210/1021- 252/1164) -83.7 E61-E56 N07-N16
      4  970  943 (N00E04 L=196)   263     63        2 ( 935/ 918-1005/ 963)   0.0 E05-E01 S01-N01
      5 1465  920 (S05W27 L=227)   672    175        1 (1390/ 898-1547/ 951) -20.0 W22-W33 S06-S02
      6 1729  869 (S11W48 L=248)   743    124        0 (1686/ 837-1785/ 914) -31.0 W44-W54 S13-S07
      7  318  707 (S12E54 L=146)   970    231        2 ( 279/ 611- 364/ 789)  70.0 E59-E48 S18-S07
      8 1509  587 (S27W32 L=231)   610    198        1 (1481/ 504-1547/ 679)  90.0 W29-W33 S33-S21
      9  867  639 (S18E13 L=187)   301     75        0 ( 841/ 615- 913/ 655)  11.3 E14-E09 S19-S17
     10  769  512 (S27E21 L=178)   515     82        2 ( 734/ 486- 813/ 540) -21.8 E23-E18 S28-S24
     11  712  405 (S35E29 L=170)   364     54        1 ( 693/ 378- 739/ 430)  76.0 E30-E26 S36-S32
     12  683  249 (S49E43 L=156)   469     66        0 ( 664/ 225- 709/ 280) -76.0 E46-E39 S51-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