Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  554 1698 (N24E54 L=222)  2588    370        4 ( 441/1597- 678/1779)  70.6 E58-E48 N15-N32
      2  706 1723 (N30E44 L=232)   317     43        0 ( 685/1702- 723/1749) -76.0 E45-E42 N27-N31
      3 1145 1496 (N26E05 L=271)  1877    453        3 (1034/1288-1198/1636) -48.2 E16-W01 N14-N33
      4 1853 1194 (N29W56 L=332)  1641    240        1 (1798/1086-1896/1292) -39.8 W49-W62 N23-N33
      5  351 1141 (S13E46 L=230)   388    107        0 ( 336/1089- 368/1198) -56.3 E47-E43 S15-S09
      6  173 1070 (S21E64 L=212)  1195    175        0 ( 132/ 998- 247/1123)  34.7 E70-E55 S25-S16
      7  485  968 (S20E32 L=244)   399    116        0 ( 466/ 928- 519/1027) -51.3 E33-E28 S21-S16
      8  200  955 (S27E59 L=217)   938    118        1 ( 145/ 927- 267/ 995)  20.6 E67-E51 S29-S23
      9 1658  880 (N05W42 L=318)  1016    212        1 (1589/ 810-1706/ 977)  74.7 W39-W44 N00-N11
     10  507  505 (S50E25 L=251)  1726    255        2 ( 410/ 431- 615/ 591) -15.3 E35-E13 S53-S45
     11 1685  467 (S16W66 L=342)  1746    262        1 (1574/ 409-1787/ 528)  40.9 W58-W76 S22-S09
     12 1058  226 (S52W38 L=314)   327     40        0 (1040/ 202-1075/ 254) -20.6 W33-W41 S54-S50
  

    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