Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  916 1829 (N40E38 L= 26)  1409    247        3 ( 841/1752-1002/1910) -15.9 E54-E26 N35-N46
      2  620 1809 (N31E60 L=  5)  1289    133        0 ( 567/1757- 652/1859) -48.0 E68-E56 N25-N35
      3  811 1668 (N25E33 L= 31)   630    179        0 ( 779/1603- 849/1760) -32.7 E41-E27 N21-N30
      4 1279 1513 (N28W01 L= 65)  1035    343        0 (1168/1369-1350/1624) -30.5 E08-W09 N21-N33
      5  890 1566 (N20E24 L= 41)   336     80        1 ( 856/1537- 927/1590)  68.2 E25-E21 N17-N22
      6  149 1074 (S24E69 L=356)  5758    434       11 ( 103/ 938- 224/1219) -86.6 E79-E57 S31-S14
      7  203 1083 (S22E61 L=  4)   722    122        1 ( 177/1028- 232/1124)  71.6 E63-E57 S25-S19
      8 1100  606 (S28W17 L= 82)   204     71        0 (1079/ 570-1114/ 643) -31.0 W14-W19 S29-S26
      9 1565  362 (S26W65 L=129)  1166    144        0 (1509/ 314-1628/ 421)  90.0 W62-W67 S29-S20
     10 1309  315 (S38W49 L=113)   536    123        1 (1243/ 297-1365/ 338)   9.5 W42-W54 S39-S35
  

    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