Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1102 1856 (N51E23 L= 82)   833     89        1 (1046/1835-1150/1873)  45.0 E26-E18 N47-N53
      2  675 1779 (N32E48 L= 57)   442     60        0 ( 651/1749- 698/1803) -20.6 E52-E44 N29-N34
      3 1443 1689 (N47W13 L=118)   482     69        1 (1417/1651-1472/1714) -68.2 W10-W15 N43-N49
      4  651 1648 (N22E41 L= 64)   505     99        0 ( 611/1611- 700/1671)  45.0 E42-E37 N18-N24
      5  993 1541 (N22E14 L= 91)   627    271        0 ( 923/1486-1078/1606)  -5.2 E20-E09 N18-N25
      6  510 1466 (N08E43 L= 62)   339     53        0 ( 487/1438- 538/1500)  90.0 E44-E41 N06-N10
      7 1496  434 (S28W51 L=156)  1117    249        1 (1424/ 349-1604/ 497)  71.6 W48-W56 S33-S21
      8  932  240 (S60W25 L=131)   790    166        0 ( 881/ 188-1001/ 290)  -8.5 W16-W37 S62-S56
      9 1126  158 (S58W56 L=161)   976     94        1 (1102/ 118-1161/ 198) -30.3 W46-W67 S61-S54
  

    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