Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1208 1870 (N59W54 L= 49)  1344    173        2 (1139/1819-1284/1913)   0.0 W39-W80 N56-N60
      2  768 1508 (N41E10 L=345)   573    166        1 ( 743/1417- 797/1577)  70.0 E12-E07 N34-N45
      3 1755 1398 (N14W61 L= 57)   431     97        1 (1727/1378-1791/1420) -36.9 W57-W64 N12-N15
      4  473 1310 (N31E35 L=320)  2092    412        4 ( 317/1208- 644/1415) -25.9 E54-E21 N22-N39
      5  293 1014 (N15E51 L=304)   607    156        2 ( 250/ 958- 351/1051)  39.8 E55-E45 N11-N16
      6  183  889 (N08E65 L=290)  1242    215        0 ( 143/ 807- 215/ 994) -57.5 E70-E62 N02-N14
      7 1904  958 (S15W67 L= 62)  1178    113        1 (1857/ 932-1950/ 992)   0.0 W59-W77 S16-S12
      8  943  733 (S10E09 L=346)   741     64        4 ( 914/ 672- 975/ 805)  63.4 E11-E06 S13-S05
      9 1161  696 (S15W03 L=359)   822    348        2 (1022/ 641-1285/ 764) -29.7 E04-W10 S19-S09
     10  482  565 (S14E44 L=311)   827    259        1 ( 421/ 513- 549/ 622) -36.9 E48-E38 S18-S10
     11  407  540 (S16E52 L=303)   535     69        0 ( 368/ 512- 435/ 575) -45.0 E54-E50 S17-S12
     12  657  334 (S33E42 L=314)  1620    321        3 ( 543/ 267- 763/ 421)  17.4 E60-E28 S37-S27
  

    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