Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  898 1424 (N30E07 L=121)   628    223        2 ( 832/1346- 936/1511) -57.5 E11-E04 N25-N36
      2  377 1383 (N29E49 L= 79)   311    127        1 ( 343/1351- 402/1410)  36.9 E51-E45 N26-N30
      3 1275 1337 (N22W19 L=146)   796    238        0 (1199/1276-1360/1371)   5.2 W12-W23 N18-N24
      4  256 1299 (N23E60 L= 67)   950    110        2 ( 233/1244- 285/1351)  90.0 E62-E56 N19-N26
      5 1830 1273 (N13W65 L=192)   735     89        2 (1799/1242-1863/1311) -53.1 W60-W69 N11-N15
      6  186 1232 (N19E68 L= 60)  1158    117        2 ( 154/1175- 211/1281)  60.3 E71-E63 N15-N22
      7 1573 1196 (N11W38 L=165)   281     54        1 (1554/1169-1604/1225)  63.4 W35-W40 N09-N13
      8 1849  951 (S06W61 L=188)   500     71        2 (1836/ 921-1866/ 990) -68.2 W59-W63 S08-S03
      9 1527  838 (S10W32 L=159)   328     78        2 (1510/ 794-1547/ 878)  78.7 W30-W33 S12-S07
     10 1647  818 (S12W41 L=168)   258     37        0 (1632/ 792-1661/ 845)  76.0 W40-W42 S14-S10
     11  979  744 (S12E04 L=123)   430    201        0 ( 932/ 666-1040/ 815)  52.1 E07-W00 S17-S08
     12 1529  745 (S16W32 L=159)   250    101        0 (1497/ 712-1569/ 764)  18.4 W29-W35 S17-S14
     13  506  731 (S12E36 L= 91)   252     81        2 ( 482/ 687- 526/ 758)  53.1 E37-E34 S14-S10
     14 1811  692 (S23W61 L=188)  1145    127        2 (1762/ 661-1858/ 738)  18.4 W55-W67 S24-S20
     15 1530  647 (S23W33 L=160)   550    174        1 (1486/ 593-1571/ 698)  90.0 W29-W36 S26-S19
     16  617  534 (S25E32 L= 96)  1028    166        4 ( 594/ 457- 655/ 600)  51.3 E35-E27 S30-S20
     17  781  505 (S28E20 L=107)   242     62        0 ( 758/ 468- 803/ 531) -68.2 E21-E19 S30-S25
     18  512  431 (S33E45 L= 82)  1197    202        2 ( 422/ 411- 602/ 455)   0.0 E54-E35 S34-S30
  

    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