Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1383 1631 (N44W31 L=130)  1265    212        2 (1328/1564-1422/1683)  33.7 W24-W36 N38-N48
      2  909 1472 (N32E09 L= 91)   384    107        0 ( 892/1421- 926/1520)  90.0 E09-E07 N28-N35
      3 1715 1459 (N30W55 L=154)   461    100        2 (1693/1434-1746/1480) -26.6 W52-W58 N28-N31
      4  550 1306 (N20E33 L= 66)   270     78        0 ( 538/1262- 563/1352)  90.0 E33-E31 N17-N23
      5 1537 1312 (N21W35 L=135)   488    106        1 (1484/1273-1588/1349) -26.6 W31-W39 N18-N23
      6 1538 1226 (N16W34 L=134)   555    173        0 (1467/1205-1600/1258)   0.0 W29-W39 N14-N17
      7 1077 1165 (N12W03 L=102)   379     66        1 (1055/1114-1103/1224) -80.5 W01-W04 N09-N15
      8  228 1190 (N12E60 L= 39)   385     76        1 ( 206/1171- 250/1219) -71.6 E62-E57 N10-N13
      9 1541 1096 (N08W34 L=133)   391    159        0 (1479/1076-1611/1130) -16.7 W29-W39 N06-N09
     10  475  997 (N01E36 L= 63)   410    130        0 ( 422/ 968- 529/1034)  20.6 E40-E32 N00-N03
     11  765  744 (S14E16 L= 83)   264     65        1 ( 746/ 718- 792/ 774)  71.6 E17-E14 S15-S12
     12  376  710 (S18E46 L= 53)   445    172        2 ( 359/ 663- 390/ 762)  63.4 E48-E44 S20-S14
     13 1562  629 (S22W38 L=137)   695     90        2 (1521/ 593-1598/ 665)  21.8 W35-W41 S24-S19
     14  691  629 (S22E22 L= 77)   621    124        1 ( 630/ 605- 754/ 662) -20.6 E26-E17 S23-S19
     15  496  607 (S24E38 L= 62)   783    244        3 ( 415/ 579- 577/ 655)  19.7 E44-E30 S26-S20
     16  968  352 (S43E04 L= 95)  1192    164        4 ( 899/ 319-1034/ 413)   5.2 E10-W01 S45-S37
     17  700  377 (S41E27 L= 73)   399     83        0 ( 670/ 357- 722/ 409)  38.7 E29-E24 S42-S38
  

    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