Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1003 1829 (N59W38 L= 50)   460     44        1 ( 972/1814-1032/1849) -18.4 W33-W42 N57-N60
      2  950 1785 (N58W27 L= 39)   477     57        1 ( 928/1763- 978/1812)  63.4 W23-W31 N54-N60
      3 1571 1657 (N27W64 L= 76)  1084    117        2 (1521/1624-1630/1685) -45.0 W58-W67 N23-N29
      4 1440 1614 (N29W48 L= 61)   348     57        1 (1408/1590-1469/1636) -63.4 W46-W49 N26-N30
      5  625 1603 (N54E15 L=357)   290     43        0 ( 610/1577- 647/1623)  36.9 E16-E12 N52-N55
      6 1502 1497 (N21W45 L= 58)  1171    250        0 (1409/1460-1604/1534) -33.7 W38-W52 N15-N24
      7 1281 1331 (N18W23 L= 36)   444     64        0 (1239/1309-1321/1353) -45.0 W21-W25 N15-N20
      8 1097 1229 (N17W09 L= 22)   415     96        3 (1065/1195-1132/1276) -90.0 W07-W11 N14-N20
      9  542 1198 (N28E27 L=345)   534    127        1 ( 486/1165- 597/1224) -18.4 E31-E22 N25-N30
     10 1349 1044 (N00W19 L= 31)  1292    291        3 (1291/ 930-1408/1181)  42.5 W13-W26 S05-N06
     11  599 1024 (N17E26 L=346)   319    115        1 ( 576/ 978- 642/1063) -71.6 E27-E23 N13-N19
     12 1493  781 (S19W22 L= 35)   358    101        0 (1446/ 757-1539/ 798) -33.7 W19-W25 S21-S16
     13 1391  589 (S28W12 L= 24)   720    199        0 (1312/ 551-1476/ 650)  -4.8 W06-W18 S30-S26
     14 1654  592 (S37W34 L= 46)   726     72        2 (1619/ 555-1686/ 619)  45.0 W30-W37 S40-S35
     15 1580  508 (S41W26 L= 38)   558     75        0 (1541/ 491-1619/ 526) -14.0 W22-W30 S43-S39
     16 1459  426 (S43W14 L= 26)   784    117        2 (1401/ 386-1514/ 462)   0.0 W09-W19 S45-S40
     17  789  284 (S34E41 L=332)   716     92        1 ( 761/ 250- 826/ 324)  74.1 E42-E37 S37-S30
     18  494  287 (S28E70 L=302)  1993     88        0 ( 452/ 260- 529/ 311) -20.6 E78-E64 S30-S25
  

    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