Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1302 1504 (N24W31 L=214)  3949    868        9 ( 966/1447-1664/1535) -23.1 W10-W57 N13-N33
      2  782 1338 (N29E06 L=177)   435    110        3 ( 751/1296- 822/1386) -36.9 E07-E04 N25-N31
      3 1248 1237 (N11W18 L=201)   240     70        0 (1231/1201-1262/1282)  53.1 W16-W20 N08-N12
      4 1905 1125 (S15W65 L=248)  5583    622       11 (1831/ 997-1959/1267)  86.8 W59-W70 S23-S05
      5  285 1216 (N34E46 L=137)   420    109        1 ( 243/1185- 317/1232)  -8.1 E50-E43 N31-N35
      6 1814 1135 (S11W54 L=238)   607    101        1 (1781/1103-1841/1172)   8.1 W51-W58 S13-S09
      7 1596  745 (S27W28 L=211)  2748    376        3 (1413/ 533-1706/1037)  36.3 W11-W42 S36-S13
      8 1864  996 (S21W56 L=239)   569     72        0 (1849/ 958-1881/1028)  53.1 W54-W58 S23-S18
      9  477  817 (N06E38 L=145)   323     62        2 ( 452/ 790- 500/ 851)  53.1 E40-E36 N04-N08
     10  637  754 (S01E29 L=154)   337    111        0 ( 589/ 703- 675/ 795) -71.6 E31-E27 S04-N02
     11 1428  702 (S24W15 L=199)   642    185        1 (1358/ 606-1492/ 791)  31.6 W08-W21 S28-S20
     12  697  640 (S09E29 L=154)   209     46        0 ( 666/ 625- 725/ 657) -33.7 E30-E27 S10-S07
     13 1114  543 (S25E08 L=175)   656    233        2 (1064/ 442-1153/ 642)  50.7 E14-E03 S30-S19
     14  756  423 (S23E34 L=149)  1616    507        3 ( 663/ 282- 848/ 599)  31.6 E48-E22 S31-S14
     15  490  401 (S19E54 L=129)   459    106        0 ( 466/ 364- 516/ 449)  80.5 E56-E52 S21-S15
     16  883  259 (S38E36 L=148)   626    119        1 ( 831/ 228- 922/ 285)   5.7 E41-E31 S40-S36
     17  755  242 (S36E46 L=138)   444    134        1 ( 715/ 212- 781/ 270)  20.6 E50-E42 S37-S33
  

    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