Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1282 1802 (N58E07 L=217)   814     74        2 (1260/1756-1303/1844) -70.0 E10-E03 N52-N63
      2 1009 1642 (N34E19 L=206)   678    214        1 ( 966/1549-1065/1725) -35.5 E25-E11 N28-N39
      3 1534 1549 (N44W23 L=248)   726     76        1 (1499/1504-1568/1585) -38.7 W19-W26 N40-N46
      4 1134 1480 (N26E05 L=219)   246     36        1 (1116/1456-1158/1505) -56.3 E06-E03 N24-N28
      5 1537 1338 (N28W25 L=250)   373    131        0 (1513/1279-1565/1398) -45.0 W22-W28 N25-N31
      6 1831 1185 (N28W55 L=280)  1978    386        2 (1720/1027-1929/1307) -20.1 W42-W73 N20-N33
      7 1194  791 (S12W16 L=240)   211     49        1 (1176/ 761-1216/ 814) -76.0 W14-W16 S14-S10
      8  209  770 (S37E62 L=163)   855     74        1 ( 181/ 738- 239/ 801)  24.4 E67-E56 S39-S34
      9 1428  606 (S17W36 L=261)   242    142        0 (1356/ 589-1473/ 640)  33.7 W32-W38 S19-S13
     10  271  575 (S51E61 L=163)   679     74        1 ( 244/ 548- 297/ 607)  17.5 E72-E53 S53-S47
     11 1324  445 (S29W38 L=262)  2590    497        4 (1099/ 379-1528/ 514)  37.2 W24-W49 S38-S19
     12  422  428 (S59E43 L=182)   434     27        0 ( 406/ 405- 438/ 455)  32.0 E47-E38 S61-S56
  

    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