Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  613 1672 (N22E45 L=238)  1497    186        5 ( 531/1620- 705/1705)  24.4 E51-E39 N17-N25
      2 1519 1513 (N35W19 L=303)  1123    204        1 (1413/1462-1606/1562)   3.4 W10-W27 N33-N37
      3  284 1502 (N06E65 L=218)   785    104        5 ( 239/1457- 313/1541)  51.3 E70-E62 N02-N08
      4  505 1332 (S01E38 L=245)   222     95        0 ( 464/1311- 557/1354)  -8.1 E41-E34 S01-N00
      5 1789 1223 (N25W46 L=329)   867    244        3 (1758/1108-1815/1339) -66.8 W42-W49 N18-N32
      6  188 1260 (S10E65 L=218)   474     58        1 ( 165/1231- 212/1285) -23.2 E69-E61 S11-S08
      7  795 1226 (S01E18 L=266)   312     92        2 ( 747/1200- 830/1249)   0.0 E20-E15 S02-N00
      8 1708 1149 (N18W39 L=322)   264     44        0 (1688/1134-1735/1166)   0.0 W37-W41 N16-N19
      9 1850  827 (N05W61 L=344)   770    159        1 (1771/ 797-1904/ 870)  31.6 W53-W66 N01-N09
     10  965  765 (S23W03 L=286)   190     45        0 ( 942/ 742- 988/ 785) -18.4 W01-W04 S23-S21
     11  490  573 (S47E28 L=255)   821    148        2 ( 431/ 532- 565/ 614)  -8.1 E34-E20 S48-S44
     12  616  544 (S46E16 L=267)   334     70        1 ( 589/ 518- 643/ 578) -56.3 E18-E13 S47-S43
  

    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