Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1643 (N38E46 L=328)   426     59        1 ( 551/1606- 587/1670) -39.8 E49-E42 N35-N40
      2  729 1621 (N38E30 L=344)   948    157        0 ( 657/1586- 801/1647)  -4.1 E36-E22 N36-N40
      3  310 1524 (N27E68 L=305)   582     49        1 ( 283/1508- 345/1536)   5.7 E74-E63 N25-N28
      4 1298 1396 (N27W16 L= 30)  1259    342        1 (1220/1296-1352/1504) -56.3 W10-W19 N20-N34
      5  924 1465 (N29E11 L=  3)   234     37        1 ( 900/1450- 949/1487)   0.0 E12-E08 N28-N30
      6  467 1362 (N18E42 L=331)   524    166        1 ( 411/1291- 509/1401) -26.6 E48-E36 N14-N20
      7 1742 1266 (N22W52 L= 66)   424    106        0 (1718/1212-1758/1309) -80.5 W50-W53 N18-N24
      8 1590 1270 (N21W38 L= 51)   324     77        1 (1568/1240-1623/1302) -63.4 W35-W40 N19-N23
      9  938 1122 (N07E06 L=  8)   837    259        3 ( 855/1045-1044/1187)  45.0 E10-W00 N01-N12
     10 1300  562 (S25W23 L= 37)   394     94        0 (1254/ 545-1347/ 583)  -8.1 W19-W26 S26-S24
  

    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