Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1084 1638 (N39W25 L=161)   800    221        1 (1045/1557-1142/1704)  30.3 W19-W32 N33-N41
      2 1493 1607 (N23W49 L=185)  1236    212        4 (1426/1564-1578/1635) -56.3 W43-W53 N18-N27
      3  555 1255 (N30E23 L=113)  1075    407        0 ( 394/1225- 733/1297) -23.7 E35-E10 N26-N37
      4 1557 1176 (S01W35 L=171)  1344    274        6 (1452/1132-1655/1216) -15.9 W27-W41 S05-N02
      5  815 1082 (N14E10 L=126)   267     72        1 ( 784/1058- 853/1117) -53.1 E11-E08 N12-N16
      6 1713  996 (S16W41 L=177)   595    115        2 (1666/ 966-1752/1031)  -7.1 W37-W45 S17-S14
      7 1443  848 (S16W19 L=155)   418    145        1 (1428/ 793-1468/ 940)  77.5 W17-W21 S19-S10
      8 1106  705 (S14E04 L=132)   715    220        0 ( 986/ 668-1199/ 729) -42.0 E10-W00 S19-S10
      9  815  396 (S26E32 L=104)   292    142        2 ( 785/ 341- 857/ 450)  77.5 E33-E30 S30-S21
     10 1053  227 (S46E28 L=108)  1081    136        3 (1021/ 182-1081/ 272)  31.0 E33-E23 S49-S41
  

    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