Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1697 (N34E16 L=307)  1127    185        4 ( 983/1638-1148/1740)  -7.6 E23-E08 N32-N36
      2  330 1544 (N09E62 L=261)   655    123        2 ( 293/1484- 356/1584)  90.0 E65-E59 N05-N12
      3  330 1470 (N05E58 L=265)   352     56        1 ( 311/1438- 351/1500)  90.0 E59-E56 N02-N06
      4 1324 1421 (N23W08 L=331)   252     69        0 (1308/1382-1338/1460) -78.7 W06-W09 N20-N25
      5 1178 1410 (N18E00 L=322)   321     87        0 (1140/1387-1227/1431)   0.0 E03-W02 N16-N19
      6  284 1391 (S01E59 L=264)   909    120        3 ( 233/1370- 354/1418)   0.0 E66-E51 S02-N00
      7 1850 1108 (N20W54 L= 17)  2697    357        7 (1818/1005-1934/1261) -40.2 W49-W66 N14-N30
      8 1910 1186 (N28W62 L= 25)   666     85        1 (1876/1141-1930/1216) -26.6 W56-W65 N25-N29
      9 1636 1125 (N14W34 L=357)   402    108        1 (1619/1070-1662/1173) -59.0 W32-W36 N11-N16
     10 1371 1076 (N04W18 L=341)   509    189        1 (1283/1047-1455/1112)  37.9 W13-W22 N00-N07
     11  489  896 (S26E31 L=292)   250     48        0 ( 461/ 873- 512/ 915) -14.0 E33-E28 S27-S24
     12 1797  768 (S00W57 L= 20)   606     95        0 (1757/ 754-1852/ 790)   4.8 W52-W64 S01-N01
     13 1717  640 (S09W56 L= 19)   587    116        3 (1690/ 603-1750/ 690) -20.6 W51-W61 S10-S07
     14  291  645 (S46E52 L=271)  2005    287        4 ( 170/ 604- 429/ 683)   0.0 E76-E35 S47-S43
     15 1698  549 (S14W60 L= 23)   960    121        3 (1669/ 487-1728/ 605) -54.5 W56-W62 S18-S11
     16 1634  508 (S19W57 L= 20)   330     66        1 (1617/ 472-1651/ 538) -45.0 W53-W60 S20-S16
     17 1458  491 (S25W44 L=  7)   759    110        2 (1402/ 465-1509/ 516)  11.3 W38-W48 S26-S22
     18 1469  337 (S33W57 L= 20)  2007    275        2 (1340/ 307-1593/ 378)  32.9 W46-W66 S38-S26
  

    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