Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  848 1492 (N20E18 L=305)   403    119        0 ( 826/1431- 867/1549) -54.5 E20-E15 N16-N23
      2 1459 1304 (N17W23 L=346)   807    194        0 (1382/1236-1528/1375) -24.4 W17-W28 N13-N20
      3 1839 1275 (N23W57 L= 20)   808     97        0 (1817/1232-1869/1319) -45.0 W54-W60 N20-N25
      4 1783 1248 (N20W50 L= 13)   560    107        2 (1753/1203-1813/1300) -71.6 W46-W53 N17-N23
      5  358  952 (S17E44 L=279)   675    213        0 ( 289/ 913- 470/ 976)  17.4 E50-E34 S20-S15
      6  163  848 (S25E67 L=256)  1023    133        1 ( 123/ 792- 202/ 895)  19.4 E77-E60 S27-S21
      7  997  676 (S27W03 L=326)   722    303        0 ( 896/ 599-1069/ 807)  65.0 E02-W06 S33-S18
      8 1403  610 (S25W32 L=355)   286     81        0 (1358/ 583-1440/ 638)  51.3 W29-W33 S27-S22
      9 1488  409 (S35W49 L= 11)   397     52        0 (1458/ 386-1516/ 433)  63.4 W46-W50 S37-S33
     10  512  323 (S61E47 L=275)   728     49        1 ( 488/ 300- 545/ 346)  13.2 E56-E39 S62-S58
  

    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