Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  314 1401 (S02E57 L=204)  2981    437        9 ( 164/1255- 428/1556)  54.5 E70-E53 S12-N09
      2 1517 1520 (N38W17 L=279)   532     99        0 (1475/1494-1567/1543)   0.0 W13-W21 N36-N39
      3  452 1247 (S08E40 L=221)   682    205        0 ( 427/1159- 506/1315)  84.3 E42-E38 S13-S03
      4 1136 1270 (N10E01 L=261)   246     71        0 (1116/1231-1155/1305) -45.0 E02-W01 N08-N11
      5  305 1030 (S24E48 L=214)  1016    194        0 ( 265/ 944- 379/1097)  54.0 E51-E41 S29-S18
      6 1619  801 (S01W41 L=303)   309     81        0 (1583/ 777-1651/ 837)  53.1 W39-W42 S02-N02
      7 1089  813 (S16W09 L=271)   199     39        0 (1060/ 796-1113/ 827)   0.0 W07-W11 S16-S14
      8 1213  772 (S14W18 L=279)   306     75        2 (1178/ 747-1246/ 795)  63.4 W16-W19 S16-S12
      9  446  588 (S48E31 L=231)   672    133        0 ( 389/ 569- 497/ 610)   8.7 E37-E24 S49-S45
     10 1444  553 (S19W40 L=302)   366     95        0 (1400/ 538-1497/ 565)  33.7 W37-W43 S20-S16
     11 1321  494 (S26W36 L=297)   512    122        1 (1258/ 454-1372/ 514)  12.5 W31-W40 S27-S24
     12 1019  323 (S47W28 L=290)  1335    323        1 ( 890/ 261-1167/ 386)   6.1 W15-W43 S48-S44
  

    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