Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  869 1766 (N34E36 L=264)   651    104        2 ( 840/1724- 916/1815) -69.4 E37-E33 N30-N38
      2 1288 1612 (N36E01 L=298)   707    146        0 (1239/1562-1330/1675)  80.5 E02-W00 N30-N42
      3 1166 1584 (N30E08 L=292)   484    129        1 (1101/1554-1224/1604)   6.3 E12-E03 N27-N30
      4 1201 1382 (N18W00 L=300)   352    140        0 (1154/1346-1243/1437)  82.9 E01-W01 N14-N22
      5 1554 1257 (N21W25 L=325)   949    236        0 (1477/1203-1640/1290)  35.0 W21-W32 N15-N23
      6 1895  940 (N15W64 L=  4)  2231    217        2 (1868/ 850-1937/1011)  79.7 W60-W71 N09-N20
      7 1075  916 (S12W06 L=305)   318    121        0 (1032/ 879-1138/ 946)   0.0 W02-W10 S12-S10
      8 1016  728 (S24W08 L=308)   869    301        1 ( 911/ 671-1118/ 816)  58.4 W03-W11 S29-S16
      9 1218  685 (S20W21 L=321)   363    114        1 (1170/ 652-1261/ 711)   9.5 W17-W23 S21-S19
     10 1553  622 (S13W45 L=344)   574    155        0 (1491/ 583-1628/ 653)  48.8 W41-W49 S17-S09
  

    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