Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1194 1590 (N36W28 L=293)  1445    362        1 (1027/1534-1338/1619) -14.0 W14-W38 N32-N39
      2 1522 1599 (N27W53 L=317)   845    136        1 (1449/1576-1589/1619) -17.1 W46-W59 N24-N29
      3 1643 1528 (N19W58 L=323)   554     69        2 (1608/1494-1679/1560) -68.2 W55-W60 N16-N21
      4 1304 1322 (N18W24 L=289)   221     83        0 (1260/1305-1348/1336)  -9.5 W20-W27 N16-N18
      5  203 1307 (N39E66 L=198)   515     48        0 ( 182/1289- 232/1323)   7.1 E70-E61 N38-N40
      6 1444 1273 (N11W31 L=296)   260     81        0 (1409/1240-1479/1297) -59.0 W29-W32 N08-N13
      7 1053 1200 (N17W06 L=270)   448    155        1 (1017/1119-1075/1249)  69.4 W03-W07 N11-N19
      8  521 1070 (N20E31 L=234)   306     74        1 ( 496/1049- 553/1100) -53.1 E32-E28 N17-N21
      9  519  911 (N10E34 L=231)   377    169        0 ( 486/ 839- 561/ 984) -78.7 E34-E31 N05-N15
     10 1879  831 (S28W61 L=326)  2171    217        4 (1813/ 746-1933/ 914)  19.4 W51-W72 S31-S23
     11 1573  793 (S20W29 L=294)   244     56        0 (1548/ 769-1599/ 817)  21.8 W26-W31 S20-S18
     12  637  752 (S01E29 L=236)   981    324        0 ( 570/ 636- 747/ 816)  18.4 E36-E21 S06-N02
     13  466  610 (S07E46 L=219)  1811    251        6 ( 401/ 527- 566/ 687) -32.5 E52-E37 S11-S01
     14  301  519 (S11E68 L=197)  1847    136        4 ( 250/ 466- 351/ 563) -36.9 E75-E61 S14-S07
     15  587  513 (S15E40 L=224)   327     66        0 ( 564/ 482- 608/ 552) -90.0 E41-E39 S16-S11
     16  515  473 (S16E48 L=217)  1226    202        6 ( 474/ 387- 560/ 546)  42.3 E56-E41 S21-S11
     17 1441  422 (S42W15 L=279)  1211    217        2 (1357/ 357-1526/ 503)  17.4 W06-W22 S45-S36
     18  863  437 (S25E25 L=240)   248     85        0 ( 830/ 409- 897/ 460)  15.9 E28-E21 S26-S23
     19  594  405 (S22E45 L=219)   564     96        1 ( 553/ 364- 618/ 438) -90.0 E47-E42 S25-S19
  

    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