Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1170 1484 (N32W17 L=320)  1249    293        3 (1033/1432-1304/1526) -25.2 W07-W25 N27-N36
      2 1705 1469 (N22W60 L=  3)   468    104        0 (1659/1445-1748/1492) -29.7 W56-W64 N20-N24
      3  613 1290 (N26E25 L=278)   284    106        0 ( 564/1278- 667/1302)  -7.1 E28-E20 N25-N27
      4  723 1168 (N18E18 L=285)   217     81        1 ( 682/1147- 755/1192)  21.8 E20-E15 N16-N18
      5  397 1115 (N17E42 L=261)   255     39        0 ( 371/1098- 420/1130)  14.0 E44-E39 N15-N17
      6 1510 1000 (S02W30 L=333)   215     80        1 (1475/ 966-1537/1025) -63.4 W28-W32 S04-N00
      7 1779  801 (S19W52 L=355)  4065    530        9 (1675/ 614-1876/ 985)  52.7 W43-W64 S30-S09
      8 1382  832 (S10W20 L=323)   928    191        4 (1331/ 744-1425/ 895)  60.9 W16-W23 S15-S06
      9 1540  850 (S12W31 L=334)   374    104        2 (1500/ 803-1575/ 880) -68.2 W28-W33 S14-S09
     10  177  753 (S05E71 L=231)  4725    273        8 ( 133/ 626- 207/ 878)  59.0 E88-E65 S13-N02
     11 1475  719 (S19W26 L=329)   664    216        2 (1391/ 666-1563/ 764)  13.0 W19-W32 S21-S16
     12  855  694 (S13E14 L=289)   373     94        0 ( 811/ 673- 900/ 716) -38.7 E16-E11 S14-S10
     13  811  533 (S23E20 L=283)  1526    362        2 ( 747/ 452- 932/ 673) -63.4 E23-E12 S29-S13
     14  614  441 (S28E38 L=265)   665    125        1 ( 590/ 385- 641/ 498)  45.0 E42-E33 S31-S23
     15 1598  383 (S48W45 L=348)  1215     92        2 (1548/ 350-1633/ 420)  12.5 W38-W51 S51-S44
     16 1053  311 (S43E09 L=294)   600     65        0 (1021/ 288-1090/ 333)   9.5 E11-E05 S45-S41
  

    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