Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1002 1794 (N45E30 L=320)   994    122        2 ( 948/1760-1061/1820)   5.2 E35-E23 N42-N47
      2  586 1557 (N17E43 L=306)   623    204        1 ( 495/1525- 676/1589)   7.1 E51-E35 N15-N17
      3  437 1505 (N11E53 L=297)   553    127        0 ( 404/1450- 472/1558)  90.0 E53-E51 N06-N14
      4  691 1463 (N14E32 L=318)   238     64        1 ( 661/1432- 714/1489) -26.6 E34-E29 N12-N15
      5 1546 1127 (N16W29 L= 19)   354    151        0 (1492/1092-1600/1184) -12.5 W24-W33 N14-N17
      6  748  815 (S22E11 L=338)   381    212        0 ( 692/ 747- 807/ 897) -31.0 E16-E06 S25-S18
      7  987  801 (S17W03 L=353)   217     66        3 ( 959/ 769-1013/ 831) -36.9 W01-W05 S18-S15
      8  919  778 (S20W00 L=350)   313     94        1 ( 902/ 735- 938/ 819)  90.0 E01-W01 S22-S17
      9 1586  760 (S04W42 L= 31)   261    100        2 (1551/ 725-1633/ 779)  51.3 W39-W44 S06-S01
     10 1134  564 (S27W20 L= 10)   609    186        2 (1036/ 548-1214/ 588)  31.0 W13-W24 S30-S24
  

    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