Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1391 1727 (N53W02 L= 79)   435     61        1 (1361/1697-1417/1754)  76.0 W00-W04 N48-N56
      2 1107 1479 (N24E08 L= 69)   220    110        0 (1059/1459-1141/1507)  53.1 E11-E06 N21-N26
      3 1697 1402 (N39W37 L=114)   420     52        0 (1677/1377-1719/1426) -21.8 W34-W39 N36-N40
      4 1254 1313 (N18W05 L= 82)   243    121        0 (1216/1259-1291/1343) -26.6 W02-W08 N15-N18
      5 1477 1148 (N15W23 L=100)   222     48        2 (1461/1124-1496/1179) -45.0 W21-W24 N13-N16
      6  211 1085 (S21E60 L= 17)   641    100        2 ( 190/1035- 241/1137) -81.9 E62-E56 S24-S17
      7 1593 1078 (N14W32 L=109)   275     57        0 (1575/1052-1609/1104) -33.7 W30-W33 N12-N15
      8 1514  376 (S26W59 L=136)   467    120        1 (1459/ 363-1568/ 396)  45.0 W54-W62 S28-S22
      9  988  250 (S53W33 L=110)   711     92        0 ( 958/ 219-1026/ 292) -18.4 W27-W39 S55-S50
     10 1117  209 (S51W48 L=125)   455     58        1 (1093/ 181-1143/ 239)  -9.5 W42-W54 S52-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