Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1353 1844 (N55W66 L=219)   612     61        2 (1329/1809-1366/1879)  11.8 W58-W82 N51-N56
      2  601 1437 (N37E26 L=127)   539     72        3 ( 570/1399- 636/1464)  31.0 E28-E23 N34-N39
      3 1564 1367 (N18W42 L=196)   287    102        0 (1513/1354-1607/1382)  -8.1 W38-W46 N16-N19
      4  677 1170 (N19E21 L=133)   359    154        2 ( 625/1107- 720/1203) -54.5 E24-E17 N14-N21
      5  611 1093 (N15E26 L=128)   264    136        0 ( 586/1038- 653/1147) -66.8 E27-E23 N11-N18
      6  450  926 (N06E39 L=115)   921    204        5 ( 400/ 855- 493/1029)  48.0 E43-E34 N02-N12
      7 1704  652 (S29W46 L=199)   815    113        1 (1661/ 613-1757/ 689)  11.3 W41-W51 S31-S26
      8  308  575 (S15E62 L= 91)   610     93        2 ( 288/ 541- 341/ 616)  78.7 E64-E58 S17-S12
      9  401  460 (S24E58 L= 95)  1367    205        3 ( 370/ 359- 438/ 548)  25.9 E84-E50 S33-S17
     10  630  366 (S32E42 L=112)   738    137        1 ( 593/ 312- 660/ 437)  32.7 E49-E35 S36-S27
  

    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