Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  774 1810 (N63W11 L= 41)   517     36        1 ( 753/1780- 793/1836)  29.7 W07-W15 N59-N65
      2  611 1738 (N63E12 L= 18)  2001    256        6 ( 502/1667- 724/1799)   0.0 E28-W04 N58-N65
      3  474 1442 (N43E29 L=  2)   309     76        0 ( 437/1428- 516/1454) -14.0 E32-E24 N40-N44
      4 1509 1295 (N07W35 L= 65)   509     87        1 (1483/1256-1536/1334) -78.7 W33-W37 N04-N09
      5  686 1256 (N25E14 L= 16)   566    127        1 ( 636/1216- 737/1296)   7.1 E18-E10 N23-N27
      6 1827 1265 (S04W60 L= 90)   371     47        1 (1803/1241-1845/1285) -71.6 W57-W62 S05-S02
      7 1828 1118 (S12W55 L= 85)   713    153        0 (1798/1048-1857/1191)  31.0 W50-W60 S14-S08
      8 1413 1101 (S01W24 L= 54)   240     14        0 (1391/1078-1431/1124)  45.0 W22-W25 S02-N00
      9  936  654 (S14E14 L= 16)   373    107        1 ( 908/ 607- 959/ 697)  63.4 E15-E11 S17-S11
     10  601  482 (S17E41 L=349)  1047    231        2 ( 548/ 363- 641/ 586)  32.9 E51-E34 S23-S10
     11 1639  513 (S44W31 L= 61)   624    117        1 (1577/ 498-1698/ 533) -12.1 W24-W38 S46-S41
     12 1370  419 (S42W06 L= 36)  1214    264        2 (1257/ 367-1486/ 447)  -9.0 E03-W16 S44-S38
  

    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