Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  485 1538 (N51E34 L=246)   326     70        0 ( 466/1507- 515/1561)  36.9 E36-E30 N48-N52
      2 1241 1511 (N30W27 L=307)  2634    736        4 ( 930/1412-1548/1547) -28.9 W06-W44 N16-N39
      3 1853 1193 (S05W63 L=343)  1821    172        4 (1802/1129-1907/1261) -56.3 W57-W71 S09-N00
      4  318 1093 (N24E48 L=232)  1357    261        1 ( 279/1006- 357/1209) -77.0 E51-E44 N18-N31
      5  819 1083 (N15E11 L=269)   871    233        3 ( 788/ 955- 850/1175)  81.9 E12-E08 N06-N20
      6 1657 1076 (S05W41 L=321)   282     81        0 (1633/1044-1680/1100) -76.0 W38-W42 S07-S03
      7 1772  792 (S26W47 L=327)  1383    240        4 (1696/ 686-1824/ 885)  31.6 W39-W53 S31-S21
      8  424  708 (S00E46 L=234)   987    247        0 ( 388/ 612- 519/ 763)  19.7 E51-E37 S05-N02
      9  353  531 (S11E60 L=220)  7149    490       14 ( 253/ 385- 454/ 663)  61.2 E73-E51 S21-S01
     10 1299  401 (S39W03 L=283)   344     48        0 (1273/ 372-1326/ 429)  31.0 W00-W05 S40-S37
  

    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