Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  702 1770 (N31E45 L=316)   428     56        2 ( 669/1743- 729/1798)   0.0 E48-E41 N29-N33
      2  810 1622 (N23E29 L=332)   251     64        0 ( 794/1591- 826/1661) -63.4 E29-E26 N21-N25
      3  687 1596 (N19E35 L=325)   470    134        0 ( 625/1575- 756/1614)  12.5 E39-E30 N17-N19
      4 1639 1562 (N43W32 L= 32)   404     36        1 (1611/1547-1662/1585)  18.4 W28-W34 N41-N45
      5 1416 1438 (N26W14 L= 14)   353    117        0 (1386/1396-1463/1472)  59.0 W12-W16 N23-N29
      6  826 1334 (N06E18 L=342)   861    383        1 ( 638/1281- 967/1396)  35.2 E28-E11 N00-N12
      7 1471 1272 (N18W20 L= 21)   283    140        0 (1440/1245-1529/1320)  -8.1 W17-W24 N16-N19
      8  199 1179 (S14E62 L=299)   488    121        0 ( 160/1139- 220/1230)  50.2 E65-E59 S16-S10
      9 1719 1062 (N13W42 L= 42)   533    131        1 (1705/ 983-1732/1126) -76.0 W40-W44 N08-N16
     10 1343 1009 (S01W18 L= 19)   284     96        1 (1301/ 989-1403/1036)   0.0 W15-W22 S01-N00
     11  145 1012 (S24E67 L=293)   453     64        0 ( 119/ 987- 170/1032) -26.6 E71-E62 S25-S22
     12 1303  727 (S17W24 L= 24)   221     45        0 (1274/ 705-1324/ 753)  56.3 W22-W24 S18-S15
     13  998  724 (S25W06 L=  6)   249     60        0 ( 972/ 706-1037/ 747)  26.6 W04-W08 S25-S23
     14 1355  598 (S23W32 L= 32)  2050    580        2 (1121/ 533-1622/ 696)  39.2 W18-W46 S31-S09
     15 1700  542 (S15W60 L= 61)   835    104        1 (1644/ 525-1752/ 570)  20.0 W55-W67 S17-S12
     16  949  351 (S50W17 L= 17)   499    100        1 ( 914/ 326- 994/ 390)   0.0 W13-W21 S51-S46
     17 1321  262 (S43W52 L= 53)   769    112        2 (1292/ 209-1349/ 311) -14.0 W45-W61 S45-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