Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  833 1800 (N40E35 L=328)   803    104        1 ( 787/1771- 885/1824) -15.3 E40-E29 N37-N41
      2 1073 1721 (N39E13 L=350)   644    105        1 (1014/1681-1117/1754) -11.3 E18-E08 N35-N40
      3 1179 1700 (N40E05 L=358)  1095    138        2 (1129/1636-1214/1754) -41.6 E09-W00 N35-N43
      4  343 1603 (N18E66 L=297)   835     94        0 ( 314/1555- 368/1641) -74.1 E68-E62 N14-N21
      5 1422 1614 (N39W16 L= 18)   325     78        1 (1387/1589-1452/1635)  45.0 W12-W18 N36-N41
      6  495 1413 (N08E41 L=322)   369     90        0 ( 457/1379- 530/1449)  59.0 E43-E39 N05-N10
      7 1066 1371 (N14E04 L=359)   412    218        1 ( 952/1344-1153/1384)  10.3 E10-W01 N12-N15
      8 1534 1343 (N23W26 L= 29)   247     53        0 (1508/1319-1559/1369) -26.6 W24-W28 N21-N23
      9  210 1136 (S12E60 L=303)  2490    340        6 ( 109/1023- 276/1261)  36.4 E74-E53 S18-S04
     10  307  983 (S20E48 L=314)   408    133        1 ( 293/ 924- 317/1054)  90.0 E49-E47 S23-S15
     11  270  779 (S33E54 L=309)   745    139        1 ( 227/ 730- 302/ 849)  36.0 E60-E49 S36-S28
     12 1618  790 (S08W42 L= 45)   290     74        2 (1598/ 751-1639/ 821)  90.0 W40-W43 S10-S05
     13 1693  678 (S12W53 L= 55)   646    130        2 (1657/ 628-1732/ 739) -20.0 W47-W58 S13-S09
     14 1299  624 (S24W26 L= 29)   869    259        1 (1271/ 509-1341/ 732) -58.4 W23-W31 S30-S17
     15 1541  488 (S26W50 L= 53)  1851    347        4 (1425/ 401-1675/ 585)  59.5 W44-W58 S34-S17
     16 1590  392 (S30W62 L= 65)  1106    165        1 (1514/ 351-1655/ 437)  52.1 W57-W66 S34-S25
  

    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