Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  604 1790 (N54E40 L= 38)   539     75        1 ( 573/1768- 634/1805) -15.3 E46-E35 N51-N55
      2 1684 1558 (N29W56 L=134)  1675    226        1 (1591/1489-1767/1628) -39.8 W49-W63 N24-N34
      3 1020 1509 (N28W02 L= 80)   368    157        0 ( 983/1455-1092/1544)  40.6 W00-W07 N24-N30
      4 1486 1436 (N21W33 L=112)   283     65        0 (1448/1416-1513/1459)  21.8 W30-W35 N20-N22
      5 1399 1387 (N18W26 L=104)   245     62        2 (1382/1359-1423/1412) -56.3 W24-W27 N16-N19
      6  180 1247 (N16E64 L= 15)   533     84        1 ( 160/1208- 203/1279)  45.0 E66-E61 N13-N18
      7 1454  923 (S10W26 L=105)   569    177        2 (1392/ 857-1518/ 994)  45.0 W22-W30 S13-S05
      8 1407  690 (S24W24 L=102)   424    138        1 (1388/ 621-1426/ 755) -82.9 W22-W25 S28-S20
      9 1717  650 (S28W51 L=130)   316     43        1 (1692/ 626-1732/ 674) -71.6 W48-W53 S29-S26
     10  475  609 (S25E42 L= 37)   339     97        1 ( 439/ 568- 504/ 648)  23.2 E45-E38 S27-S22
     11 1744  488 (S39W65 L=143)  1029    131        1 (1712/ 433-1771/ 558) -56.3 W61-W70 S42-S33
     12  662  292 (S49E41 L= 37)   807     94        3 ( 621/ 263- 706/ 328) -51.3 E45-E36 S51-S45
  

    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