Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1075 1353 (N23W12 L=109)   287    168        0 (1057/1281-1111/1427)  77.5 W10-W15 N17-N26
      2  757 1211 (N23E10 L= 86)  1766    585        4 ( 507/1044- 941/1333) -51.8 E25-E02 N09-N37
      3 1937 1226 (S12W76 L=173)  1718     74        1 (1910/1186-1956/1256) -39.8 W70-W86 S14-S09
      4 1676  961 (S16W37 L=134)  2216    762        3 (1474/ 706-1800/1244)  29.1 W17-W54 S25-S04
      5 1869  955 (S24W56 L=153)  1297    170        4 (1802/ 884-1917/1010)  18.4 W47-W63 S25-S20
      6 1357  751 (S18W11 L=108)   809    282        1 (1304/ 621-1408/ 862)  40.2 W04-W17 S24-S13
      7 1087  757 (S10E04 L= 93)   371    161        1 (1062/ 669-1121/ 828)  48.8 E07-W00 S14-S06
      8  373  607 (S02E53 L= 43)   309     57        1 ( 346/ 587- 396/ 629)  11.3 E56-E51 S02-N00
      9  616  506 (S13E39 L= 57)   667    207        4 ( 594/ 415- 632/ 614)  53.1 E44-E34 S18-S06
     10  395  469 (S10E59 L= 38)  4129    481        5 ( 262/ 323- 525/ 580) -73.3 E72-E50 S22-S01
     11  527  535 (S09E44 L= 52)   594     81        4 ( 493/ 499- 549/ 578)  33.7 E47-E41 S10-S06
  

    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