Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  789 1844 (N37E48 L=343)  1495    191        2 ( 740/1763- 842/1921) -48.8 E58-E42 N29-N45
      2  553 1688 (N20E53 L=339)   759     98        3 ( 512/1665- 593/1721)  78.7 E56-E49 N17-N22
      3  804 1654 (N24E33 L=358)   279    103        0 ( 766/1618- 844/1685) -14.0 E36-E28 N22-N24
      4  770 1600 (N19E32 L=359)   739    183        3 ( 685/1539- 824/1635)  -9.5 E38-E26 N16-N21
      5  518 1477 (N06E43 L=348)  1475    364        3 ( 348/1427- 639/1561)  35.2 E56-E37 N01-N13
      6 1031 1418 (N15E10 L= 21)   218     94        0 (1006/1399-1073/1452) -11.3 E12-E07 N14-N16
      7 1150 1309 (N12E01 L= 30)   286    146        0 (1086/1296-1218/1321)  23.2 E04-W03 N10-N14
      8 1744 1245 (N27W40 L= 71)   897    152        1 (1694/1183-1793/1282)  -6.3 W35-W45 N25-N29
      9 1383 1214 (N13W15 L= 46)   332     92        0 (1362/1169-1410/1259) -38.7 W12-W17 N11-N15
     10  238 1196 (S15E58 L=333)  1149    140        2 ( 189/1145- 293/1245)  31.0 E62-E52 S18-S11
     11 1088 1188 (N04E01 L= 30)   454    107        1 (1030/1163-1138/1211)  26.6 E03-W02 N02-N05
     12  209 1019 (S26E58 L=333)   958    191        0 ( 186/ 932- 251/1102) -45.0 E61-E52 S30-S21
     13 1328  891 (S06W20 L= 52)   485    175        0 (1299/ 819-1349/ 989) -66.0 W17-W22 S09-N00
     14  437  884 (S29E34 L=357)  1865    382        6 ( 325/ 832- 648/ 934)  15.1 E44-E18 S33-S22
     15  823  850 (S22E07 L= 24)   405    136        0 ( 754/ 834- 883/ 870)  20.6 E11-E03 S22-S19
     16 1050  774 (S20W08 L= 40)   333     98        2 (1015/ 736-1093/ 815) -31.0 W06-W11 S21-S17
     17 1241  709 (S18W21 L= 53)   457     93        1 (1199/ 683-1269/ 755)  78.7 W19-W23 S19-S14
  

    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