Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  784 1708 (N54W06 L=350)  2572    442        3 ( 607/1604- 952/1805)   5.6 E14-W27 N50-N56
      2 1487 1731 (N31W61 L= 44)   805    106        2 (1442/1682-1536/1764) -81.9 W58-W63 N26-N34
      3 1487 1646 (N26W53 L= 36)   333     59        1 (1456/1632-1518/1663) -45.0 W50-W54 N24-N28
      4 1425 1616 (N27W46 L= 29)   266     82        0 (1403/1572-1442/1646)  59.0 W43-W47 N23-N28
      5  717 1474 (N39E07 L=336)   301     47        2 ( 688/1453- 746/1495) -36.9 E09-E05 N37-N40
      6 1520 1452 (N15W43 L= 26)   631    140        2 (1462/1417-1584/1488) -54.5 W39-W46 N11-N18
      7  512 1139 (N23E29 L=315)  1181    303        3 ( 376/1078- 645/1192) -26.6 E39-E19 N17-N29
      8  270  967 (N19E52 L=291)  1789    446        5 ( 215/ 758- 347/1137) -87.5 E57-E48 N05-N28
      9  350 1037 (N21E43 L=300)   483    154        0 ( 294/ 992- 401/1097) -45.0 E48-E39 N17-N25
     10 1562 1010 (S10W31 L= 15)   751    159        0 (1518/ 942-1598/1091)  40.6 W27-W35 S13-S06
     11 1246  960 (S04W11 L=354)   276     86        1 (1221/ 921-1271/1008)  38.7 W08-W13 S06-S02
     12  412  933 (N14E40 L=303)   277     59        2 ( 383/ 911- 444/ 965)  11.3 E42-E37 N12-N15
     13 1666  818 (S25W35 L= 19)  1499    270        7 (1557/ 739-1736/ 916)  13.2 W25-W42 S27-S19
     14 1269  811 (S13W08 L=352)   586    151        2 (1228/ 746-1317/ 874) -83.7 W06-W10 S18-S09
     15 1461  822 (S18W20 L=  4)   308     96        0 (1430/ 787-1487/ 864)  31.0 W17-W22 S19-S16
     16  859  416 (S27E28 L=316)   511    112        0 ( 817/ 387- 908/ 452)  12.5 E32-E23 S28-S25
     17  643  283 (S30E53 L=291)  1524    263        2 ( 573/ 214- 668/ 388)  29.7 E65-E43 S35-S23
  

    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