Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1367 1592 (N31W32 L= 58)   422    149        0 (1298/1552-1428/1620) -32.0 W27-W35 N27-N32
      2 1720 1454 (N18W56 L= 83)   344     38        0 (1697/1434-1738/1475)  21.8 W53-W58 N16-N19
      3  892 1403 (N24E04 L= 22)   464    166        0 ( 820/1347- 956/1433) -36.9 E08-W00 N19-N26
      4 1780  773 (S24W53 L= 79)   357    102        0 (1752/ 730-1797/ 814)  59.0 W50-W54 S26-S21
      5 1725  648 (S31W49 L= 76)   336     33        2 (1704/ 627-1743/ 674) -56.3 W47-W51 S32-S29
      6  292  499 (S24E68 L=318)  1091    141        3 ( 257/ 452- 316/ 556)  54.5 E76-E62 S27-S20
      7 1563  377 (S50W43 L= 69)   519     60        2 (1534/ 355-1594/ 394)   0.0 W37-W47 S52-S48
      8 1381  331 (S51W23 L= 49)   909     88        1 (1331/ 297-1426/ 381) -22.6 W17-W29 S55-S47
      9 1475  331 (S53W34 L= 61)   779     79        2 (1436/ 305-1515/ 361) -23.2 W27-W41 S56-S49
  

    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