Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1351 1737 (N47W43 L=131)  1825    313        6 (1297/1642-1410/1812)  55.3 W38-W51 N39-N52
      2  894 1483 (N37E03 L= 85)   734    209        1 ( 852/1398- 949/1594)  70.3 E06-W03 N30-N44
      3 1119 1499 (N34W14 L=102)   996    196        2 (1042/1419-1195/1563) -63.4 W09-W17 N27-N39
      4 1700 1252 (N09W50 L=138)  1469    356        3 (1575/1165-1817/1334) -37.4 W41-W60 N01-N16
      5 1279 1203 (N13W18 L=106)   288    103        0 (1226/1171-1324/1228)   8.1 W14-W21 N11-N14
      6  330 1082 (N16E48 L= 40)  1136    314        4 ( 253/1002- 406/1196) -42.7 E54-E41 N10-N23
      7  776 1116 (N14E14 L= 74)   304    123        0 ( 752/1054- 794/1181) -90.0 E15-E13 N10-N18
      8 1694  924 (S10W44 L=131)   267     68        0 (1668/ 881-1711/ 955)  53.1 W41-W45 S12-S08
      9 1556  576 (S31W32 L=120)  2275    386        6 (1385/ 489-1706/ 705)  19.1 W18-W44 S34-S24
     10  781  490 (S25E24 L= 64)   795    158        0 ( 714/ 458- 841/ 544)  20.0 E29-E18 S27-S21
  

    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