Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1627 1540 (N15W56 L=128)   506    111        0 (1591/1485-1651/1584)  71.6 W54-W58 N12-N18
      2 1155 1366 (N21W17 L= 89)  1851    456        1 ( 957/1239-1320/1442) -56.3 W07-W23 N09-N30
      3 1715 1171 (S06W46 L=117)   914    196        1 (1615/1142-1786/1199) -11.3 W37-W52 S08-S04
      4 1692  929 (S18W38 L=109)  1602    292        0 (1648/ 793-1736/1058)  47.3 W32-W44 S25-S12
      5  849  937 (N06E12 L= 59)   251     83        0 ( 816/ 899- 884/ 971) -78.7 E13-E10 N03-N08
      6 1775  891 (S24W45 L=116)   450     51        0 (1757/ 860-1791/ 928)  63.4 W42-W46 S25-S21
      7  841  645 (S10E21 L= 51)   290     25        1 ( 818/ 620- 866/ 665)   0.0 E22-E18 S11-S08
      8  958  583 (S17E16 L= 55)   258     47        0 ( 932/ 561- 984/ 604)  18.4 E17-E14 S17-S15
      9  745  474 (S18E33 L= 39)  1285    321        1 ( 637/ 363- 810/ 600) -90.0 E35-E30 S26-S07
     10  258  537 (S04E71 L=  0)  1012    105        1 ( 241/ 485- 281/ 590)  20.6 E80-E63 S07-S01
     11  512  470 (S13E49 L= 22)  1057    173        2 ( 430/ 437- 576/ 495) -20.0 E57-E44 S15-S10
     12  841  188 (S41E47 L= 24)  1449    150        2 ( 804/ 119- 870/ 243)  39.1 E59-E40 S47-S34
  

    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