Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  581 1765 (N30E55 L=312)  1647    176        1 ( 522/1701- 625/1822) -20.6 E71-E46 N25-N34
      2  822 1576 (N21E26 L=342)   729    245        1 ( 697/1555- 930/1606)  10.6 E34-E18 N18-N24
      3 1299 1533 (N30W05 L= 13)   414    178        0 (1230/1486-1381/1564)  -5.2 W00-W11 N28-N30
      4 1585 1320 (N24W29 L= 36)  1566    384        1 (1486/1187-1704/1455) -28.8 W19-W39 N18-N30
      5  161 1031 (S22E66 L=302)  1374    287        3 ( 105/ 908- 207/1154)  46.8 E78-E60 S29-S13
      6  247  997 (S23E55 L=313)   314     48        0 ( 223/ 977- 267/1024)  31.0 E57-E52 S24-S21
      7 1376  995 (S01W21 L= 28)   212     58        0 (1351/ 970-1404/1014)  45.0 W19-W22 S02-N00
      8  350  900 (S27E43 L=324)  1025    256        0 ( 265/ 826- 409/ 994)  45.0 E52-E38 S33-S20
      9  760  793 (S26E11 L=357)   209     84        1 ( 742/ 760- 778/ 832) -78.7 E12-E09 S28-S23
     10  467  744 (S35E32 L=335)  1178    289        2 ( 359/ 656- 580/ 808)  35.0 E43-E23 S42-S28
  

    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