Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1885 (N60E38 L=118)  1086    115        1 ( 932/1851-1012/1926) -13.5 E56-E30 N56-N64
      2  645 1580 (N27E39 L=118)   461    153        0 ( 601/1557- 708/1625)  -5.2 E44-E33 N25-N29
      3  319 1552 (N20E70 L= 86)  1014     99        2 ( 279/1525- 359/1570)   0.0 E79-E63 N18-N20
      4  544 1469 (N18E41 L=115)   476     71        1 ( 524/1438- 570/1502) -38.7 E43-E38 N16-N20
      5 1555 1254 (N23W31 L=187)   349    126        0 (1539/1197-1569/1314) -90.0 W29-W32 N19-N26
      6 1878 1282 (N31W70 L=226)  1495     94        3 (1842/1243-1906/1313)  -8.7 W63-W78 N27-N33
      7 1816  932 (N08W58 L=215)   586    140        0 (1787/ 871-1841/ 988)  69.4 W56-W60 N03-N11
      8  231  854 (S24E58 L= 98)   852    148        0 ( 210/ 780- 268/ 911)  52.1 E60-E52 S29-S20
      9  360  810 (S25E43 L=113)   388    146        1 ( 305/ 788- 411/ 868) -11.3 E48-E38 S26-S22
     10 1035  614 (S25W08 L=164)  2574    576        7 ( 806/ 512-1247/ 750)  43.7 E05-W18 S35-S13
     11 1536  623 (S16W42 L=198)   618    145        1 (1490/ 577-1597/ 662)  54.5 W39-W45 S19-S12
  

    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