Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1514 1522 (N29W45 L=222)   382     73        0 (1481/1493-1553/1543) -53.1 W42-W47 N26-N30
      2 1430 1485 (N28W36 L=213)   309     87        0 (1410/1448-1458/1524) -90.0 W35-W37 N25-N30
      3  366 1420 (N38E50 L=127)   651     97        2 ( 313/1399- 408/1441)  -8.7 E57-E44 N36-N39
      4  318 1074 (N16E49 L=128)   484    174        0 ( 295/1002- 357/1150) -63.4 E51-E45 N10-N20
      5  200  913 (N06E63 L=114)   986    126        2 ( 168/ 858- 238/ 963)  35.0 E69-E58 N02-N09
      6 1531  721 (S20W30 L=207)   243     63        1 (1507/ 694-1553/ 752)  45.0 W27-W31 S21-S18
      7 1475  359 (S48W30 L=207)   905    106        1 (1427/ 335-1523/ 393)  11.3 W24-W34 S49-S45
      8  532  293 (S39E60 L=117)   926    100        1 ( 497/ 262- 568/ 331)  15.6 E75-E50 S42-S35
  

    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