Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  799 1763 (N49E35 L=293)   454     68        0 ( 763/1745- 839/1782)  33.7 E37-E31 N46-N50
      2  363 1462 (N20E57 L=271)  1055    197        0 ( 311/1395- 444/1503)  45.0 E61-E49 N15-N23
      3 1474 1368 (N28W28 L=356)  1153    228        0 (1450/1273-1524/1455)  77.0 W25-W32 N21-N34
      4  325 1389 (N15E57 L=271)   530     95        1 ( 289/1351- 360/1434)  51.3 E60-E54 N12-N18
      5  288 1317 (N11E58 L=270)   694    115        1 ( 244/1284- 334/1341)  12.5 E62-E52 N07-N12
      6 1284 1198 (N14W14 L=343)  1431    301        5 (1163/1147-1402/1248)   3.6 W06-W22 N10-N18
      7 1009  691 (S20W03 L=331)   905    343        1 ( 942/ 551-1066/ 838) -66.0 E02-W08 S28-S10
      8  318  726 (S26E50 L=278)   324     93        0 ( 307/ 686- 332/ 775)  71.6 E51-E47 S28-S22
      9  382  731 (S25E43 L=285)   368     81        1 ( 368/ 697- 396/ 774) -90.0 E44-E42 S26-S21
     10 1650  528 (S23W55 L= 23)  2291    372        2 (1547/ 435-1752/ 662)  60.6 W49-W60 S30-S14
     11 1055  411 (S39W11 L=339)   529     68        1 (1022/ 379-1092/ 440)  31.0 W08-W13 S41-S36
     12 1022  247 (S54W14 L=342)  1126    161        0 ( 978/ 183-1070/ 305)  90.0 W10-W17 S62-S47
  

    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