Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  619 1713 (N48E42 L=344)   501     88        0 ( 598/1683- 639/1741) -35.5 E44-E37 N45-N50
      2  564 1637 (N41E42 L=344)   337     36        0 ( 544/1618- 582/1657)  71.6 E43-E40 N39-N43
      3  882 1489 (N32E12 L= 14)   480     96        0 ( 836/1457- 937/1523) -26.6 E15-E07 N30-N34
      4  462 1383 (N23E42 L=344)   643    125        1 ( 439/1325- 480/1450) -82.9 E43-E39 N19-N27
      5 1146 1414 (N28W07 L= 33)   219     45        0 (1124/1387-1164/1443) -45.0 W05-W08 N26-N29
      6 1410 1321 (N22W25 L= 51)  1722    376        2 (1239/1228-1575/1423) -27.6 W13-W36 N16-N28
      7  942 1307 (N20E06 L= 20)   863    180        2 ( 896/1219- 988/1383) -59.0 E09-E02 N14-N25
      8  740 1085 (N06E18 L=  8)  1339    324        2 ( 654/ 978- 869/1217) -47.3 E24-E09 N00-N13
      9  604  940 (S04E27 L=360)   231     64        1 ( 566/ 923- 634/ 954)   0.0 E29-E24 S04-S02
     10 1322  256 (S51W34 L= 60)   430    152        0 (1283/ 230-1357/ 282) -26.6 W27-W39 S54-S48
  

    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