Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1118 1829 (N59W35 L= 57)  1404    258        2 (1006/1777-1222/1869)   1.7 W18-W52 N56-N60
      2 1452 1596 (N33W44 L= 65)   471    146        0 (1383/1558-1512/1615) -36.9 W39-W47 N29-N35
      3 1370 1494 (N29W33 L= 54)   528    135        1 (1312/1473-1444/1517) -29.1 W29-W38 N26-N31
      4  468 1343 (N33E36 L=345)  1887    374        5 ( 412/1170- 522/1501)  76.6 E39-E32 N22-N43
      5  262 1241 (N28E56 L=325)  2668    295        6 ( 167/1138- 348/1348) -23.4 E77-E45 N21-N35
      6  727 1215 (N22E16 L=  5)   295     81        0 ( 689/1199- 768/1235)   0.0 E18-E13 N20-N23
      7 1578 1056 (S01W35 L= 56)   603    108        0 (1518/1035-1629/1084)   0.0 W30-W39 S02-N00
      8 1706  865 (S15W43 L= 64)  1296    193        3 (1625/ 816-1775/ 923)   8.7 W36-W50 S17-S12
      9 1587  849 (S13W33 L= 54)   243     28        0 (1567/ 827-1604/ 872) -56.3 W31-W34 S14-S11
     10  877  625 (S16E15 L=  6)  1138    336        2 ( 754/ 550- 964/ 742)  35.5 E23-E08 S19-S09
     11  558  640 (S11E36 L=345)  1190    190        2 ( 539/ 539- 592/ 720)  63.4 E40-E32 S17-S05
     12  709  647 (S12E25 L=356)   277     68        0 ( 672/ 633- 748/ 662) -26.6 E27-E23 S13-S11
     13 1651  474 (S42W43 L= 64)   496     72        1 (1614/ 451-1685/ 496)   8.1 W38-W46 S43-S40
     14 1503  367 (S47W29 L= 50)   877     91        2 (1462/ 345-1547/ 395)   6.3 W24-W34 S49-S44
     15 1360  307 (S49W15 L= 36)   754     93        2 (1334/ 261-1400/ 342) -77.5 W11-W19 S54-S45
     16 1284  254 (S53W07 L= 28)  1263    131        2 (1237/ 223-1335/ 294) -35.0 W01-W12 S57-S48
     17 1283  190 (S61W07 L= 28)   692     62        1 (1241/ 176-1316/ 206) -26.6 W01-W11 S64-S58
  

    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