Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1481 1730 (N31W59 L=156)   452     39        0 (1456/1712-1511/1744)  -9.5 W55-W63 N29-N32
      2  427 1615 (N59E37 L= 60)   370     60        0 ( 412/1582- 448/1637)  68.2 E40-E34 N55-N60
      3 1682 1359 (N05W50 L=147)   871    237        0 (1624/1242-1719/1447)  80.5 W47-W52 S01-N11
      4  722 1195 (N22E13 L= 84)  1013    294        1 ( 589/1165- 863/1221) -25.2 E21-E04 N18-N26
      5 1621  925 (S16W33 L=130)   735    276        0 (1545/ 840-1716/1029)  17.4 W26-W42 S19-S13
      6 1653  787 (S25W33 L=130)  1749    519        1 (1543/ 590-1757/1026)  41.0 W21-W45 S34-S14
      7  196  915 (N17E60 L= 37)  2355    309        4 (  97/ 869- 334/ 995) -17.2 E76-E47 N12-N22
      8 1329  977 (S05W16 L=113)   183     78        0 (1296/ 955-1367/ 990) -14.0 W13-W17 S06-S04
      9  534  834 (N06E33 L= 64)   529    139        1 ( 510/ 780- 568/ 894)  56.3 E35-E30 N02-N09
     10  388  767 (N05E46 L= 51)   402     82        1 ( 352/ 740- 429/ 796) -51.3 E48-E43 N02-N07
     11 1821  745 (S34W50 L=147)   568     72        1 (1789/ 711-1850/ 778)  18.4 W46-W54 S35-S32
     12  806  703 (S08E20 L= 77)   374    104        0 ( 762/ 666- 856/ 739)  15.9 E23-E16 S08-S06
     13  391  604 (S04E51 L= 45)   856    142        1 ( 340/ 557- 440/ 677)  21.0 E57-E44 S06-N00
     14 1476  418 (S45W14 L=111)   661    112        1 (1422/ 396-1538/ 438) -35.0 W09-W19 S49-S41
     15 1157  316 (S42E14 L= 83)   413    114        0 (1100/ 288-1205/ 346)   0.0 E18-E09 S43-S40
     16  913  201 (S43E39 L= 58)  1549    286        3 ( 794/ 152-1056/ 251)  -4.1 E55-E25 S46-S40
  

    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