Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  759 1806 (N35E47 L= 19)  2061    471        3 ( 727/1720- 791/1917) -38.7 E58-E38 N28-N45
      2  777 1682 (N26E37 L= 29)   306     34        0 ( 756/1663- 802/1703) -14.0 E38-E34 N25-N27
      3  775 1532 (N17E30 L= 36)   264     46        1 ( 749/1509- 810/1547)  45.0 E31-E27 N15-N18
      4 1695 1421 (N40W35 L=100)  2209    359        4 (1567/1293-1769/1497) -25.2 W22-W42 N32-N44
      5  810 1261 (N02E18 L= 47)   240    112        0 ( 753/1251- 875/1272)  23.2 E21-E14 N00-N04
      6 1394 1131 (N11W18 L= 83)   328     97        1 (1374/1073-1416/1180) -63.4 W15-W19 N07-N13
      7 1494 1135 (N14W24 L= 89)   539    221        0 (1455/1063-1551/1180) -81.9 W21-W28 N09-N17
      8  265 1070 (S21E52 L= 14)   577    154        2 ( 242/1018- 280/1125) -56.3 E53-E49 S24-S18
      9  469  930 (S25E32 L= 33)   426    180        0 ( 403/ 910- 547/ 958)  31.0 E36-E26 S27-S21
     10  669  804 (S27E15 L= 50)   222     29        0 ( 644/ 787- 688/ 827)  56.3 E16-E13 S28-S25
     11  453  694 (S40E30 L= 35)   774    206        2 ( 371/ 669- 530/ 728)  23.2 E37-E23 S44-S37
     12 1347  298 (S36W50 L=115)   423     78        2 (1310/ 284-1389/ 314)  31.0 W47-W52 S37-S34
  

    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