Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1598 1511 (N15W52 L=313)   355    121        0 (1550/1474-1640/1533) -45.0 W48-W54 N11-N17
      2 1654 1352 (N05W49 L=310)  1172    246        3 (1582/1315-1739/1381) -22.6 W43-W56 N02-N08
      3 1237 1282 (N14W19 L=280)   672    232        4 (1147/1211-1339/1317)   0.0 W11-W26 N12-N14
      4  338 1021 (N22E44 L=217)  3138    698       12 ( 147/ 964- 549/1121) -28.0 E61-E29 N13-N31
      5  849 1020 (N10E10 L=251)   274     72        1 ( 831/ 985- 879/1064) -90.0 E11-E08 N07-N12
      6 1707  836 (S25W38 L=299)  2701    467        3 (1586/ 651-1789/1025)  39.0 W25-W48 S33-S16
      7  926  676 (S12E15 L=246)   691    257        8 ( 860/ 617- 989/ 733)  15.3 E20-E09 S15-S08
      8  966  554 (S20E17 L=244)   392    110        3 ( 915/ 524-1016/ 580)   7.1 E20-E12 S21-S18
      9  493  431 (S16E53 L=208)  1523    406        5 ( 451/ 295- 537/ 555)  54.2 E62-E48 S25-S07
     10  634  361 (S24E46 L=215)  1966    324        6 ( 582/ 243- 693/ 473)  66.0 E51-E41 S33-S15
  

    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