Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  888 1896 (N52E40 L=292)   861     99        1 ( 848/1873- 934/1916)  -4.4 E46-E33 N49-N53
      2  640 1775 (N35E48 L=283)   415     52        2 ( 619/1754- 657/1803) -33.7 E51-E45 N33-N38
      3 1210 1452 (N22W04 L=335)   301    113        1 (1171/1415-1254/1481) -23.2 W00-W07 N20-N23
      4 1852 1091 (N15W57 L= 28)   673     84        2 (1817/1058-1877/1132)  51.3 W53-W60 N12-N17
      5  415  946 (S20E39 L=293)   479    167        1 ( 366/ 885- 459/ 998)  53.1 E42-E35 S24-S16
      6 1705  653 (S14W55 L= 26)   312     55        2 (1674/ 635-1732/ 673)  21.8 W52-W57 S15-S12
      7 1696  564 (S19W59 L= 30)  1101    198        3 (1617/ 508-1775/ 619)  33.7 W51-W64 S22-S13
      8  330  585 (S45E53 L=279)   627     41        1 ( 310/ 554- 351/ 617)  33.7 E56-E48 S46-S42
      9 1132  548 (S33W17 L=349)   321    114        1 (1071/ 530-1182/ 573)  23.2 W13-W20 S34-S30
     10 1629  478 (S25W58 L= 30)   608     92        0 (1597/ 447-1676/ 504)  51.3 W55-W61 S27-S22
     11 1034  248 (S57W26 L=358)   419     44        1 (1002/ 230-1061/ 268)   9.5 W23-W29 S58-S54
  

    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