Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1338 1815 (N57E08 L=271)   695     72        3 (1295/1794-1366/1843)  26.6 E12-E05 N53-N62
      2  394 1462 (N02E53 L=226)   324     45        0 ( 367/1445- 415/1484)  33.7 E55-E51 N01-N04
      3 1237 1399 (N20W01 L=281)   651    305        0 (1104/1313-1294/1444)   0.0 E05-W06 N16-N23
      4  274 1375 (S05E61 L=218)   694     95        1 ( 237/1345- 308/1402)  -5.7 E66-E56 S06-S02
      5 1833 1150 (N26W52 L=331)   454     69        1 (1803/1131-1859/1172)   0.0 W47-W55 N24-N28
      6 1898 1040 (N23W62 L=342)  2334    217        2 (1872/ 945-1929/1153) -72.9 W58-W68 N16-N29
      7  178 1064 (S25E64 L=215)   515     67        0 ( 156/1028- 202/1096)  45.0 E67-E60 S26-S22
      8 1792 1055 (N19W49 L=328)   352     57        0 (1775/1021-1810/1082)  78.7 W46-W50 N16-N21
      9 1773  968 (N13W49 L=328)   313     55        0 (1745/ 950-1802/ 985)  38.7 W46-W51 N11-N15
     10  362  926 (S29E41 L=238)  1482    338        5 ( 210/ 884- 529/ 960)   7.1 E58-E26 S31-S26
     11 1016  639 (S28W12 L=291)   262     83        0 ( 983/ 611-1049/ 672)  68.2 W10-W12 S30-S25
     12 1717  636 (S06W57 L=337)   328     72        0 (1700/ 609-1732/ 669) -21.8 W54-W60 S07-S04
     13 1606  523 (S15W55 L=334)   390     59        1 (1590/ 498-1633/ 553) -26.6 W51-W57 S16-S13
     14 1403  469 (S25W43 L=323)   689    189        0 (1320/ 452-1498/ 491)  38.7 W37-W48 S28-S20
     15 1611  394 (S21W66 L=346)  1912    145        4 (1563/ 342-1658/ 449) -63.4 W60-W75 S25-S17
     16 1560  419 (S22W59 L=338)   374     64        1 (1538/ 391-1578/ 448) -71.6 W56-W60 S23-S20
     17  609  301 (S68E07 L=272)   507     50        0 ( 585/ 270- 628/ 332) -45.0 E11-E03 S71-S65
  

    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