Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  960 1752 (N48W04 L= 13)   294     45        0 ( 945/1723- 975/1778)  68.2 W02-W05 N45-N50
      2 1453 1545 (N27W36 L= 46)  1462    304        3 (1291/1513-1570/1584) -10.8 W24-W45 N24-N31
      3  627 1510 (N33E24 L=346)   357     97        0 ( 583/1494- 679/1525) -14.0 E27-E19 N31-N34
      4  865 1411 (N24E07 L=  3)  1054    272        0 ( 761/1345- 994/1458) -29.7 E13-W01 N18-N27
      5 1357 1441 (N21W26 L= 36)   412     89        0 (1310/1426-1408/1456)  -8.1 W22-W29 N19-N22
      6  490 1435 (N29E34 L=335)   449     99        0 ( 441/1414- 540/1454) -18.4 E38-E29 N27-N30
      7 1790 1403 (N15W62 L= 72)   700     84        2 (1758/1377-1828/1429)   0.0 W57-W67 N13-N17
      8 1418 1211 (N06W26 L= 36)   469    127        1 (1375/1166-1465/1263) -54.5 W23-W28 N03-N10
      9 1656  967 (S10W41 L= 51)   431    104        0 (1635/ 914-1686/1017)  59.0 W39-W43 S12-S06
     10 1874  883 (S17W64 L= 73)  4571    271       10 (1810/ 779-1946/1001)  47.3 W55-W75 S22-S09
     11 1787  804 (S21W54 L= 64)  1964    457        6 (1723/ 588-1851/ 981)  76.0 W51-W59 S34-S10
     12  300  632 (S18E58 L=311)   601     78        0 ( 271/ 602- 321/ 675)  33.7 E62-E55 S19-S15
  

    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